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ABSTRACT 


Estimates  of  probable  maximum  precipitation  (PMP)  and 
probable  maximum  flood  (PMF)  are  necessary  considerations  in 
the  design  and  operation  of  multipurpose  structures  in  which 
failure  can  have  catastrophic  consequences.  Climatological 
and  meteorological  studies  can  be  used  to  gain  an 
understanding  of  the  hydrological  phenomena  that  result  in 
major  rainstorms  causing  extreme  flooding. 

Meteorological  and  statistical  techniques  were  used  to 
obtain  estimates  of  PMP  for  the  six  major  river  basins  in 
Alberta.  Using  three  different  hydrological  models,  PMP 
estimates  were  employed  to  obtain  estimates  of  PMF.  Using 
the  statistical  technique  developed  by  Hershfield,  it  was 
found  that  point  estimates  of  PMP  were  well  correlated 
within  a  river  basin  and  were  dependent  on  the  position  and 
elevation  of  the  location  for  which  an  estimate  is  required. 
The  resulting  relationship  provides  an  easy  method  of 
computing  point  estimates  of  PMP  for  locations  lacking 
meteorological  data. 

The  rain-on-snow  event  had  the  best  potential  for 
producing  the  maximum  water  loading.  As  a  result,  the 
contribution  of  snowmelt  to  PMF  was  described,  and  an 
equation  was  developed  to  estimate  maximum  snowmelt  for 
Alberta  river  basins. 

To  examine  the  effect  of  basin  shape  on  PMF  the  author 


devised  the  Time-Area  Probable  Maximum  Flood  Model.  This 
numerical  model  uses  simulation  techniques  with  minimum  data 
requirements  and  computational  efficiency.  Its  main  features 
allow  subdivision  of  the  drainage  area,  accomodation  of  the 
shape  of  a  watershed,  consideration  of  flood  travel  time, 
and  production  of  flood  hydrographs  for  a  watershed. 

Two  well-known  models  were  also  examined  for  PMF 
estimates,  parameterization  and  sensitivity  in  the  Red  Deer 
River  Basin:  the  Streamflow  Synthesis  and  Reservoir 
Regulation  Model  and  the  HYMO  Model.  The  normalized  PMF 
discharges  from  the  three  models  were  compared  to  the 
maximum  recorded  discharges. 

The  spatial  distribution  and  Q-A  relationships  for  PMF 
estimates  are  also  presented  for  the  major  river  basins  in 
Alberta.  Combining  the  results  from  all  the  river  basins  in 
Alberta,  a  Q-A  relationship  with  exponent  equal  to  0.29 
seems  to  best  represent  the  data,  while,  for  individual 
river  basins  a  square-root  (exponent  equal  to  0.50) 
describes  best  the  relationship  for  the  enveloping  curve  of 
PMF.  For  many  river  basins  in  Alberta  the  coefficient  of  the 
Q-A  relationship  with  a  square-root  exponent  was  found  to 
vary  logarithmically  with  slope  and  intensity. 
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CHAPTER  1 


INTRODUCTION 


1.1.0  Background 

For  many  years  scientists  have  recognized  the 
importance  of  meteorological  and  climatological  phenomena  to 
hydrological  problems.  Major  floods  due  to  severe 
rainstorms,  for  example,  usually  are  a  result  of 
meteorological  conditions  and  hence  require  an  understanding 
of  the  climatology  of  a  given  region.  Rainfall  studies 
frequently  are  directed  toward  estimation  of  the  physical 
upper  limits  of  storm  rainfall  in  a  basin,  termed  the 
Probable  Maximum  Precipitation  (PMP).  The  American 
Meteorological  Society  (Huschke,  1970)  defines  the  PMP  as 
"the  theoretical  greatest  depth  of  precipitation  for  a  given 
duration  that  is  physically  possible  over  a  particular 
drainage  area  at  a  certain  time  of  year.  In  practice  this  is 
derived  over  flat  terrain  by  storm  transposition  and 
moisture  adjustment  to  observed  storm  patterns".  Another 
more  operational  definition  emphasizing  application 
(Huschke,  1970)  states  that  "PMP  is  that  magnitude  of 
rainfall  over  a  particular  basin  which  will  yield  the  flood 
flow  of  which  there  is  virtually  no  risk  of  being  exceeded." 

PMP  estimate  can  be  converted  into  flows  by  empirical 
methods  or  simulation  models  to  produce  a  theoretical  flood 
that  is  known  as  the  Probable  Maximum  Flood  (PMF).  Estimates 
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of  this  flood  incorporate  meteorological  and  hydrological 
parameters  in  an  attempt  to  describe  the  physical  upper 
limits  to  the  flooding  that  could  occur  over  a  certain  basin 
during  a  specified  time  interval. 

For  hydrological  applications,  estimates  of  PMP  and  PMF 
should  be  considered  in  the  design  and  operation  of 
multipurpose  structures.  In  some  cases,  the  lack  of  these 
estimates  has  caused  underestimation  of  the  necessary  design 
and  has  resulted  in  disastrous  consequences  when  a  rainfall 
or  flood  occurred  that  exceeded  the  one  used  in  the  design. 

Although  a  number  of  studies  of  PMP  have  been  conducted 
during  the  past  40  years  in  the  United  States,  only  a  few 
have  addressed  PMF,  and  these  were  only  for  specific  basins. 
In  Canada  the  number  of  studies  on  this  topic  is  even 
smaller.  Only  after  the  disastrous  hurricane  storm  of  14-15 
October  1954  (called  Hazel)  crossed  into  Ontario  from  the 
United  States  and  caused  extensive  property  damage  did 
hydrologists  in  Canada  initiate  the  much-needed  program  of 
PMP  studies  for  Canadian  river  basins.  Research  advances 
have  been  slow  in  the  last  20  years,  with  only  a  few  studies 
published  for  specific  basins.  In  Alberta,  there  are  very 
few  publications  on  PMP  and  PMF,  and  these  are  for  basins 
where  specific  projects  were  needed.  One  of  the  main  reasons 
for  limited  research  on  this  topic  is  that  these  estimates 
require  extensive,  time-consuming  studies  that  the  user  may 
not  have  the  resources  to  undertake.  Secondly,  the  required 
climatological  and  meteorological  data  are  time  consuming 
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and  expensive  for  the  potential  user  or  researcher  to 
obtain.  Thirdly,  the  user  often  is  not  familiar  with 
existing  hydrological  models,  especially  the  more 
sophisticated  ones.  The  user  must  Know  what  models  are 
available  and  where  to  obtain  them,  as  well  as  how  to  use 
these  models  and  recognize  their  limitations.  These 
limitations  are  of  Key  importance,  for  some  models  have  been 
misused  even  though  they  were  designed  for  specific  purposes 
and  regions.  A  final  reason  for  the  minimal  effort  on  this 
topic  has  been  the  lacK  of  a  simplified  approach  for 
estimating  PMP  and  PMF.  It  was  this  deficiency,  namely,  the 
1 acK  of  a  simplified  approach  for  estimating  PMF,  together 
with  an  interest  in  hydrological  research  into  severe  storms 
that  led  the  author  to  investigate  this  topic. 

1.1.1  Objectives  of  Study 

The  objectives  of  this  study  can  be  divided  into  two 
main  categories:  (a)  those  dealing  with  PMP  estimates,  and 
(b)  those  dealing  with  PMF  estimates. 

The  objectives  associated  with  PMP  estimates  are: 

(1)  to  obtain  a  better  understanding  of  the  spatial  and 
temporal  distribution  of  PMP  in  Alberta,  including 
climatological  and  meteorological  analyses  of  rainstorms  and 
estimates  of  PMP, 

(2)  to  obtain  PMP  estimates  for  the  six  main  river 
basins  in  Alberta  using  the  two  most  common  approaches:  (a) 
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the  meteorological  or  traditional  approach,  and  (b)  the 
statistical  technique,  and 

(3)  to  derive  a  model  for  the  determination  of  PMP  for 
any  location  in  Alberta. 

The  objectives  associated  with  PMF  estimates  are: 

(1)  to  examine  the  rain-on-snow  event  and  to  define  the 
relevant  parameters  and  processes  associated  with  snowmelt 
necessary  for  PMF  estimates, 

(2)  to  develop  a  conceptual  model  for  estimating  PMF  in 
Alberta  (in  this  thesis  termed  the  TAPMF  Model  1  )  and  to 
compare  the  results  of  this  model  with  estimates  obtained 
from  two  generally  used  hydrological  models  ( SSARR  2  and 
HYMO  3  )  4  i n  an  Alberta  watershed,  and 

(3)  to  examine  the  spatial  distribution  of  PMF 
estimates  in  Alberta  and  to  obtain  a  relationship  with  which 
to  estimate  PMF  for  ungauged  watersheds. 

The  PMP  objectives  of  this  work  have  been  researched 
previously  and  were  presented  in  two  publications  by  the 
author:  (1)  Verschuren  and  Wojtiw  (1980)  and  (2)  Wojtiw  and 
Verschuren  (1981).  A  summary  of  the  published  results  is 


1  TAPMF  Model: Time  Area  Probable  Maximum  Flood  Model. 

2  SSARR  Model : St reamf low  Synthesis  and  Reservoir  Regulation 
Model  (U.S.  Army  Corps  of  Engineers,  1972). 

3  HYMO : Hydrologica 1  Model  Computer  Language  (Williams  and 
Hann,  1973). 

4  The  SSARR  and  HYMO  models  are  poplular  models  used  by 
government  agencies  and  consultants  in  Alberta. 
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presented  in  Chapter  2  of  this  thesis.  The  PMF  objectives 
are  examined  in  Chapters  3  to  7. 

Before  the  results  of  the  research  work  are  presented, 
a  review  of  the  available  literature  is  given.  The  review 
has  been  divided  into  a  number  of  parts  to  give  the  reader 
better  understanding  of  the  history  and  techniques  involved 
in  PMP  and  PMF  studies:  (1)  climatological  studies;  (2) 
synoptic  studies;  (3)  PMP  studies;  (4)  historical  flood 
studies  in  Alberta;  (5)  snowmelt  studies;  (6)  flood 
estimates  and  PMF  studies;  and  (7)  computer  simulation 
mode  1 1 i ng  of  PMF . 

1.2.0  Literature  Review 

A  large  number  of  precipitation  studies  have  been 
conducted  in  Alberta.  Some  of  these  have  examined  the 
climatology  of  surface  rainfall,  while  others  attempted  to 
relate  the  synoptic  weather  condi t ions  to  precipitation  on 
the  surface. 

1.2.1  Climatological  Studies 

In  most  studies  of  the  climatology  of  surface 
precipitation,  point  measurements  are  examined  in  terms  of 
the  regional  and  temporal  variations  of  rainfall.  The 
results  of  the  point  measurement  analyses  indicate 
variations  in  the  point  measurements  of  precipitation  and 
not  variations  in  the  rainstorms,  which  are  areal  in 
dimension.  Comparisons  are  usually  made  of  the  heaviest 
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precipitation  between  one  point  and  another,  which  usually 
are  produced  by  two  independent  (in  time)  events.  A  good 
example  of  such  a  study  is  the  work  by  Storr  (1963,  1967), 
who  examined  the  precipitation  data  for  the  heaviest  1-,  2- 
and  3-day  rainfalls  from  81  stations  for  each  month  (May  to 
August)  from  1921  to  1960.  The  results  were  expressed  as 
maps  of  heaviest  rainfall  for  the  5-,  10- ,  and  25-year 
return  period.  These  maps  provide  the  spatial  distribution 
of  the  heaviest  rainfall,  with  each  point  representing  the 
maximum  observed  value  recorded  in  most  cases  from 
independent  events. 

The  first  study  to  examine  the  climatology  of 
rainstorms  in  Alberta  was  the  work  of  Verschuren  and  Wojtiw 
(1980)  in  which  the  authors  identified  and  analyzed  611 
rainstorms  (between  1921  and  1978)  that  produced  a  minimum 
depth  of  50  mm  (2  in.)  precipitation.  In  addition  to 
presenting  the  yearly  frequency  of  rainstorms  (Figure  1.1), 
the  authors  showed  that  an  average  of  about  11  storms 
ooccurred  each  year  with  depths  greater  than  50  mm  (2  in.) 
and  that  the  number  of  rainstorms  decreased  logarithmically 
with  increased  depth.  Over  50%  of  the  rainstorms  in  Alberta 
occurred  in  June  (23.6%)  and  July  (26.7%),  with  only  small 
percentages  occurring  in  April  (5.6%)  and  September  (10.1%) 
The  authors  found  that  the  greatest  frequency  was  in  the 
Waterton  Lakes  National  Park  area,  2:1  year  event  for  depth 
50  mm  (2  in.)  and  more  (Figure  1.2);  approximately  a  1:3 
year  event  for  depths  100  mm  (4  in.)  or  more;  and 


Number  of  Storms  with  Depth 
50mm  and  more 
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FIGURE  1.1.  Yearly  Frequency  of  Rainstorms  with  Maximum 
Depth  50  mm  (2  in.)  and  more  in  Alberta  (Verschuren  and 

Wojt iw ,  1 980  )  . 

approximately  a  1:10  year  event  for  depths  150  mm  (6  in.)  or 
more  (Figure  1.3).  The  frequency  decreased  along  the 
Continental  Divide,  with  a  number  of  pockets  of  maxima  in 
central  Alberta.  Severe  storms  5  occur  in  four  main  regions 
(or  belts)  of  the  province:  the  first  extends  through 
southern  Alberta  just  south  of  Calgary;  the  second  is  in 
central  Alberta  from  south  of  Edson  to  the  Edmonton  region; 
the  third  is  from  Lesser  Slave  Lake  to  the  Fort  McMurray 


5  Storms  that  produce  depths  150  mm  (6  in.)  and  more. 


8 


FIGURE  1.2. 
mm  (2  in.) 


Isopleths  of  the  Number  of  Storms  with  Depths  50 
and  more  in  100  years  (Verschuren  and  Wojtiw, 

1980) . 


K  I 


FIGURE  1.3. 
150  mm  (6  in 


Isopleths  of  the  Number  of  Storms  with  Depths 
)  and  more  in  100  years  (Verschuren  and  Woitiw 

1980) . 
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area;  and  the  fourth  is  around  the  Fort  Vermilion  area. 
Storms  with  depths  greater  than  150  mm  (6  in.)  usually  cause 
f 1 ood i ng . 

A  number  of  other  studies  have  also  been  published,  and 
although  they  are  listed  in  the  bibliography  they  are  not 
discussed  here  since  they  are  not  directly  relevant  to  the 
objectives  in  this  dissertation. 

1.2.2  Synoptic  Studies 

In  western  Canada  the  500  millibar  (mb)  cold  low  has 
been  recognized  by  forecasters  (Thompson,  1950)  as  a 
synoptic  feature  responsible  for  extensive  precipitation  on 
the  surface  in  Alberta.  Studies  of  the  synoptic  conditions 
of  severe  rainstorms  have  also  been  published  ( Mokosch , 

1961;  McKay,  1965b;  Burrows,  1966;  Warner  and  Thompson, 

1974;  and  Thompson,  1976).  Thompson's  1976  report  summarizes 
the  characteristics  of  storm  events  of  different  synoptic 
scale  to  identify  the  storm  types  most  likely  to  produce 
extreme  rainfall  events  over  the  Saskatchewan  Basin.  In 
Alberta  two  categories  seem  to  emerge  based  on  the  events' 
duration  and  synoptic  scale. 

The  first  category  contains  shor t -durat i on  and  small- 
scale  intense  rainfall  events;  these  are  produced  by 
thunderstorms.  The  typical  life  cycle  of  a  thunderstorm  is  3 
to  8  hours,  and  the  area  exposed  to  rainfall  is 
approximately  several  hundred  square  kilometres.  Most  of 
these  storms  result  from  daytime  heating  of  the  earth's 
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surface  so  that  although  land-use  patterns,  type  of 
vegetation,  and  aspect  may  play  roles  in  determining  the 
point  of  storm  formation,  topography  cannot  otherwise  be 
related  to  the  rainfall  intensity.  Many  of  the  1-day  extreme 
annual  rainfalls  at  stations  in  Alberta  are  of  this  type.  As 
a  rule,  these  usually  do  not  result  in  extensive  flooding  in 
Alberta . 

The  second  category  contains  medium-scale  events;  these 
are  associated  with  low-pressure  and  frontal  systems.  Most 
of  the  3-day  annual  extreme  rainfalls  and  probably  also 
several  of  the  1-day  annual  extreme  rainfalls  in  the  Eastern 
Slopes  and  plains  are  from  storms  of  this  scale.  Extreme  wet 
periods  longer  than  about  3  days  would,  in  most  instances, 
be  due  to  disturbances  on  a  planetary  scale  that  cause  a 
breakdown  in  the  normal  eastward  movement  of  the  large-scale 
pressure  configurations.  Under  such  conditions  the  low- 
pressure  generation  areas  and  the  upper- level  winds  that 
steer  the  storm  systems  remain  stationary  and  cause 
successive  medium-scale  storms  to  track  over  the  same 
general  area  for  an  extended  period.  Such  a  situation  has 
been  called  a  cold  low  and  is  classed  as  a  medium-scale 
storm.  The  500  mb  map  (roughly  at  5.5  km  or  18  000  ft.)  and 
the  surface  weather  map  are  used  to  identify  such  phenomena; 
examples  of  these  maps  are  shown  in  Figures  1.4  and  1.5, 
respectively. 

In  Figure  1.4,  pressure  isopleths  of  a  500  mb  low  are 
shown  over  south-central  Alberta,  there  is  cold  air  over 
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FIGURE  1.4.  The  500  mb  Analysis  of  a  Typical  Cold  Low  Over 

Alberta  (After  Thompson,  1976). 

southern  Alberta  and  British  Columbia,  and  warm  air  is 
spiralling  around  the  east,  north,  and  northwest  sides  of 
the  low.  In  Figure  1.5,  pressure  isopleths  of  about  1000  mb 
are  depicted  in  the  figure,  and  a  surface  low  shown 
southwest  of  Edmonton  is  causing  a  nor theaster ly  upslope 
flow  of  warm  moist  air  over  the  foothills  of  the  Peace  River 
Basin.  The  band  of  warm  air  on  the  east,  north,  and 
northwest  sides  of  the  upper  low  is  often  the  region  of 
heaviest  precipitation  in  a  cold  low  system,  even  in  the 
absence  of  an  upslope  condition.  The  existence  of  the 
appropriate  combination  of  topography  and  low-level  winds 
provides  additional  lift  to  the  already  ascending  air  to 
further  increase  the  intensity  of  the  precipitation. 

Thompson  (1976)  found  that  the  cold  low  was  the  most 
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FIGURE  1.5.  The  1000  mb  Analysis  of  a  Typical  Cold  Low  Over 

Alberta  (After  Thompson,  1976). 

frequent  producer  of  rainfalls  of  sufficient  intensity  to 
pose  a  flood  threat  in  the  South  Saskatchewan  River  Basin 
and  that  this  type  of  storm  occurs  most  frequently  during 
May  or  June. 

Burrows  (1966)  examined  daily  precipitation  data  for  16 
years  (1950  to  1965)  for  the  Edmonton  region  and  classified 
the  data  into  two  basic  types  of  storms:  cold  low  and  cold 
trough  storms.  In  meteorology  a  cold  trough  is  an  elongated 
area  of  relatively  low  atmospheric  pressure  (the  opposite  of 
a  ridge).  Of  the  60  storms  analyzed  (with  rainfall  19  mm 
(0.75  in.)  or  greater),  42  storms  were  identified  as  cold 
low,  while  18  were  classified  as  cold  trough.  In  their 
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study,  Verschuren  and  Wojtiw  (1980)  found  that  all  34  severe 
rainstorms  producing  150  mm  (6  in.)  or  more  in  precipitation 
at  the  surface  could  be  classified  according  to  Burrows' 
scheme  and  that  85%  of  the  storms  were  of  the  cold  low  type, 
with  the  remainder  being  of  the  cold  trough  type.  The 
authors  further  concluded  that  the  cold  low  system  is 
usually  associated  with  heavy  precipitation  in  Alberta. 

1.2.3  PMP  Studies 

In  the  United  States  a  large  number  of  PMP  studies  have 
been  conducted  over  the  past  40  years.  The  disastrous  floods 
during  the  late  1930s  directed  investigations  by  a  number  of 
researchers  (Bailey  and  Schneider,  1939;  Craeger ,  1939;  and 
Hathaway,  1939)  into  the  subject  of  PMF  and  its  -relationship 
to  spillway  capacity  in  considering  the  design  of  dams.  With 
these  studies  the  concept  of  PMP  gained  favor  and  a  number 
of  comprehensive  research  programs  were  initiated.  By  1939, 
through  advanced  studies  and  a  greater  accumulation  of  data, 
the  probability  method  had  been  proven  to  be  entirely 
inadequate.  Because  of  these  circumstances,  the  importance 
of  meteorological  studies  in  the  design  of  flood  control 
structures  was  realized  (Creager,  1939).  During  this  period 
the  idea  of  mass  rainfall  curves  was  developed  and  used  for 
storm  analysis  (Bailey  and  Schneider,  1939;  Hathaway,  1939). 
The  next  important  development  was  the  use  of  dew  points  as 
a  measure  of  the  precipi table  water  in  the  moist  air. 
Showalter  and  Solot  (1942)  concluded  that  the  maximum  amount 
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of  rain  that  can  occur  over  a  given  basin  can  be  determined 
from  the  amount  of  moist  air  that  can  flow  over  the  basin 
and  the  maximum  amount  of  moisture  that  can  be  precipitated 
from  that  moist  air. 

It  was  not  until  1953  (Paulhus  and  Gilman,  1953)  that 
Probable  Maximum  Precipitation  was  defined  and  a  distinction 
drawn  between  it  and  Probable  Maximum  Storm.  Until  then,  the 
terms  were  used  liberally  and  interchangeably  in  the 
literature.  The  distinction  is,  that  the  PMP  for  a  specific 
area  is  usually  determined  by  several  types  of  storms, 
whereas  the  Probable  Maximum  Storm  is  determined  by  one 
specific  storm.  Paulhus  and  Gilman  outlined  a  maximizing 
procedure  used  in  the  derivation  of  PMP  that  consists 
chiefly  of  moisture  adjustment  and,  when  justifiable, 
transposition  of  the  observed  storms.  This  procedure  serves 
today  as  the  basis  for  calculating  estimates  of  PMP  and  is 
found  in  the  literature  under  various  names:  traditional, 
meteorological,  or  physical. 

Moisture  adjustment  of  a  storm  involves  estimation  of 
the  increased  precipitation  that  could  be  expected  if 
maximum  atmospheric  moisture  were  available.  The  surface  dew 
point  is  used  as  an  index  of  the  moisture  in  the  storm  and 
of  the  maximum  moisture  in  the  basin.  Paulhus  and  Gilman 
further  mentioned  that  their  tests  indicated  that  the  lowest 
dew  point  in  the  12-hour  period  cor respondi ng  most  nearly  to 
the  12-hour  period  of  greatest  rainfall  was  most 
representative  of  the  average  moisture  in  the  storm  center. 
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This  dew  point  was  labelled  the  representative  12-hour 
persisting  dew  point.  The  maximum  12-hour  persisting  dew 
point  to  which  the  storm  values  are  adjusted  is  the  highest 
value  below  which  the  dew  point  did  not  drop  during  any  12- 
hour  recorded  period.  All  dew  points  were  reduced 
pseudoadiabatical ly  to  the  1000  mb  level  to  make  the  dew 
points  for  stations  at  different  elevations  comparable  and 
to  permit  the  use  of  tables  and  charts  of  moisture  in 
atmospheric  columns  with  base  at  1000  mb.  Maps  of  maximum 
dew  points  and  tables  of  moisture  content,  expressed  in 
terms  of  precipi table  water,  have  been  published  by  the  U.S. 
Department  of  Commerce  (1948,  1951). 

For  many  years  this  traditional  approach  continued  to 
be  popular  and  gained  general  acceptance.  In  1961  a 
different  way  of  calculating  PMP  was  introduced  by 
Hershfield  (1961),  who  proposed  a  statistical  method  for  the 
systematic  analysis  of  precipitation  data  in  estimating  PMP; 
this  method  was  modified  several  years  later  (Hershfield, 
1965).  He  proposed  that  the  PMP  at  an  observable 
precipitation  point  could  be  estimated  from  generalized 
frequency  equation  (Chow,  1951)  of  the  form 

XR  ■  1  +  K  Sn  (1.1) 

where  is  the  return  period  rainfall,  X  and  Sn  are  the 
mean  and  standard  deviation  of  a  series  of  n  annual  maxima 
for  a  specific  duration,  and  K  is  the  frequency  factor 


. 


17 


dependent  upon  the  frequency  distribution  chosen  to  fit  the 
data.  In  the  original  works  by  Hershfield  (1961,  1965)  the 
frequency  factor  was  found  to  vary  from  about  1.2  to  14.5 
and  was  postulated  to  have  a  maximum  value  of  15.  This  value 
of  15  has  been  used  by  many  researchers  to  estimate  PMP 
without  verification  or  question  of  applicability.  In  1977, 
Hershfield  (1977)  published  results  showing  that  the  PMP  at 
any  point  was  also  a  function  of  duration  and  the  mean 
annual  extreme  rainfall;  for  a  duration  of  i  hours  this  can 
be  expressed  as  follows; 


Ki  =  19  (10)‘b  Xi 


(1.2) 


where  K.  is  the  frequency  factor,  X.  is  the  mean  annual 
i  l 

extreme  rainfall  (inches)  for  the  i  hour  duration,  and  b  is 
a  constant  which  varies  with  duration.  For  durations  of  6 
and  24  hours,  Equation  1.2  can  be  expressed  as: 

(a)  for  6-hour  duration 


-0.00213 

19  (10) 


(1.3) 


(b)  for  24-hour  duration 


K 


24 


-0.000965 

19  (10) 


(  1.4) 


By  plotting  the  constant  part  of  the  empirical  power 
coefficients  against  duration  on  a  log-log  graph,  a  straight 
line  results,  thus  allowing  the  interpolation  or 
extrapolation  to  other  durations. 

This  method  is  now  known  as  the  statistical  method  of 
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estimating  the  PMP.  Similar  to  other  schemes  dependent  on 
empirical  coefficients,  the  equation  proposed  by  Hershfield 
is  a  concise  and  convenient  way  of  achieving  an  approximate 
answer  when  little  is  Known  about  the  estimated  quantity. 

Its  main  advantage  is  that  it  can  be  used  in  remote  areas 
where  little  meteorological  data  are  available  and  where  it 
is  virtually  impossible  to  make  a  reasonable  estimate  using 
the  physical  approach. 

The  PMP  obtained  by  the  physical  or  statistical  method 
is  an  enveloping  value.  Even  if  an  excellent  theoretical 
model  were  available,  the  analyst  would  still  have  to 
maximize  the  input.  At  the  present  time,  the  maximization  of 
the  parameters  is  imperfectly  understood. 

Comparisons  of  the  physical  and  statistical  methods 
have  been  made  in  a  number  of  studies.  Most  studies  in 
Canada  which  have  used  both  methods,  found  that  the  two 
methods  produced  similar  results.  During  the  late  1950s  and 
throughout  the  1960s  a  number  of  studies  were  conducted  to 
estimate  PMP  for  various  Canadian  river  basins.  One  of  the 
preliminary  studies  using  PMP  estimates  for  severe  storms 
was  researched  by  Bruce  (1957),  who  performed  an  analysis  of 
the  rainfall  for  the  most  severe  storm  on  record  in  Ontario 
i.e.,  hurricane  Hazel  (14-15  October  1954).  At  that  time, 
the  only  Depth- Area-Durat ion  (DAD)  analyses  available  for 
Canadian  storms  were  for  storms  that  crossed  the  border  from 
the  United  States  and  were  documented  in  "Storm  Rainfall  in 
the  U.S.".  Employing  the  storm  maximization  and 
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transposition  methods  of  Paulhus  and  Gilman  (1953),  Bruce 
obtained  preliminary  estimates  for  available  DAD  analyses 
for  southern  Ontario.  With  this  study  a  need  for  DAD 
analyses  became  apparent,  and  the  Department  of  Transport's 
Meteorological  Branch  (now  Atmospheric  Environment  Service) 
initiated  a  program  that  resulted  in  the  production  of  the 
"Storm  Rainfall  in  Canada"  series  (Atmospheric  Environment 
Service,  1 96 1  - ) .  This  is  a  continuing  series  with  a  number 
of  rainstorm  analyses  being  added  each  year  from  throughout 
Canada.  Although  the  series  mainly  emphasizes  rainstorms  in 
eastern  Canada,  slightly  over  60  such  analyses  are  for 
Alberta  storms. 

Canadian  research  advanced  in  the  1960s  when  a  number 
of  studies  were  conducted  to  examine  the  critical 
meteorological  conditions  for  maximum  floods  for  specific 
river  basins.  One  such  study  was  done  by  Bruce  and  Sporns 
(1963)  for  the  St.John  River  Basin,  New  Brunswick,  to 
provide  design  engineers  with  the  basic  information  needed 
to  estimate  maximum  flood  flows.  This  included  data  on  the 
physical  upper  limits  to  storm  rainfall,  winter  snow 
accumulation,  snowmelt  rates,  and  optimum  combinations  of 
snowmelt  and  rain.  DAD  analyses  of  63  storms  were  carried 
out  and  were  used  to  estimate  PMP  by  the  physical  approach. 
The  results  from  this  method  showed  very  good  agreement  with 
results  obtained  using  the  Hershfield  method  with  K  equal  to 
15  in  the  frequency  equation. 

A  second  study  was  done  by  Bruce,  Richards,  and  Sporns 
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(1965)  for  the  Portage  Mountain  Reservoir  on  the  Peace  River 
in  British  Columbia.  As  in  the  previous  study,  maximization 
of  the  major  recorded  rainstorms  was  carried  out  using  the 
physical  approach.  When  Hershfield' s  method  was  used  to 
corroborate  the  estimates  of  PMP,  it  was  found  that  the 
statistical  method  gave  only  slightly  lower  values  than 
those  obtained  by  the  physical  analysis. 

Another  study  was  performed  for  the  St.  Francois  and 
Chaudiere  river  basins,  Quebec,  by  Gagnon,  Pollock,  and 
Sparrow  (1970).  Similar  to  the  previous  two  studies,  this 
work  considered  only  the  meteorological  aspects  of 
rainstorms,  snow  accumulations,  and  snowmelt.  The  physical 
approach  was  also  used  in  examining  PMP  for  the  two  basins. 

A  total  of  77  large  rainstorms  (storms  with  point  rainfalls 
of  more  than  50  mm  (2  in.)  in  24  hours)  occurring  between 
1912  and  1964  were  analyzed.  Two  families  of  curves  or 
isochrones  were  obtained.  Figure  1.6  shows  the  maximum 
observed  rainfall,  and,  Figure  1.7  shows  the  maximum 
physically  possible  rainfall  as  functions  of  area.  The  five 
largest  storms  were  used  to  obtain  the  maximum  observed 
curves.  Figure  1.7  gives  the  maximum  possible  rainfall  for 
areas  of  130  to  12  950  km2  (50  to  5  000  sq .  mi.)  and  for 
durations  of  6  to  108  hours. 

In  Alberta,  most  of  the  work  on  estimating  PMP  has  been 
conducted  by  McKay  (1962,  1964,  1965b,  1966,  and  1968);  his 
efforts  have  concentrated  on  the  statistical  approach.  In 
one  of  the  first  of  these  studies  ( McKay , 1 965b ) ,  the  author 
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FIGURE  1.6.  Maximum  Observed  Rainfall  as  a  Function  of  Area 
for  the  Chaudiere  and  St.  Francois  Watersheds  in  Quebec 

(After  Gagnon  et  a  1 . ,  1970). 

obtained  estimates  of  PMP  for  the  prairie  provinces  by 
examining  24-hour  extreme  precipitation  values  for  191 
weather  stations  with  long-term  rainfall  records  from  1916 
to  1960.  The  characteristics  of  the  frequency  factor  K  were 
presented  along  with  maps  of  the  coefficient  of  variation 
and  the  mean.  The  author  also  provided  graphs  of  point  for 
conversion  of  24-hour  extreme  rainfall  to  other  rainfall 
durations  (i.e.,  1,  6,  and  12  hours).  The  author  concluded 
that  during  the  late  spring  and  throughout  the  summer  the 
foothills  of  southwestern  Alberta  have  the  greatest  storm 
potential  within  Alberta. 

In  another  study,  McKay  (1968)  examined  the 
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FIGURE  1.7.  Maximum  Possible  Rainfall  as  a  Function  of  Area 
for  the  Chaudiere  and  St.  Francois  Watersheds  in  Quebec 

(After  Gagnon  et  al.,  1970). 

meteorological  conditions  that  influenced  the  project  design 
and  PMF  on  the  Paddle  River,  Alberta.  Again  the  statistical 
approach  was  used  to  examine  the  15  available  severe  storms 
analyses,  including  those  storms  within  Alberta  that  the 
author  believed  could  be  transposable  to  the  Paddle  River 
area.  In  addition,  the  author  computed  snowmelt  rates  using 
generalized  snowmelt  equations  and  coefficients  and  combined 
these  estimates  with  the  snowmelt  associated  with  a  major 
spring  storm.  It  should  be  pointed  out  that  in  this  study  no 
critical  analysis  was  made  of  the  coefficients  in  the 
generalized  snowmelt  equations,  in  particular  of  their 
validity  for  the  basin  being  examined.  Estimates  of  the 
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water  yield  from  combined  rain  and  snowmelt  were  also 
computed,  as  shown  in  Figure  1.8.  The  solid  lines  in  the 


FIGURE  1.8.  Snowmelt  and  Precipitation  Computed  for  Various 
Combinations  of  Snowmelt  (After  McKay,  1968). 


figure  indicate  the  total  water  released  from  the  storm,  the 
storm  plus  1-day  melt,  the  storm  plus  2-day  melt,  and  so  on. 
The  broken  line  indicates  how  much  snow  can  be  expected  to 
be  on  the  ground  under  average  ablation  rate  conditions. 

Another  study  that  examined  critical  meteorological 
conditions  for  maximum  flow  was  that  presented  by  Buckler 
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and  Quine  (1971)  for  the  Pembina  Valley  north  of  Entwistle, 
Alberta.  The  authors  estimated  the  probable  maximum  snow 
accumulation  in  the  winter  season  and  calculated  the  melt 
caused  by  the  highest  possible  warm  spell  followed  by  a  3- 
day  rainstorm.  The  melt  rates  were  calculated  using  the 
generalized  snowmelt  equations  and  coefficients  developed  by 
the  U.S.  Army  Corps  of  Engineers  (1956).  They  estimated  the 
PMP  using  the  physical  approach  and  for  a  3~day  storm 
suggest  amounts  of  368  mm  (14.5  in.)  of  rainfall  in  the 
Pembina  Basin  (in  a  1  554  Km2  or  600  sq .  mi.  area).  The 
study's  DAD  curves  for  24,  48,  and  72  hours  are  shown  in 
Figure  1.9.  McKay's  (1968)  results  (based  on  Hershfield's 
approach)  are  also  displayed  on  the  figure. 

One  of  the  most  recent  publications  on  PMP  estimates  is 
by  Verschuren  and  Wojtiw  (1980),  in  which  both  the 
traditional  and  statistical  approaches  were  used  to  obtain 
estimates  of  PMP  for  the  six  main  river  basins  in  Alberta. 
The  results  of  this  work  are  discussed  further  under 
character i st i cs  of  PMP  in  Chapter  2. 

1.2.4  Historical  Flood  Studies  in  Alberta 

A  number  of  publications  have  documented  individual 
storms  that  have  caused  extensive  flooding  in  Alberta. 
Although  estimates  of  PMP  or  PMF  were  not  discussed  in  these 
studies,  it  is  important  to  recognize  the  work  done  in  this 
area,  for  some  of  these  storms  play  an  important  part  in  the 
PMP  estimates  and  also  demonstrate  that  rainstorms  do  cause 
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FIGURE  1.9.  DAD  Curves  for  24-, 48-, and  72-hour  Durations  for 
the  Pembina  Basin,  Alberta  (After  Buckler  and  Quine,  1971; 

McKay ,  1 968 )  . 

extensive  flooding.  One  such  study  was  done  by  Thompson 
(1976),  who  used  archive  records  to  reconstruct  three  severe 
storms  (June  1897,  July  1902,  and  June  1915)  on  the  Bow  and 
North  Saskatchewan  rivers.  The  author  concluded  that  the 
storm  type  that  produced  heavy  rainfall  in  all  three  cases 
was  a  cold  low.  The  rainfall  that  caused  record  floods  from 
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the  1897  and  1915  storms  was  of  extraordinary  intensity  6  , 
whereas  the  third  storm  (1902)  was  not  as  intense  7  .  The 
latter  event,  however,  was  preceded  by  an  exceptionally  wet 
period  during  which  the  soil  very  likely  became  saturated 
and  thereby  enhanced  storm  runoff.  The  resulting  peak 
discharge  from  the  1897  storm  is  believed  to  have  produced 
the  highest  flood  on  the  Bow  River  at  Calgary  dating  back  to 
at  least  1884.  The  storm  that  preceded  the  1915  flood  on  the 
North  Saskatchewan  River  is  believed  to  have  produced 
extremely  heavy  rainfall  amounts  over  a  wide  section  of  the 
Eastern  Slopes  of  the  Rockies  and  the  foothills  of  Alberta. 
The  instantaneous  discharge  recorded  on  the  North 
Saskatchewan  River  at  Edmonton  was  5  806  m3  per  sec  (205  000 
cubic  feet  per  second  (cfs))  on  28  June  1915,  and  the  daily 
discharge  of  4  644  m3  per  sec  (164  000  cfs)  recorded  on  the 
following  day  is  by  far  the  highest  on  record  for  that 
location  since  1911.  The  same  storm  produced  the  highest 
instantaneous  discharge  on  the  Red  Deer  River  at  Red  Deer 
since  records  began  in  1913,  the  fifth  highest  on  the  Bow 
River  at  Calgary  since  1884,  and  very  high  discharges  on  the 
Macleod  and  Athabasca  rivers. 

Another  example  of  a  study  where  rare  floods  were 
examined  is  the  report  by  Neill  (1965),  in  which  the  author 


6  The  4-day  rainfall  amount  in  the  1897  storm  exceeded  a 
1:100  year  event,  while  in  the  1915  storm  the  2-day  rainfall 
amount  exceeded  a  1:50  year  value  at  Nordegg  and  Lovett. 

7  The  heaviest  1-  and  3-day  amounts  recorded  during  the 
storm  are  less  than  a  1:10  year  event. 
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tabulated  rare  flood  discharge  data  from  47  selected 
hydrometric  stations  in  Alberta.  In  addition,  the  author 
plotted  the  100-year  mean  daily  flows  against  drainage  areas 
and  obtained  a  general  empirical  relationship  of  the  form 


Q 


100 


F  A 


0.8 


(1.5) 


where  Qiqo^s  the  100'vear  mean  daily  flow  (cfs),  F  is  a 
coefficient,  and  A  is  the  drainage  area  ( sq .  mi.).  The 
author  also  calculated  the  values  of  F  for  each  station  and 
depicted  this  on  a  map  of  Alberta.  A  plot  of  maximum 
recorded  flows  against  drainage  areas  indicated  an  upper 
envelope  curve  with  the  equation 
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A  more  recent  example  of  a  rainstorm  that  produced 
severe  flooding  occurred  on  7-8  June  1964  (Warner,  1973)  in 
the  headwaters  of  many  streams  in  Montana  and  southwestern 
Alberta.  The  main  flood  area  in  Alberta  occurred  in  Waterton 
Lakes  National  Park.  Flooding  was  primarily  attributed  to 
extreme  rainfall  (unofficially,  over  250  mm  (10  in.)  in 
Waterton  Park)  influenced  by  the  right  combination  of 
antecedent  conditions:  below-normal  temperatures  that 
delayed  the  usual  snowmel t -runoff  pattern  (from  March  to 
May),  and  above-normal  precipitation  in  May  that  resulted  in 
large-scale  melting  of  the  snowpack  in  the  latter  part  of 
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May,  which  continued  into  June  at  a  sustained  high  rate. 
Although  this  storm  occurred  mainly  in  United  States,  with 
Alberta  receiving  only  a  small  segment,  the  possible 

occurrence  of  such  a  storm  further  north  cannot  be  ruled 
out . 

Flood-producing  rainstorms  are  not  restricted  to  the 
southern  part  of  the  province.  For  example,  Warner  and 
Thompson  (1974)  presented  a  report  on  the  11-12  June  1972 
rainstorm  in  which  more  than  150  mm  (6  in.)  of  rain  fell 
over  parts  of  the  Peace  River  Basin  southwest  of  Grande 
Prairie,  resulting  in  record  flows  in  nearly  all  streams  in 
that  area.  Other  documented  examples  of  rainstorms  causing 
flooding  conditions  in  northern  Alberta  include  studies  by 
McKay  (1966)  of  the  30  July  to  1  August  1953  storm  on  the 
Paddle  River;  by  Froelich  (1967)  of  the  June  1965  storm  in 
the  Peace,  Athabasca,  and  North  Saskatchewan  basins;  and  by 
Mustapha  ( 1970b!  of  the  27  June  to  1  July  1970  storm  in  the 
Lesser  Slave  Lake  and  Lac  La  Biche  forestry  regions. 

A  more  recent  publication  on  floods  is  the  study  by 
Mustapha  et  al.  (1981)  in  which  the  authors  presented  a 
compilation  of  historical  flood  data  and  information 
concerning  recorded  floods  within  the  North  Saskatchewan 
River  Basin.  Data  from  18  selected  hydrometric  stations  in 
the  basin  with  varying  streamflow  records  between  1911  to 
1978  were  selected  as  representat i ve  of  floods  in  the 
mountains,  foothills,  plains,  and  mainstream  of  the  North 
Saskatchewan  River .  In  the  report  the  authors  examined  the 
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maximum  annual  flood  discharge,  the  causes  of  floods,  the 
effect  of  ice  on  recorded  stages,  the  historical  flood 
levels,  flood  damages,  and  flood  frequency  analyses.  In  the 
mountain  region  most  of  the  volume  of  annual  runoff  and  the 
annual  flood  peaks  are  due  to  snowmelt  in  the  spring  and 
summer.  In  the  foothills,  although  snowmelt  contributes  to 
runoff,  the  major  flood  peaks  are  generated  by  heavy 
rainfall  in  this  region.  Runoff  from  mountain  snowmelt 
combined  with  runoff  from  major  storms  in  the  foothills 
generally  produced  the  largest  flood  peaks  on  the  North 
Saskatchewan  River.  In  this  region,  cold  lows  produce  the 
major  floods.  In  the  plains  region,  snowmelt  is  the  major 
contributing  factor  to  the  maximum  annual  flood  peak,  and 
this  usually  occurs  in  April  or  early  May.  On  occasion,  the 
early  spring  runoff  from  this  region  is  due  to  a  combination 
of  rain  and  snowmelt.  In  some  years  the  maximum  annual  flood 
peaks  result  primarily  from  heavy  summer  rainfall. 

Mustapha  et  al.  (1981)  also  discuss  urban  storms,  in 
particular  the  storm  that  occurred  in  the  Edmonton  area  on 
10-11  July  1978  causing  over  $1.3  million  damage.  Although 
this  storm  did  not  produce  very  large  recorded  maximum  point 
precipitation,  it  did  produce  extensive  flooding  in  the  city 
of  Edmonton  and  surrounding  towns.  Urban  flooding  is 
generally  caused  by  highly  localized  rainstorms  whose 
intensity  and  subsequent  runoff  exceed  the  design  capacity 
of  the  urban  storm  sewer  system. 
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1.2.5  Snowmelt  Studies 

As  was  discussed  earlier,  some  major  floods  result  from 
melting  snowpacks  or  from  snowmelt  combined  with  rain.  To 
estimate  the  maximum  floods  in  these  regions  it  is  necessary 
to  consider  the  contributions  of  snowmelt  water  to  major 
floods.  This  requires  (1)  determination  of  the  maximum 
seasonal  snow  accumulation  and  (2)  estimation  of  the 
critical  melting  rate  of  the  snowpacK. 

Three  methods  have  been  used  to  estimate  the  upper 
limits  of  snow  accumulation  on  watersheds:  (a)  partial 
season  (or  synthetic  season)  method,  (b)  snowstorm 
maximization  method,  and  (c)  statistical  method. 

In  the  partial  season  method  an  estimate  of  the 
physical  upper  limit  of  seasonal  snow  accumulation  is  made 
by  assuming  that  the  largest  snowfalls  for  discrete 
intervals  (i.e.,  2,  4  days,  or  1  month)  during  the  recorded 
period  can  be  combined,  regardless  of  their  year  of 
occurrence,  to  produce  a  synthetic  year  of  very  high 
snowfall.  The  validity  of  such  a  synthetic  year  has  been 
questioned,  since  different  years  are  being  combined.  If  a 
long  recorded  period  is  available,  the  upper  limit  obtained 
by  this  method  would  rarely  be  exceeded;  however,  for  data 
bases  with  a  short  duration  record  this  is  not  necessarily 
true,  and  hence  caution  must  be  applied  in  the 
i nterpretat i on  of  such  results. 

In  the  second  method,  the  maximum  snowfall  is  obtained 
by  snowstorm  maximization  of  all  the  snowfalls  over  the 
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basin  for  the  entire  winter.  A  factor  is  used  to  estimate 
the  maximum  moisture  flow  causing  the  snowfall.  This  factor, 
similar  to  that  used  for  rainfall,  is  equal  to  the  ratio 
between  the  maximum  possible  moisture  content  for  the  time 
of  year  when  the  snow  fell  and  the  moisture  observed  in  the 
airmass  that  actually  produced  the  snowfall. 

The  last  approach,  the  statistical  method,  estimates 
the  snowpack  on  the  basis  of  a  long  return  period,  such  as 
100  or  1000  years,  by  applying  the  extreme  value  theory  of 
Gumbel  (1954)  to  the  maximum  annual  water  equivalent  of  snow 
cover.  This  approach  is  applied  to  point  values,  which  are 
then  translated  into  an  areal  estimate  based  on  the 
assumption  that  there  is  little  variation  in  the  water 
equivalent  of  the  snow  cover.  This  assumption  is  usually 
correct  for  the  plain  areas;  however,  for  mountainous 
regions  it  may  not  be  applicable.  The  approach  is  also 
highly  dependent  on  the  recorded  period,  and  caution  should 
be  used  when  obtaining  estimates  from  stations  with  a  short 
record  length. 

Some  studies  of  eastern  Canadian  basins  (Gagnon  et  al . , 
1970)  suggest  that  the  result  obtained  by  the  snowstorm 
maximization  method  corresponds  fairly  closely  to  the  value 
derived  by  the  partial  season  method;  however,  none  of  these 
methods  is  entirely  satisfactory.  Perhaps  the  second 
approach,  the  winter  snowstorm  maximization,  has  the  most 
credibility,  but  there  is  the  problem  of  compounding 
unlikely  events  by  assuming  that  all  winter  storms  in  a 
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season  occur  with  maximum  water  vapor  content  in  the  snow- 
producing  air  mass.  Research  on  new  and  better  methods  of 
estimating  maximum  snow  accumulation  is  needed. 

For  estimating  critical  snowmelt  rates,  two  main 
approaches  are  commonly  used.  The  first  is  the  degree-day 
method,  while  the  second  is  to  apply  generalized  snowmelt 
equations  based  on  energy  balance  considerations  (U.S.  Army 
Corps  of  Engineers,  1956). 

In  the  degree-day  method  the  melt  (M)  is  expressed  as  a 
linear  relationship  of  the  daily  air  temperature  and  can  be 
wr i t  ten  as 


M  ■  CIX  (1.7) 

where  C  i  s  an  empirically  determined  coefficient  and  X!  Ta  is 
the  sum  of  positive  daily  air  temperatures  ( °C )  for  a 
designated  period.  Either  maximum  or  mean  temperatures  can 
be  used  for  .  There  is  some  physical  justification  for 
this  approach,  since  air  temperature  is  reasonably  well 
correlated  at  a  particular  time  and  place  with  atmospheric 
factors  that  affect  melt  rates,  such  as  solar  radiation  and 
vapor  pressure.  This  is  a  simple  and  useful  approach  for 
calculating  snowmelt  but  has  some  limitations  because  there 
is  considerable  variability  in  the  coefficient  C  from  basin 
to  basin,  year  to  year,  and  time  to  time  within  the  snowmelt 
season,  depending  upon  factors  that  cannot  be  represented  by 
air  temperature  alone  (i.e.,  albedo  of  the  snow  surface, 
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wind,  dew  point,  and  radiation).  To  date,  little  effort  has 
been  expended  to  justify  the  usage  of  the  mean  or  maximum 
temperature  in  the  interpretation  of  these  quantities. 

The  other  and  more  rigourous  approach  is  to  apply 
generalized  snowmelt  or  energy  balance  equations,  which  can 
be  writ  ten  as 
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where  M  is  the  total  snowmelt,  Mrg  is  the  short-wave 
radiation  melt.M^  is  the  long -wave  radiation  melt,  ^ce  is 
the  melt  due  to  convective  heat  transfer  from  the  atmosphere 
and  to  latent  heat  of  water  vapor  condensing  onto  the  snow 
surface,  M  is  the  melt  due  to  heat  of  raindrops,  and  M  is 
the  melt  by  heat  conduction  from  the  ground.  Approximation 
of  each  of  the  various  melts  have  been  described  by  the  U.S. 
Army  Corps  of  Engineers  (1956).  Because  these  approximation 
were  developed  for  a  specific  basin,  the  user  must  calibrate 
the  coefficients  for  the  basin  being  studied. 

For  both  the  energy  balance  equations  and  the  degree- 
day  method,  curves  of  highest  maximum  (or  mean)  temperatures 
for  the  snowmelt  season  need  to  be  developed  for  a 
representat i ve  long-term  record  station  or  stations  within 
the  basin.  In  using  the  degree-day  method,  the  temperature 
sequence  is  all  that  is  required.  For  the  energy  balance 
procedure,  critical  values  of  other  meteorological  factors 
(insolation  and  albedo,  dew  point  temperatures,  and  wind) 
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are  needed.  The  energy  balance  equations  have  been  employed 
in  Alberta  by  various  authors  (McKay,  1965b,  1968;  Buckler 
and  Quine,  1971;  Storr ,  1978)  in  computing  snowmelt.  In  many 
of  these  studies,  the  coefficients  given  by  the  U.S.  Army 
Corps  of  Engineers  were  used,  and  coefficients  were  not 
calibrated  for  the  Alberta  river  basins. 

1.2.6  Flood  Estimates  and  PMF 

To  determine  character i st i cs  of  a  flood  that  would 
result  from  a  given  rainstorm  and/or  snowmelt  period,  it  is 
necessary  to  estimate  the  percentage  of  basin  water  that 
will  appear  as  surface  runoff.  Some  of  the  rain  or  snowmelt 
water  will  infiltrate  the  soil  and  will  not  contribute 
directly  to  -the  flood  rise;  instead,  it  may  recharge 
groundwater  or  be  stored  in  the  soil.  The  water  that  flows 
quickly  into  the  stream  channels  mainly  overland  or  as 
interflow  is  known  as  direct  runoff.  It  is  this  volume  that 
must  be  estimated.  Direct  runoff  varies  from  time  to  time 
within  a  basin  and  from  basin  to  basin.  Each  basin  has  a 
characteristic  response  depending  on  factors  such  as  the 
permeability  of  the  soils,  the  vegetation,  the  slopes  of 
main  land  areas  of  the  basin,  the  amount  of  the  basin  in 
swamp  area  or  lakes,  and  the  amount  of  small  depression 
storage  in  the  basin.  Within  a  given  basin  the  volume  of 
direct  runoff  from  a  given  amount  of  rain  varies  with  the 
season,  the  antecedent  conditions,  and  the  duration  and 
intensity  of  storm  rainfall.  Since  all  of  these  factors  are 
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complex  and  some  are  i nter re  1 ated ,  it  is  very  difficult  to 
estimate  runoff  volumes.  For  a  particular  river  basin  with 
records  of  streamflow  and  precipitation,  a  common  procedure 
has  been  to  develop  multiple  variable  rai nf a  1 1 -runoff 
correlations  either  graphically  or  by  empirical  formulas 
(World  Meteorological  Organization,  1975). 

Graphical  relationships  have  been  published  by  the  U.S. 
Geological  Survey  (Chow,  1964)  in  f lood- frequency  reports 
for  various  states.  From  these  the  magnitude  and  probable 
recurrence  interval  of  floods  may  be  determined  at  any  place 
on  a  stream,  within  the  limits  of  the  basic  data.  There  are 
two  di scharge- frequency  curves  that  are  commonly  used:  the 
first  expresses  the  flood  di scharge- t i me  relationship  and 
shows  variation  of  peak  discharge,  which  is  expressed  as  a 
ratio  of  the  mean  annual  flood,  as  well  as  recurrence 
interval  8 9  ;  the  second  curve  relates  the  mean  annual  flood 
9  to  the  area  alone  or  to  the  size  of  area  and  other 
significant  basin  characteristics. 

From  about  1913  to  the  1930s,  what  is  now  Known  as  the 
first  flood  frequency  formula  was  developed  by  Fuller  (1914) 
from  the  discharge  records  of  hundreds  of  streams.  The 


8  The  average  interval  of  time  within  which  the  given  flood 
will  be  equalled  or  exceeded.  A  flood  having  a  recurrence 
interval  of  10  years  is  one  that  has  a  10%  chance  of 
recurring  in  any  year.  Likewise,  a  50-year  flood  has  a  2% 
chance,  and  a  100-year  flood  has  a  1%  chance,  of  recurring 
in  any  year. 

9  The  Gumbel  distribution  is  used  to  express  the  mean  annual 
flood,  which  is  defined  as  the  flood  having  a  recurrence 
interval  of  2.33  years. 
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formula  is: 

Qt  =  Q  (1  +  0.8  log10  T)  (1>9) 

where  QT  is  the  flood  of  return  period  T  (in  years),  and  Q 
is  the  mean  highest  annual  flood.  Fuller  suggested  using  the 
above  formula  to  estimate  floods  with  as  long  a  return 
period  as  1000  years.  As  late  as  1932,  T  =  1000  was 
recommended  as  the  spillway  design  flood  for  certain  classes 
of  earth  dams,  together  with  a  generous  freeboard  safety 
factor.  This  suggestion  (i.e.,  T  =  1000)  is  based  not  only 
on  extrapolation  of  what  appeared  to  be  a  law,  but  also  on 
transposition  and  envelopment.  It  is  now  Known  that 
extrapolation  to  such  return  periods  is  uncertain. 

With  more  discharge  records,  directly  enveloping  the 
records  of  peak  discharges  (normalized  for  area)  rather  than 
assigning  specific  frequencies  became  more  popular  in 
deriving  the  maximum  flood.  The  most  famous  is  the  Myers' 
equation  (Chow,  1964)  in  which  peak  discharge  is 
proportional  to  the  square  root  of  the  area  drained  and  to  a 
coefficient  that  varies  with  region  and  geology.  Myers' 
formula  is: 

0.5 

Qm  =  c  A-  (1.10) 

where  Qm  is  the  maximum  discharge  (cfs);  C  is  an  empirical 
coefficient  varying  with  climate  and  geological 
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characteristics;  and  A  is  the  drainage  area  ( sq .  mi.). 

The  second  emphasis  in  the  period  from  about  1914  to 
1935  was  on  the  use  of  probability  theory  to  estimate 
magnitude  of  the  flood.  Curves  were  derived  from  past 
records  on  a  stream,  and  the  frequency  was  determined  for 
the  occurrence  of  a  given  flood.  From  data  with  a  few  years 
of  records  (i.e.,  seldom  exceeding  20  years  and  rarely 
exceeding  30  or  40  years),  probability  curves  were 
extrapolated  to  estimate  the  average  expected  discharge  for 
longer  time  durations  (i.e.,  once  in  1  000,  5  000,  or  10  000 
years).  This  theory  proved  to  be  inadequate  in  that  some  of 
the  floods  that  occurred  in  the  1930s,  were  estimated  by  the 
probability  method  to  have  a  return  period  of  once  in 
millions  or  even  billions  of  years,  and  not  the  usual 
assumed  1  000  to  10  000  years.  One  such  example  is  the  flood 
in  1935  on  the  Republican  River  in  Nebraska,  which  was  over 
10  times  larger  than  had  occurred  on  that  river  during  the 
previously  recorded  40  years.  The  probability  method  would 
not  have  made  provision  for  an  adequate  spillway  for  that 
f 1 ood . 

The  nature  of  flood  analysis  changed  in  the  1930s,  with 
the  increased  knowledge  of  floods  gained  by  superimposing 
the  enveloping  curves  of  record-breaki ng  floods  that 
occurred  between  1890  to  1939  (Figure  1.10)  (Craeger,  1939). 
With  more  gauging  stations  and  longer  records,  the  curves 
depict  an  upward  trend  in  flood  flows  (there  is  no  evidence 
that  these  changes  were  climatological).  The  increased  data 
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Drainage  Area  (sq.  mi.) 


FIGURE  1.10.  Enveloping  Curves  of  Known  Floods  in  the  United 
States  up  to  a  Given  Time  Period  (After  Craeger,  1939). 

significantly  raised  the  records  of  flood  flows  for  all 
drainage  areas.  For  instance,  the  maximum  flood  for  an  area 
of  259  Km2  (100  sq .  mi.)  that  in  1890  had  a  limit  of  2.2  m3 
per  sec  per  Km2  (200  cfs  per  sq .  mi.)  was  increased  about 
ninefold  in  1939  to  19.7  m3  per  sec  per  Km2  (1800  cfs  per 
sq.  mi.).  For  larger  areas  the  increase  between  the  1890  and 
1939  results  was  less,  for  example,  the  flood  for  25  900  Km2 
(10  000  sq .  mi.)  being  about  double. 

Many  rainfal 1 -runoff  formulas  have  been  derived  for 
computing  flood  discharges  since  the  1940s,  but  most  are 
considered  inadequate  for  engineering  design.  Difficulties 
in  the  application  of  empirical  relationships  arise  not  so 
much  from  the  empiricism  of  equations  but  more  from  the  lack 
of  knowledge  of  the  exact  conditions  under  which  they  may  be 
applied.  Peak  flow  formulas  can  be  classified  according  to 
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the  number  of  variables  used  in  the  equation.  Listed  below 
are  examples  of  the  functional  relationships  of  the  numerous 
peak  formulas  compiled  by  Gray  (1970)  in  which  flood  flows 
are  related  to 
(  1  )  drainage  area 

(2)  drainage  area  and  frequency 

(3)  area,  rainfall,  and  time  paramaters 

(4)  basin  width  and  rainfall  parameters 

(5)  basin  width,  rainfall  parameters,  and  slope 

(6)  area,  slopes,  and  rainfall  parameters 

(7)  area,  elevation,  length,  and  lakes. 


A  comparison  of  some  of  the  flood  formulas  is  presented  in 
Figure  1.11  (Gray,  1970).  No  correlation  relation  between 
the  various  formulas  is  given  and  hence  it  is  difficult  to 
assess  their  reliability. 

Another  important  character i st i c  of  a  flood  to  consider 
in  spillway  design  is  the  hydrograph  shape.  The  numerous 
papers  published  on  the  time  distribution  of  runoff  not  only 
indicate  the  importance  of  hydrograph  shape,  but  suggest 
that  the  techniques  are  far  from  being  completely 
definitive.  Where  sufficient  data  on  streamflow  are  not 
available  to  enable  the  shape  of  the  flood  hydrograph  to  be 
determined,  simplified  geometric  forms  have  been  used  for 
approximation.  It  is  obvious  that  simple  geometrical  forms 
are  only  first  approximations  and  that  more  attention  should 
be  directed  to  the  character i st ics  of  the  drainage  area. 
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CATCHMENT  AREA  (  SQ,  KM.) 


FIGURE  1.11.  Comparison  of  Some  of  the  Flood  Formulas  in 

General  Use  (After  Gray,  1970). 

The  unit  hydrograph  method  originally  presented  by 
Sherman  (1932)  has  been  the  most  widely  used  method  of  flood 
wave  form  computation.  This  method  has  been  modified  by  many 
authors,  but  most  of  the  basic  principles  have  not  changed. 
The  unit  hydrograph  is  defined  as  a  hydrograph  derived  from 
storm  rainfall  of  a  specified  duration,  where  the  volume  of 
surface  runoff  accounted  for  is  of  unit  depth  on  the  basin. 
The  unit  hydrograph  method  and  its  modifications  are  based 
on  assumptions  not  entirely  fulfilled  in  nature.  Johnstone 
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and  Cross  (1949)  give  a  fitting  characterization  of  this 
situation:  All  these  propositions  (unit  hydrograph 

assumptions)  are  empirical.  It  is  not  possible  to  prove  them 
mathematically.  In  fact,  it  is  a  rather  simple  matter  to 
demonstrate  by  hydraulic  analysis  that  not  a  single  one  of 
them  is  mathematically  accurate.  Fortunately,  Nature  is  not 
aware  of  this." 

The  difficulties  encountered  with  collection  of  data 
required  for  derivation  of  unit  hydrographs  caused  many 
workers  to  seek  further  modifications  of  the  method.  Lack  of 
sufficient  precipitation  and  runoff  data  inspired  Snyder 
(1938)  to  develop  a  synthetic  hydrograph  method.  In  this 
method  the  individual  elements  of  the  hydrograph  are 
determined  by  computation.  Peak  discharge  for  a  given 
duration  was  derived  by  Snyder  from  the  equation 

%  =  640  CP  Ad  /  S  (i.ii) 

where  Qm  i s  the  peak  discharge  (cfs);  the  coefficient 

accounting  for  flood  wave  and  storage  conditions;  Ad  is  the 

watershed  area  ( sq .  mi.);  and  t  the  lag  time  (hr).  Values 

P 

for  range  from  0.4  to  0.8  and  generally  indicate 

retention  or  storage  capacity  of  the  watershed.  Application 

to  streams  in  other  parts  of  the*  world  requires  different 

values  of  the  Cp  coefficient.  The  relation  between  the  time 

of  the  center  of  precipitation  period  and  the  peak  of  the 

hydrograph  rise  ( t  )  is  given  by  the  formula 

P 
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,  ,  0.3 

tp  =  ct  (L  Lc)  (1.12) 

where  ct  is  a  coefficient  from  1.8  to  2.2  (dependent  on 
slope  of  the  drainage  area);  L  is  the  total  length  of  the 
main  stream  (miles);  and  Lc  is  the  length  of  main  stream 
from  gauging  station  to  cross  section  nearest  to  the 
centroid  of  the  drainage  area  (miles). 

The  U.S.  Department  of  Agriculture  Soil  Conservation 
Service  (SCS,  1969)  has  developed  a  dimensionless  unit 
hydrograph  where  ordinate  and  abscissa  values  are  expressed 
as  a  nond i mens i on  1  ess  relationship.  Peak  runoff  and  rise 
time  depend  on  basin  area,  duration  of  rainfall  excess, 
length  of  the  main  channel  and  slope  of  the  watercourse.  The 
peak  flow  of  the  unit  hydrograph  can  be  calculated  using  the 
f o 1  lowing  f ormu 1  a : 


Qp  =  4-84  A  /  PR 


(1.13) 


where  Qp  is  the  peak  runoff,  A  is  the  area  of  the  watershed 
( sq .  mi.),  and  pR  is  rise  time  (hours). 

The  Tennessee  Valley  Authority  (1973)  has  developed  a 
unit  hydrograph  using  a  double  triangle  defined  at  four 
points.  Major  elements  correspondi ng  to  these  four  values 
are  basin  area,  shape,  length  of  the  main  watercourse, 
drainage  density,  percentage  wooded  area  and  soil  type. 

Hoang  (1977)  has  modified  and  applied  this  method  to  small 
watersheds  in  southern  Quebec.  The  peak  flow  is  calculated 
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using  the  following  formula: 


where  q  is  the  peak  flow  of  a  unit  hydrograph  resulting 
from  a  net  rainfall  of  2  hours,  k  is  a  regional  constant, 
A  is  the  area  of  the  watershed,  and  x  is  a  regional 
constant.  Hoang  (1977)  has  developed  by  regression  the 
values  of  the  regional  constants  for  southern  Quebec. 

It  is  interesting  to  compare  the  different  formulas 
used  to  calculate  peak  flow.  In  all  of  these  formulas,  area 
is  the  common  element.  In  the  Snyder  and  SCS  formulas,  the 
flow  varies  linearly  according  to  area,  while  the  Hoang 
model  contains  an  exponent  of  area  that  approaches  unity. 

Empirical  formulas  have  been  applied  to  watersheds  in 
western  Canada.  McKay  and  Stichling  (1961)  obtained  a 
relationship  for  the  envelope  curve  of  extreme  flooding  on 
the  prairies,  for  which  peak  discharge  and  drainage  area 
(Figure  1.12)  are  determined  by  Equation  1.10,  with  C  equal 
to  400. 

Another  example  is  the  work  of  Godwin  (1975),  who 
carried  out  regionalized  estimates  of  PMF  for  72  sites  in 
the  prairie  provinces,  primarily  for  the  Saskatchewan  Nelson 
Basin.  The  estimates  were  made  for  areas  with  floods  from 
mountain  and  prairie  streams.  The  results  showed  that  for 
both  types  of  streams  the  maximum  mean  daily  peak  flow  per 
unit  area  (Figure  1.13)  decreased  exponentially  with  an 
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FIGURE  1.12.  Envelope  Curves  of  Extreme  Floods  on  Prairies 
(After  McKay  and  Stichling,  1961). 

increase  in  the  drainage  area,  with  the  mountain-fed  streams 
having  a  higher  magnitude.  The  relationships  for  these 
curves  also  take  the  form  of  Equation  1.10. 

1.2.7  Computer  Simulation  Modelling  of  PMF 

The  concept  of  a  mathematical  model  for  simulation  of 
synthesized  streamflow  has  become  prominent  in  the  past  20 
years,  primarily  with  the  advancement  of  electronic  computer 
technology.  The  basic  concept  involves  simulating  a 
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FIGURE  1.13.  Mean  Daily  Peak  Flow  Versus  Drainage  Area  for 
Probable  Maximum  Floods  (After  Godwin,  1975). 


mathematical  model  comprising  the  physical  processes 
involved  in  the  hydrologic  cycle.  Such  a  model  may  include 
rainfall  or  snowmelt,  effects  of  tr anspi r at  ion ,  interception 
losses,  infiltration  rates,  soil  moisture,  depression 
storage,  and  groundwater  storage.  Usually  the  flow  is  then 
separated  into  its  components  of  surface,  subsurface,  and 
base  flow.  Various  versions  of  computer  models  have  been 
developed  with  additions  or  updates  of  the  various  runoff 


processes . 


■ 


46 


At  present,  there  are  more  than  20  major  hydrologic 
simulation  models  (U.S.  Army  Waterways  Experiment  Station, 
1974)  available  for  representat ion  of  the  hydrological 
processes.  These  can  be  subdivided  into  three  main  classes: 
(a)  rai nf a  1 1 -runoff  event  simulation  models,  (b)  continuous 
streamflow  simulation  models,  and  (c)  urban  runoff 
simulation  models.  Some  of  the  more  common  models  (with 
their  abbreviated  notation)  are  the  following: 

a)  Ra i nf a  1 1  -  runoff  event  simulation  models 


HEC-1  Flood  Hydrograph  Package- ( HEC- 1 ) 

Computer  Program  for  Project  Hydrology- ( TR-20 ) 
USGS  Rainfal 1 -runoff  Mode  1  - ( USGS ) 

Hydrologic  Model  Computer  Language- ( HYMO ) 

Storm  Water  Management  Model  - ( SWMM ) 


b)  Conti  nous  streamflow  simulation  models 


Antecedent  Precipitation  Index  Model -(API) 

1970,  1973,  1974  Revised  Watershed  Hydrology- ( USDAHL ) 
Stanford  Watershed  Model  IV-(SWM-IV) 

Kentucky  Watershed  Model -(KWM) 

Se 1 f -opt imi zi ng  Hydrologic  Simulation  Mode  1  - ( OPSET ) 
Hydrocomp  Simulation  Program- ( HSP ) 

Texas  Watershed  Model  - ( TWM ) 

National  Weather  Service  Runoff  Forecast  System- ( NWSRFS ) 
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Streamflow  Synthesis  and  Reservoir  Regu 1  at i on- ( SSARR ) 
Ohio  State  University  Version  of  SWM- 1 V- ( OSUSWM ) 

c)  Urban  runoff  simulation  models 

University  of  Cincinnati  Urban  Runoff  Model- (UCUR) 
Chicago  Hydrograph  Method- ( NERO ) 

Quantity  and  Quality  of  Urban  Runof f- ( STORM ) 

Road  Research  Laboratory  Model - ( RRL ) 

MIT  Catchment  Model  - (MITCAT ) 

A  review  of  these  models  (U.S.  Army  Waterways 
Experiment  Stations,  1974;  World  Meteorological 
Organization,  1975;  Viessman  et  al.,  1977)  indicated  a 
tremendous  diversity  in  scope  and  purpose,  mathematical 
detail,  system  elements,  and  hydrologic  phenomena  being 
modelled,  size  of  the  system  that  can  be  handled,  data 
requi rements ,  and  computer  output.  These  diversities  are  the 
result  of  the  varying  conditions  and  objectives  that  govern 
the  design  and  evaluation  of  individual  drainage  systems. 

The  state  of  development  of  these  models  also  varies 
significantly.  Some  models  have  been  developed  and  verified 
extensively,  others  have  been  developed  but  not  verified, 
and  some  have  been  developed  conceptually  but  not  carried  to 
the  point  of  application.  In  addition,  because  no  standards 
exist  for  evaluating  and  comparing  models,  differing 
criteria  have  been  used  in  the  evaluations.  The  rigor  with 
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which  the  models  have  been  tested  varied  greatly,  from 
intuitive  judgments  to  graphical  comparisons  and  more 
demanding  statistical  analyses. 

To  provide  a  quick  overview  of  the  features  of  these 
models,  a  summary  indicating  the  mathematical  formulations 
used  in  the  models  is  presented  in  Figure  1.14.  It  must  be 
remembered  that  most  of  the  models  are  being  revised, 
expanded,  and  improved  continuously  by  the  model  developers. 
A  major  drawback  with  many  of  the  existing  models  is  that 
they  have  not  been  developed  for  large  rural  drainage  basins 
or  have  not  incorporated  snowmelt  computations.  The  latter 
is  important  for  estimating  the  PMP  and  PMF  for  watersheds 
in  Alberta.  Other  models  have  extensive  data  requi rements , 
are  difficult  to  implement  or  modify,  and  may  be  very 
expensive  for  simulation. 

The  following  are  some  of  the  above  models  and  other 
minor  models  that  have  been  employed  to  estimate  floods 
design  and  in  some  cases  PMF  estimates: 

(1)  Stanford  Watershed  Model-used  in  Quebec  for  certain 
small  watersheds, 

(2)  SSARR  Model-used  by  Alberta  Environment  for  design 
purposes  on  the  Red  Deer  and  Paddle  rivers, 

(3)  CEQUEAU  Model-used  to  determine  the  design  floods  of 
several  rivers  in  northern  Quebec  (Morin  et  al.,  1979), 

(4)  Hydro-Quebec  Model -used  to  make  real-time  forecasts  of 
floods  in  Quebec  (Bisson,  1979), 

(5)  UBC  Model  (University  of  British  Col umbi a ) -used  to 
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Infiltration  and  Losses 

1.  Hoi  tan's  Equation 

2.  Horton's  Equation 

3.  Phillip's  Equation 

k.  SCS  Curve  Ho.  Method 

5.  Variable  Loss  Rate 

6.  Standard  Capacity  Curves 

Unit  Hydrograph 

l.  Input 

2.  Clark's 

3.  Snyder's 

4.  2-Parameter  Gamma  Response 

5.  Dimensionless  Unit  Hydrograph 

River  Routing 

1.  Hydraulic 

2.  Muskingum 

3.  Tatum 

4.  Straddle-Stagger 

5.  Modified  Puls 

6.  Working  R  and  D 

7.  Variable  Storage  Coefficients 
B.  Convex  Method 

9.  Translation  Only 
10.  Lumped  Parameter  Approach 

Base  Flow 

1.  Input 

2.  Constant  Value 

3-  Recession  Equation 
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1.14.  Comparison  of  Mathematical 
Models  (After  Viessman  et  al., 


Formulations  in 
1977)  . 


calculate  the  PMF  in  10  subwatersheds  of  the  Peace  River 
( Fawkes ,  1979  )  , 

(6)  SWMM  Model _used  to  establish  the  design  flood  in  urban 
areas  for  small  basins  (Leclerc,  1979). 

Many  of  the  models  described  above  are  restricted  to 
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small  river  basin  studies  since  they  were  developed  and 
modified  for  these  needs.  The  two  models  that  appear  to  have 
been  used  in  large  river  basins  studies  are  the  SSARR  and 
UBC  models.  Mustapha  and  Ojamaa  (1975)  used  the  UBC  SASK6 
simulation  model  to  obtain  estimates  of  PMF  for  the  Red  Deer 
River  at  the  Raven.  Using  DAD  analyses  of  four  rainstorms, 
the  authors  estimated  the  PMP  to  be  292  mm  (11.5  in.)  near 
the  Raven  for  a  48-hour  duration,  and  distributed  this 
precipitation  on  2  consecutive  days,  with  229  mm  (9  in.)  in 
the  first  24  hours.  The  PMP  was  converted  to  flood  flows  by 
use  of  the  UBC  model;  the  resulting  PMF  hydrograph  for  Raven 
is  shown  in  Figure  1.15  (Muzik,  1975)  and  indicates  an 
instantaneous  peak  discharge  of  4  814  m3  per  sec  (170  000 
cf  s )  . 

The  SSARR  model  has  been  and  still  is  being  used  for 
various  studies  by  Alberta  Environment;  however,  only 
internal  reports  have  been  prepared  on  its  usage.  In  the  mid 
1970s,  the  Alberta  Flow  Forecasting  Branch  of  Alberta 
Environment  conducted  an  investigation  of  several  high- 
powered  hydrologic  computer  models  in  order  to  choose  a 
model  for  real-time  flow  forecasting  in  Alberta.  The 
watershed  models  included  in  the  investigation  were  HEC-1, 
Modified  Stanford,  UBC  Watershed  Model,  and  SSARR.  The  study 
concluded  that  while  any  of  the  models,  if  properly 
calibrated,  could  produce  consistently  better  results  than  a 
unit  hydrograph  approach,  the  SSARR  Model  was  the  most 
flexible,  comprehensive,  and  easiest  to  operate  of  the  four. 
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FIGURE  1.15.  Probable  Maximum  Flood  Hydrograph  at  Raven  on 
the  Red  Deer  River  (After  Muzik,  1975). 


Because  of  its  use  by  Alberta  Environment  and  vast 
flexibility,  the  SSARR  Model  is  used  in  this  thesis. 


1.3.0  Sources  of  Data  and  Models 

A  large  number  of  different  sources  of  data  were  used 
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in  the  analyses  presented  in  this  dissertation.  Daily 
precipitation,  maximum  and  minimum  daily  temperatures ,  dew 
point  temperatures ,  and  DAD  analyses  were  obtained  from  the 
numerous  publications  and  data  tapes  of  the  Atmospheric 
Environment  Service  of  Canada,  Downsview,  Ontario.  Daily 
streamflow  discharges  and  stream  and  river  characteristics 
came  from  publications  and  data  sources  produced  and 
collected  by  the  Water  Survey  of  Canada,  Inland  Waters 
Directorate.  River  and  hydrological  char acter i s t i cs  were 
also  extracted  from  numerous  publications,  and  these  are 
listed  in  the  bibliography  and  are  referenced  where 
applicable.  The  SSARR  model  with  appropriate  write-up  and 
descriptions  of  usage  were  obtained  from  Alberta 
Environment.  The  HYMO  model  was  obtained  from  the  City  of 
Edmonton,  Sanitary  Engineering  Department.  The  author  would 
like  to  acknowledge  the  contributions  of  data,  models  and 
information  from  all  of  the  above  sources,  for  without  this 
assistance  the  results  in  this  dissertation  together  with 
the  results  presented  in  the  author's  publications 
(Verschuren  and  Wojtiw,  1980;  and  Wojtiw  and  Verschuren, 
1981)  would  not  have  been  possible. 
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CHAPTER  2 


CHARACTERISTICS  OF  PMP 


Probable  Maximum  Precipitation  can  be  estimated  by  a 
number  of  methods,  the  two  internationally  accepted  methods 
(World  Meteorological  Organization,  1973)  are  the  physical 
(also  called  meteorological  or  traditional)  and  the 
statistical.  Both  of  these  methods  have  been  examined 
extensively  by  Verschuren  and  Wojtiw  (1980)  in  their 
publication  dealing  with  Alberta  river  basins.  In  order  not 
to  duplicate  this  effort,  only  some  of  the  results  and 
figures  are  discussed  and  presented  in  this  dissertation. 

The  reader  should  refer  to  the  above  publication  for  an  in- 
depth  look  at  this  topic.  The  results  presented  here  deal 
only  with  the  Red  Deer  River  Basin,  which  is  used  as  an 
example  throughout  this  work.  It  should  be  noted  that  the 
sequence  of  major  storm  rainfalls  was  examined  but  not 
considered  to  be  involved  in  PMP  estimation. 

2.1.0  Estimation  of  PMP  by  the  Physical  Method 

The  physical  approach  to  obtaining  an  estimate  of  PMP 
was  developed  over  a  number  of  years.  The  method  is  indirect 
and  is  based  on  the  analysis  and  maximization  of  the  largest 
rainstorms  that  occurred  during  the  period  for  which  there 
are  rainfall  records.  It  assumes  that  the  amount  of  rainfall 


53 


■ 

. 


54 


from  a  storm  depends  on  two  independent  factors:  1)  the 
moisture  content  of  the  air  mass  and  2)  the  efficiency  of 
the  rain-producing  system.  The  procedure  used  for  maximizing 
observed  storm  rainfall  to  estimate  PMP  involves  moisture 
adjustments,  storm  transposition,  and  envelopment  1  .  It 
further  assumes  that  the  largest  rainstorm  in  the  region  has 
occurred  during  the  period  of  record.  This  is  an  important 
assumption,  which  is  open  to  criticism,  especially  regarding 
regions  where  scanty  meteorological  information  is 
available.  Because  of  its  meteorological  components, 
however,  the  method  is  much  easier  to  relate  to  the  physical 
occurrence  and  effect  of  rainfall  from  rainstorms  and  hence 
has  been  readily  accepted  by  the  meteorological  and 
hydrological  communities. 

2.1.1  Estimation  of  Atmospheric  Moisture 

In  production  of  precipitation,  the  moisture  in  the 
lower  layers  of  the  atmosphere  is  most  important  (Schwarz, 
1967;  U.S.  Department  of  Commerce,  1960).  Theoretical 
computations  show  that  for  excessive  rain,  air  originally 
near  the  surface  reaches  within  an  hour  or  so  the  top  of  the 
layer  from  which  the  precipitation  is  falling.  In  the  case 
of  severe  thunderstorms ,  surface  air  may  reach  the  top  of 
the  storm  in  a  matter  of  minutes  because  of  the  high- 


1  A  process  in  which  the  largest  value  of  PMP  is  selected 
from  a  set  of  data  for  a  given  area,  and  a  smooth  curve 
drawn  is  through  the  largest  values. 
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velocity  updrafts  in  such  storms. 

The  most  realistic  assumption  seems  to  be  that  the  air 
ascends  dry-adiabatical ly  to  the  saturation  level  and  then 
ascends  moi st-adi abatical ly .  Thus  hydrometeorologists 
generally  postulate  a  saturated  pseudoadi abat i c  atmosphere 
for  extreme  storms.  To  maximize  the  moisture  of  a  storm,  two 
saturation  adiabats  are  required.  The  first  gives  a  measure 
of  the  vertical  temperature  distribution  in  the  storm  to  be 
maximized,  and  the  second  is  the  warmest  saturation  adiabat 
to  be  expected  at  the  same  location  and  time  of  year  as  the 
storm.  In  meteorology  these  two  saturation  adiabats  are 
identified  by  the  wet-bulb  potential  temperature,  which 
corresponds  with  the  dew  point  temperature  at  the  1000  mb 
level.  The  work  by  the  U.S.  Department  of  Commerce  (1960) 
has  shown  that  storm  and  extreme  values  of  precipi table 
water  may  be  approximated  by  estimates  based  on  surface  dew 
points  when  saturation  and  pseudoadi abat i c  conditions  are 
assumed.  An  approximation  based  on  surface  dew  point 
measurements  is  desirable  since  these  measurements  are  more 
readily  available  than  measurements  at  the  1000  mb  level. 
Hence  surface  dew  points  are  used  in  identifying  the  storm 
saturation  adiabat  for  maximizing  the  moisture  content  of 
the  storm.  Both  storm  and  maximum  dew  points  can  be  reduced 
pseudoadi abat i ca 1 1 y  to  the  1000  mb  level  by  use  of  Figure 
2.1,  so  that  dew  points  observed  at  stations  of  different 
elevations  are  comparable. 

Since  the  moisture  has  an  appreciable  effect  on  the 
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FIGURE  2.1.  Pseudoadi abat i c  Diagram  for  Dew  Point  Reduction 
to  1000  mb  at  Height  Zero  (After  World  Meteorological 

Organization,  1973). 

storm,  precipitation  must  be  that  which  persists  for  hours 
rather  than  minutes.  Also,  any  single  observation  of  dew 
point  may  be  considerably  in  error;  therefore,  for 
estimating  storm  and  probable  maximum  moisture  the 
conventional  procedure  is  to  use  dew  point  values  on  two  or 
more  consecutive  measurements  separated  by  a  reasonable  time 
interval.  The  adopted  procedure  is  to  use  the  so-called 
highest  persisting  12-hour  dew  point.  This  is  the  highest 
value  equalled  or  exceeded  by  all  observations  during  a  12- 
hour  period. 

Paulhus  and  Gilman  (1953)  concluded  that  the  12-hour 
persisting  dew  point  temperature  was  most  representative  of 
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the  average  moisture  in  the  storm  center.  The  results  of 
this  study  have  since  been  used  and  accepted  internationally 
in  estimating  PMP  by  the  physical  method.  This  is  the  method 
used  in  this  work. 

It  is  from  the  maximum  persisting  12-hour,  1000  mb  dew 
point  temperatures  that  the  maximum  values  of  atmospheric 
water  vapor  used  for  storm  maximization  are  estimated.  These 
dew  points  are  obtained  for  stations  in  a  given  river  basin. 
Since  numerous  estimates  of  PMP  are  required,  each  station 
must  have  an  enveloping  curve  of  semimonthly  recorded 
maximum  persisting  12-hour  dew  points.  The  author  has 
produced  such  curves  for  Alberta  (Verschuren  and  Wojtiw, 
1980),  with  an  example  for  the  city  of  Red  Deer  given  in 
Figure  2.2;  graphs  for  the  other  stations  can  be  found  in 
the  above  publication.  The  16  graphs  produced  in  that  report 
provide  a  comprehensive  source  of  analyzed  information  for 
Alberta.  Other  studies  had  used  only  a  few  stations  with 
less  than  half  the  record  period,  because  of  the  limited 
data  base  available.  With  the  aid  of  these  curves,  moisture 
adjustments  can  be  made  on  the  basis  of  the  maximum 
persisting  12-hour  dew  point  for  the  same  time  of  year  as 
the  storm  occurrence.  Thus,  for  example,  a  May  maximum  dew 
point  (and  not  one  in  September,  for  instance)  would  be  used 
to  maximize  a  May  storm.  Also  prepared  were  monthly  maps  of 
maximum  persisting  12-hour,  1000  mb  dew  points,  which  not 
only  served  as  a  convenient  source  of  maximum  dew  points  but 
aided  in  maintaining  consistency  between  estimates  for 
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RED  DEER 

NRXJHUH  12-HOUR  PER81ST1N0 
OEM  POINT  TEMP  RT  1000  ItB 
IN  16-DRY  PERIODS 


WRY  J UN  JUL 

MONTHS  OF  THE  TERR 


FIGURE  2.2.  Maximum  Recorded  12-hour  Persisting  Dew  Point 
Temperatures  at  1000  mb  for  the  city  of  Red  Deer  (Verschuren 

and  Wo j t i w ,  1980 ) . 

various  basins.  The  maps  were  prepared  by  using  the  monthly 
maximum  12-hour  dew  point  values,  adjusting  them  to  the  1000 
mb  level,  plotting  them  at  the  locations  of  the  observing 
stations,  and  drawing  smooth  isopleths.  An  example  of  such  a 
map  is  shown  in  Figure  2.3  for  the  month  of  June  in  Alberta. 
The  maps  for  other  months  can  be  found  in  Verschuren  and 
Woj  t i w  (1980). 

Before  an  estimate  of  the  moisture  maximization  is 
obtained,  it  is  necessary  to  find  the  amount  of  precipi table 
water  available  from  the  storm  as  well  as  the  potential 
maximum  amount  produced  at  the  maximum  location.  The  amount 
of  precipi table  water,  W  (cm),  can  be  computed  by  the 
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FIGURE  2.3.  Isopleths  of  Maximum  12-hour  Persisting  Dew 
Point  Temperatures  for  June  (Verschuren  and  Wojtiw,  1980). 


following  general  formula: 

q  AP 

w  =  - 

gf 


(2.1) 


_  -1 
where  q  is  the  mean  specific  humidity  in  gm  Kg  of  a  layer 
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of  moist  air.Ap  is  the  depth  of  the  layer  in  mb,  g  is  the 

_2 

acceleration  of  gravity  in  cm  sec  ,  and  j?  is  the  density 

_  j 

of  water  (equal  to  1  gm  cm  ).  For  convenience,  Equation  2.1 
can  be  pre  computed  and  is  usually  listed  in  table  or 
monogram  form.  Two  tables  are  given  in  APPENDIX  I:  the 
first,  Table  A.I.1,  presents  values  of  precipitable  water 
(mm)  between  the  1000  mb  surface  and  various  pressure  levels 
up  to  300  mb  in  a  saturated  pseudoadi abat ic  atmosphere  as  a 
function  of  the  1000  mb  dew  point.  The  second  table,  Table 
A. I. 2,  presents  values  of  precipitable  water  (mm)  for  layers 
between  the  1000  mb  surface,  assumed  to  be  at  zero 
elevation,  and  various  heights  up  to  8  Km  (5  mi.).  The  300 
mb  level  is  generally  accepted  as  the  top  of  the  storm,  but 
it  makes  little  difference  which  level  from  400  mb  on  up  is 
selected,  as  there  is  very  little  moisture  at  those  heights, 
and  the  effect  on  the  moisture  adjustment  is  negligible. 

2.1.2  Moisture  Maximization 

Moisture  maximization  of  rainstorms  in  place  2  •, 
consists  of  multiplying  the  observed  storm  rainfall  amounts 
by  the  ratio  (rm)  of  the  maximum  precipitable  water  (Wm  ) 
indicated  for  the  storm  location  over  the  precipitable  water 
(Ws  )  estimated  for  the  storm  (World  Meteorological 
Organization,  1973): 


2  Without  change  in  location.  Storm  transposition  requires  a 
further  moisture  adjustment. 
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(2.2) 


2.1.3  Storm  Transposition 

Storm  transposition  refers  to  the  transfer  of  storms 
from  locations  where  they  occurred  to  other  areas  where  they 
could  occur.  The  transposition  procedure  involves  the 
meteorological  analysis  of  the  storm  to  be  transposed,  the 
determination  of  the  limits  of  transposabi 1 i ty ,  and  the 
application  of  the  proper  adjustments  for  making  the 
modifications  required  by  the  change  in  storm  location.  Two 
schools  of  thought  seem  to  prevail  on  this  topic  in  the 
hydrological  community.  The  first  advocates  extensive 
transposition  of  storms  over  hundreds  or  thousands  of 
kilometres.  The  greatest  danger  in  this  approach  is  deciding 
to  transpose  a  storm  when  only  limited  analysis  is  planned 
or  attempted.  The  second  school  of  thought  advocates  limited 
transposition  no  more  than  a  few  hundred  kilometres.  The 
argument  is  that  if  no  storms  have  occurred  in  the  region  in 
the  last  50  years  or  so,  they  probably  will  not  occur,  since 
the  topography  is  not  conducive  to  such  occurrences.  The 
author  in  this  study  advocates  the  latter  approach  and  hence 
has  conducted  limited  storm  transposition. 

The  storms  that  were  transposed  were  adjusted  according 
to  the  geographic  features  of  that  particular  basin.  The 
limits  of  transposabi 1 i ty  were  governed  not  only  by  the 
constant  lines  of  maximum  persisting  12-hour,  1000  mb  dew 
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point  temperatures  but  also  by  the  meteorology 
storms.  After  storm  transportation  the  results 
are  usually  depicted  in  a  depth-area  curve;  an 
such  a  curve  is  shown  in  Figure  2.4. 
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FIGURE  2.4.  Depth-Area  Curves  of  Maximized  24-hour 
the  Red  Deer  River  Basin  (Verschuren  and  Wojtiw, 


PMP  for 
1980)  . 


A  single  storm  maximized  and  transposed  to  a  basin 
gives  the  precipitation  depth  that  could  fall  from  that 
particular  storm  over  that  basin.  In  employing  a  single 
storm,  there  is  no  guarantee  that  the  maximum  magnitude  of 
the  precipitation  has  been  achieved,  since  no  single  storm 
is  likely  to  yield  extreme  rainfall  values  for  all  durations 
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and  areal  sizes.  To  obtain  this  maximum,  all  storms  or  as 
many  as  are  available  are  plotted  for  the  various  durations 
and  sizes  of  area.  Some  storms  may  not  contribute  at  all  to 
the  enveloping  curve,  while  others  may  contribute  at  certain 
durations.  An  example  of  a  number  of  storms  plotted  on  a 
depth-area  curve  of  maximized  24-hour  PMP  for  the  Red  Deer 
River  Basin  is  shown  in  Figure  2.4.  In  the  final  analysis, 
the  maximum  values  are  dominated  by  only  a  few  severe 
rainstorms ;  however ,  because  it  is  difficult  to  determine 
whether  the  values  of  a  particular  rainstorm  will  be  maximum 
in  duration  and/or  area  for  the  basin,  it  is  necessary  to 
plot  all  the  most  logical  available  rainfall  analyses. 
Verschuren  and  Wojtiw  (1980),  wrote  a  computer  program  for 
rapidly  producing  such  curves  and  included  various 
combinations  of  durations  for  the  six  major  river  basins  in 
Alberta.  From  the  maximized  rainfall  data  the  largest  value 
from  any  set  of  data  is  selected,  and  a  smooth  curve  is 
drawn  through  the  largest  values  (envelopment).  A  number  of 
storms  are  plotted  in  Figure  2.5,  where  enveloping  DAD 
curves  for  each  specific  duration  are  drawn  for  the  Red  Deer 
River  Basin. 

2.2.0  Estimation  of  PMP  by  the  Statistical  Method 

Statistical  procedures  for  estimating  PMP  may  be  best 
used  wherever  sufficient  precipitation  data  are  available. 

The  method  is  useful  for  making  estimates  or  when  other 
meteorological  data,  such  as  dew  point  and  wind  records, 
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FIGURE  2.5.  Enveloping  Dep th-Area-Durat ion  Curves  of  PMP  for 
the  Red  Deer  River  Basin  (Verschuren  and  Wojtiw,  1980). 

lacking.  The  procedure  was  developed  by  Hershfield  (1961, 
1965,  and  1977),  and  though  it  is  not  the  only  approach,  it 
is  one  that  has  received  the  widest  acceptance,  perhaps 
because  it  requires  considerably  less  time  to  apply  than 
does  the  physical  approach,  and  needs  only  a  limited 
understanding  of  meteorology. 


2.2.1  Statistical  Procedure 

The  statistical  procedure,  as  developed  by  Hershfield 
and  introducted  in  Section  1.2.3,  is  based  on  the  general 
frequency  equation  (Equation  1.1)  given  by  Chow  (1951), 
which  can  be  rewritten  as 


(2.3) 
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where  Xmis  the  maximum  observed  rainfall,  X  n  and  S  n  are  the 
mean  and  standard  deviation  of  a  series  of  n  annual 
precipitation  maxima  for  a  specific  duration,  and  Kmis  the 
frequency  factor  equivalent  to  the  number  of  standard 
deviations  to  be  added  to  Xn  to  obtain  Xm. 

To  obtain  a  value  for  variable  Km,  Hershfield  first 
computed  for  each  station  the  mean  and  standard  deviation 
using  the  conventional  precedures;  however,  in  the 
computations  he  omitted  the  maximum  observed  rainfall.  Then 
in  Equation  2.3  he  substituted  these  values  for  Xn  and  Snin 
as  well  as  the  maximum  observed  rainfall  for  Xm  to  solve  for 
Km.  This  is  equivalent  to  observing  the  maximum  event  after 
the  values  of  the  basic  statistics  have  been  established. 

The  frequency  factor  has  been  observed  to  be  a  function  of 
duration  and  the  mean  annual  extreme  rainfall  (Hershfield, 
1977).  This  relationship  was  given  by  Equations  1.2,  1.3, 
and  1.4. 

Before  PMP  can  be  computed  by  the  statistical 
procedure,  a  number  of  adjustments  are  required  to  the  mean 
and  standard  deviation  for  maximum  observed  event  and  sample 
size . 

2.2.2  Adjustment  of  X n  and  Snfor  Maximum  Observed  Event 

Hershfield  (1961)  states  that  extreme  rainfall  amounts 
of  rare  magnitude  (i.e.,  with  return  periods  of  500  or  more 
years)  are  believed  to  have  occurred  during  a  much  shorter 
period  of  record,  for  example,  30  years.  Such  rare  events, 
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or  outliers,  may  have  an  appreciable  effect  on  the  mean  (Xn) 
and  standard  deviation  (SR)  of  the  annual  series.  The 
magnitude  of  the  effect  is  less  for  long  records  than  for 
short,  and  it  varies  with  the  rarity  of  the  event.  This 
effect  has  been  studied  by  Hershfield  (1961)  using 
hypothetical  series  of  varying  lengths  and  can  be  summarized 
by  the  relationship  of  the  maximum  observed  event  to  the 


FIGURE  2.6.  Adjustment  of  Mean  of  Annual  Series  for  Maximum 
Observed  Rainfall  (After  Hershfield,  1961). 
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mean,  as  shown  in  Figure  2.6.  Similiarly,  Figure  2.7  (after 


FIGURE  2.7.  Adjustment  of  Standard  Deviation  of  Annual 
Series  for  Maximum  Observed  Rainfall  (After  Hershfield, 

1961  )  . 


Hershfield)  shows  the  relationship  of  maximum  observed 
events  on  the  standard  deviation.  In  these  fiqures  X  and 
Sn_m  refer  to  the  mean  and  standard  deviation  of  the  annual 
series  computed  after  excluding  the  maximum  item  in  the 
series.  In  both  diagrams,  the  relationships  are  for  only  the 
effect  of  the  maximum  observed  event;  no  consideration  was 
given  to  other  anomalous-appearing  observations. 
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2.2.3  Adjustment  of  Xn  and  Sn  for  Sample  Size 

The  mean  ( X  of  the  annual  series  may  tend  to  increase 
with  record  length,  because  the  frequency  distribution  of 
rainfall  extremes  is  skewed  to  the  right  so  that  there  is  a 
greater  chance  of  getting  a  larger  rather  than  a  smaller 
extreme  as  record  length  increases.  Using  data  from  198  key 
weather  stations  and  adjusting  for  an  outlier,  Hershfield 
(1961)  determined  the  average  ratios  of  a  50-year  mean  to 
10- ,  15-,  20- ,  and  30-year  means;  the  adjustment  necessary 
for  the  length  of  record  is  shown  in  Figure  2.8.  In  his 
computations,  the  statistics  from  the  50-year  records  were 
used  as  a  standard,  so  that  adjustments  were  required  for 
those  with  shorter  record  periods.  A  comparison  of  a  small 
number  of  available  means  for  records  longer  than  50  years 
with  the  50-year  means  showed  only  a  negligible  difference. 
Similiarly,  the  necessary  standard  deviation  adjustment  for 
length  of  record  was  computed,  and  the  results  are  shown  in 
Figure  2.8.  The  effect  of  the  record  length  is  much  more 
pronounced  on  the  standard  deviation  than  on  the  mean.  The 
few  records  for  longer  than  50  years  indicate  the  need  for 
only  a  slight  adjustment  from  the  values  for  the  50-year 
records . 

2.2.4  Adjustment  for  Fixed  Observational  Time  Intervals 
Precipitation  data  are  usually  recorded  on  a  fixed  time 

interval,  e.g.,  hourly,  six-hourly,  or  daily.  Such  data 
rarely  yield  the  true  maximum  rainfall  amounts  for  the 
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FIGURE  2,8.  Adjustment  of  Mean  and  Standard  Deviation  of 
Annual  Series  for  Length  of  Record  (After  Hershfield,  1961). 

indicated  durations.  As  an  example,  the  annual  maximum 
observational  day  amount  is  very  likely  to  be  appreciably 
less  than  the  annual  maximum  24-hour  amount  determined  from 
intervals  of  1440  consecutive  minutes  unrestricted  by  any 
particular  observation  time.  Studies  by  Weiss  (1964) 
indicate  that  when  the  results  of  a  frequency  analysis  of 
annual  maximum  rainfall  amounts  for  a  single,  fixed-time 
interval  of  any  duration  from  1  to  24  hours  are  multiplied 
by  1.13  they  will  yield  values  closely  approximating  those 


. 
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obtained  from  an  analysis  based  on  true  maxima.  Hence,  the 
PMP  values  yielded  by  the  statistical  procedure  should  be 
multiplied  by  1.13  if  data  for  single,  fixed-time  intervals 
are  used  in  compiling  the  annual  series.  Figure  2.9  shows 
the  smaller  adjustments  necessary  when  durations  are 
determined  from  two  or  more  fixed- time  intervals.  As  an 


FIGURE  2.9.  Average  Adjustment  of  Fixed  Interval 
Precipitation  Amounts  for  Number  of  Observational  Units 
within  the  Interval  (After  Weiss,  1964). 


example,  maximum  24-hour  amounts  determined  from  24 
consecutive,  1-hour  rainfall  increments  requires  an 
adjustment  by  a  factor  of  1.01.  The  results  from  these 
studies  show  the  importance  to  the  user  of  understanding  the 
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time  collection  of  precipitation  data. 

2.2.5  Area-Reduction  Curves 

One  of  the  shortcomings  of  the  statistical  procedure  is 
that  it  yields  only  point  values  of  PMP  and  thus  requires 
area-reduction  curves  for  adjusting  the  point  values  to 
various  sizes  of  area.  A  number  of  variations  of  depth-area 
relationship  exist  (Court,  1961).  For  the  prairie  provinces, 
McKay  (1965b)  developed  the  area-reduct ion  curves  shown  in 
Figure  2.10.  These  curves  relate  point  values  to  areas  up  to 
129  500  Km2  (50  000  sq.  mi.).  For  rainfall,  point  values  are 
often  assumed  to  be  applicable  to  areas  up  to  65  Km2  (25  sq . 
mi.)  without  reduction.  These  reductions  can  be  called  areal 
in  nature,  in  the  sense  that  once  a  maximum  point  value  is 
defined  spatially  it  can  be  converted  to  an  estimate  for  a 
larger  area.  These  are  called  reductions  because  the  maximum 
point  value  has  been  observed  to  decrease  with  an  increase 
in  area  of  representat ion . 

2.2.6  Statistical  Estimates  for  Alberta  River  Basins 

Employing  the  various  concepts  described  above, 

Verschuren  and  Wojtiw  (1980)  obtained  statistical  PMP 
estimates  for  the  major  river  basins  in  Albera.  For  each 
recorded  year,  maximum  24-,  48-,  72-,  and  96-hour  rainfalls 
were  calculated  for  27  first-order  weather  stations  in 
Alberta.  Then  the  statistical  estimates  were  computed  using 
the  procedure  developed  by  Hershfield.  An  example  of  the  24- 
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FIGURE  2.10.  Depth-Area,  or  Area-Reduction,  Curves  (After 

McKay ,  1 965b ) . 

hour  enveloping  depth-area  curve  for  the  Red  Deer  station  is 
shown  in  Figure  2.11.  In  general,  for  the  six  river  basins 
the  statistical  approach  produced  results  that  were  slightly 
higher  than  those  calculated  by  the  physical  approach. 

2.2.7  Spatial  Characteristics  of  PMPs 

The  results  seem  to  suggest  that  in  Alberta  the  PMP 


73 


»  o  100  200  300  400  500  600  700 


FIGURE  2.11.  24-hour  Enveloping  Depth-Area  Curves  for  the 
Red  Deer  River  Basin  (Verschuren  and  Wojtiw,  1980). 

varies  spatially  in  each  basin,  as  can  be  seen  from  the  24- 
hour  estimates  in  Figure  2.12.  The  largest  PMP  estimates  are 
observed  in  the  South  Saskatchewan  River  Basin,  while  the 
smallest  are  in  the  eastern  and  northern  portions  of  the 
province.  Generally,  two  characteristic  variations  can  be 
defined:  the  first,  also  being  the  most  predominant,  is  the 
variation  from  west  to  east  in  the  basin,  while  the  second 
is  the  variation  northward.  The  west-east  variation  exibits 
two  main  features.  First,  there  is  a  strong  decreasing 
gradient  in  the  southern  parts  of  the  province  (namely,  the 
South  Saskatchewan  River  Basin).  PMP  values  decrease  from 
over  500  mm  (20  in.)  in  the  western  portion  to  about  100  mm 
(4  in.)  in  the  eastern  portion  of  the  river  basin.  The 
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FIGURE  2.12,  Spatial  Variation  (mm)  of  the  Maximized  24-hour 
PMP  (Verschuren  and  Wojtiw,  1980). 

second  feature  is  an  increase  (about  50%)  in  the  central  and 
northern  river  basins.  The  variation  northward  is  important 
mainly  in  the  southern  basins  of  the  province,  where  changes 
of  about  40%  are  not  uncommon.  For  basins  in  central  and 
northern  Alberta  only  slight  variations  are  observed  and 
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these  are  of  very  little  consequence  in  the  PMP  estimate. 

The  coefficient  of  variation,  which  is  the  ratio  of  the 
standard  deviation  to  the  mean,  has  been  used  (McKay,  1965b) 
to  obtain  estimates  of  the  climatological  day  24-hour 
precipitation  extreme  of  given  frequency.  To  examine  this 
aspect,  extreme  value  series  of  24-hour  precipitation  were 
used  from  all  available  long-term  weather  stations  in 
Alberta.  The  isolines  of  the  coefficient  of  variation 
(expressed  in  percentage)  for  the  annual  24-hour 
precipitation  extremes  are  depicted  in  Figure  2.13. 


FIGURE  2.13  Coefficients  of  Variation  of  Annual 

Precipitation  Extremes. 


24-hour 
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CHAPTER  3 

SNOWMELT  WITH  PMP  LOADING 

3.1.0  Introduct ion 

In  Alberta,  rainstorms  usually  occur  from  April  to 
September,  with  the  majority  of  these  occurring  during  June 
and  July  (Verschuren  and  Wojtiw,  1980).  Of  these  rainstorms, 
the  rain-on-snow  event  is  believed  by  many  authors  (McKay, 
1965b,  1968;  Buckler  and  Quine,  1971)  to  have  the  best 
potential  for  producing  the  maximum  water  loading  on  the 
watersheds;  this  usually  occurs  in  late  May  or  early  June  in 
Alberta  (Verschuren  and  Wojtiw,  1980).  The  few  studies  that 
have  examined  this  event  have  used  the  generalized  snowmelt 
equations  and  coefficients  developed  by  the  U.S.  Army  Corps 
of  Engineers  (1956)  for  basins  in  United  States,  without 
examining  their  validity  for  Alberta  basins.  An  objective  in 
the  present  study  was  to  examine  the  rain-on-snow  event  and 
determine  the  contribution  of  snowmelt  to  PMF .  It  was  also 
useful  to  define  the  relevant  processes  contributing  to 
snowmelt  for  the  PMP  loading.  This  chapter  discusses  these 
objectives,  together  with  the  development  of  an  equation  for 
computing  snowmelt  for  Alberta  river  basins.  The  snowmelt 
equation  presented  can  be  expressed  by  meteorological 
parameters;  which  in  turn  can  be  expressed  in  terms  of 
location,  elevation,  and  drainage  area.  A  number  of 
relationships  between  these  latter  parameters  can  be 


. 
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developed  to  provide  spatial  estimates  of  the  meteorological 
parameters  (Section  3.3).  These  relationships  can  then  be 
used  to  examine  the  snowmelt  for  spatial  variability  and 
maximum  melt  rate,  and  to  obtain  estimates  for  watersheds 
lacking  meteorological  data. 

3.2.0  Rain-on-Snow  Event 

The  generation  of  snowmelt  at  a  point  location  in  a 
snowpack  is  essentially  a  thermodynamic  process,  with  the 
amount  of  melt  produced  being  dependent  on  the  net  heat 
exchange  between  the  snowpack  and  its  environment.  To 
represent  physical  processes  that  are  not  easily  measured, 
indexes  expressed  in  terms  of  readily  observed  measurements 
between  the  snowpack  and  its  environment  are  usually  used. 
The  reliability  of  the  indexes  depends  upon  (a)  how  well  the 
physical  process  is  described  by  the  measured  variable,  (b) 
the  random  variability  of  the  measurement ,  (c)  the 

variability  between  point  values  and  areal  values,  and  (d) 
the  validity  of  the  measured  value.  The  index  relationship 
may  be  described  by  a  coefficient  and  may  be  either  constant 
or  variable,  depending  upon  the  variability  of  associated 
factors . 

The  atmospheric  temperature  is  a  useful  parameter  in 
snowmelt  determination  and  has  been  used  as  an  index, 
because  it  is  generally  thought  to  be  the  best  index  of  the 
heat  transfer  processes  associated  with  snowmelt  and  is  a 
reliable  and  regularly  available  meteorological  variable. 
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Snowmelt  is  the  overall  result  of  many  different 
processes  of  heat  transfer:  shortwave  and  longwave 
radiation,  condensation  of  vapor,  convection,  air  and  ground 
conduction,  and  rainfall.  What  follows  is  an  examination  of 
each  of  these  melt  contributions  in  terms  of  PMP  loading. 


3.2.1  Radiation  Melt 

The  amount  of  incoming  solar  (shortwave)  radiation 
received  by  a  snowpacK  can  be  a  very  important  source  of 
heat  energy  for  snowmelt.  The  net  amount  of  heat  absorbed  by 
the  snowpacK  depends  on  latitude,  orientation  and  slope, 
season,  time  of  day,  atmospheric  conditions  (clouds,  fog, 
rain),  forest  cover,  and  reflectivity  of  the  snow  (albedo). 
Of  these,  cloud  cover  and  albedo  have  been  identified  as 
being  significant  to  snowmelt.  The  intensity  of  insolation 
or  shortwave  solar  radiation  is  normally  expressed  in 
langleys  ( ly)  (calories  per  square  centimetre)  per  unit  of 
time  (minute,  hour,  or  day)  and  abbreviated,  for  example,  as 
ly/hr.  The  most  simple  and  accurate  measurement  of  net 
insolation  uses  pyrhe 1 i ometer s ;  these  instruments  are 
calibrated  to  give  the  shortwave  radiation  intensity  in 
ly/hr  or  ly/day. 

The  snowmelt  by  shortwave  radiation,  M  ,  can  be 

rs 

approximated  by  (U.S.  Army  Corps  of  Engineers,  1956) 
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where  M  is  the  snowmelt  in  inches  due  to  shortwave 

X  o 

radiation,  a  is  the  albedo  written  as  a  decimal  fraction, 

Rsi  is  the  effective  solar  radiation  in  ly/day,  and  B  is  the 

thermal  quality  of  the  snowpack  (defined  as  the  ratio  of 

heat  required  to  melt  a  unit  weight  of  snow  to  that  of  ice 

at  0°C  or  32°F).  The  thermal  quality  averages  from  0.97  to 

0.95,  for  a  snowpack  with  3-5%  liquid  water  and  is  usually 

assumed  to  be  0.97  for  a  melting  snowpack.  The  constant, 

203.2,  is  the  heat  input  in  langleys  required  to  produce  one 

o 

inch  of  water  from  ice  at  32  F ;  it  is  in  term  of  ly/in.. 

Since  the  snowmelt  depends  upon  the  air  temperature 
over  the  snow  and  the  temperature  of  the  cloud  base,  the 
maximum  snowmelt  rate  in  most  cases  would  probably  occur 
under  clear  skies  as  compared  to  under  skies  with  a  complete 
or  partial  cloud  cover.  The  albedo  is  commonly  considered 
80%  for  fresh  snow  and  is  assumed  to  decrease  exponentially 
to  about  40%  for  melting,  late-season  snow  (Figure  3.1). 
Equation  3.1  indicates  that  the  melt  increases  with  a 
decrease  in  albedo  and  that  the  greatest  melt  would  be 
produced  when  the  smallest  value  of  albedo  is  used. 
Insolation  during  rain  conditions  is  commonly  assumed  to  be 
close  to  40  ly/day  for  an  open  area,  and  an  albedo  of 
roughly  65%  gives  a  contribution  of  M  equal  to  1.78  mm  per 
day  (0.07  in.  per  day)  using  Equation  3.1  (Gray,  1970).  The 
albedo  value  of  65%  has  been  shown  to  be  much  too  high  for 
Alberta  (especially  from  April  to  June)  by  Hay  (1970).  Hay 
estimated  regionally  representative  monthly  mean  albedo 
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Figure  3.1.  Variation  in  Snow  Surface  Albedo  with  Time 
(After  U.S.  Army  Corps  of  Engineers,  1956). 

values  for  selected  Alberta  locations.  Mean  values  are  used 
because  rainstorms  that  may  produce  precipitation 
approaching  the  PMP  are  on  a  scale  of  days  as  compared  to 
minutes  or  a  few  hours.  For  June,  Hay  gives  mean  values  of 
20%  for  Beaver  lodge  and  23%  for  Edmonton  and  Suf field.  About 
20%  seems  as  a  reasonable  value  to  use  as  a  minimum  value 
for  Alberta,  since  there  is  little  variation  in  the  values 
between  the  northern  and  the  southern  stations.  Although 
zero  in  the  smallest  value  that  can  be  achieved  under 
specific  conditions,  it  is  not  a  practical  value  since  it  is 
rarely  recorded  in  nature  during  precipitation  events.  The 
relationship  (if  one  exists)  between  albedo  and 
precipitation  events  is  not  well  understood. 
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A  maximum  snowmelt  rate  can  be  expected  in  late  May  or 
early  June,  since  the  coincidence  of  severe  rainstorms 
(Verschuren  and  Wojtiw,  1980),  snow  cover  (Figure  3.2),  and 
air  temperatures  above  freezing  seem  to  occur  in  Alberta 
during  this  period.  The  amount  of  radiation  reaching  the 


Figure  3.2.  Date  of  last  snow  cover  of  25  mm  ( 1  in.)  or  more 
for  Alberta  (Verschuren  and  Wojtiw,  1980). 


surface  is  usually  expressed  in  terms  of  the  average 
cloudless  day's  insolation;  in  Alberta  these  values  increase 
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from  May  to  June  (Figure  3.3).  Average  values  between  760 
and  about  800  ly/day  seem  to  be  representat i ve  of  prairie 
basins  in  Alberta  (Gray,  1970),  with  the  large  insolation 
values  being  observed  at  the  head  of  the  basins  (that  is, 
near  the  mountains)  while  the  lower  values  occur  at  the 
mouth  of  the  basins  (usually  in  the  plains).  Hence  a  value 
of  800  ly/day  is  used  for  prairie  basins,  compared  to  40 
ly/day  given  in  Gray  (1970).  The  value  of  40  ly/day  may  be 
an  interpretat ion  error.  Exactly  the  same  number  appears  in 
the  U.S.  Army  Corps  of  Engineers  report  (1956,  page  196); 
however,  there  it  refers  to  the  standard  error  of  the 
estimated  insolation.  Substituting  a  value  of  800  ly/day 
into  Equation  3.1,  the  melt  component  produced  by  shortwave 
radiation  is  82.6  mm  per  day  (3.25  in.  per  day)  or  3.4  mm 
per  hour  (0.14  in.  per  hour). 

Snow  is  considered  to  be  a  nearly  perfect  black  body 
with  respect  to  longwave  radiation.  The  longwave  radiation 
emitted  by  a  snow  surface  can  be  estimated  using  Stefan's 
law  (total  radiation  is  proportional  to  the  fourth  power  of 
temperature)  (Short  ley  and  Williams,  1961;  U.S.  Army  Corps 
of  Engineers,  1956).  Computation  of  net  longwave  radiation 
involves  estimation  of  back  radiation  from  the  atmosphere 
under  clear  or  cloudy  skies  and  from  beneath  forest  cover, 
and  these  evaluations  are  difficult  to  calculate  and  usually 
not  used  in  practical  snow  hydrology.  Instead,  U.S.  Army 
Corps  of  Engineers  (1960)  recommends  the  use  of  simplified 
linear  relationships  in  terms  of  temperature.  Longwave 
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Figure  3.3.  Average  Cloudless  Day  Insolat 
May  15  and  dune  15  for  Canada  (After 


on  in  Langleys  for 
Gray ,  1970). 


radiation  melt 


for 


rain  condi t ions 


can  be  considered 


to  be 


under  complete  cloud  cover;  also,  the  cloud  base  temperature 
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is  approximately  equal  to  air  temperature.  The  relationship 
between  longwave  radiation  and  air  temperature  was  found 
experimentally  to  be  linear  (U.S.  Army  Corps  of  Engineers, 
1956)  and  can  be  expressed  as 

Mrl  =  °'029  <Tm  -  32)  (3. 2) 

where  Mrl  is  the  daily  melt  in  inches  resulting  from 

longwave  radiation  exchange  and  is  the  mean  of  the  daily 

o 

maximum  air  temperatures  in  F  during  the  specific  month. 

There  is  no  indication  in  the  literature  whether  the 

coefficient  in  Equation  3.2  is  dependent  upon  parameters 

other  than  the  sky  conditions  and  the  canopy.  Hence  it  is 

probably  reasonable  to  assume  that  the  coefficient  can  be 

applied  to  Alberta  watersheds.  Above- freezi ng  air 

temperatures  make  up  the  most  important  component  in 

determining  melt  due  to  longwave  radiation,  with  an 

approximation  representative  of  maximum  conditions.  The  mean 

of  the  daily  maximum  air  temperatures  is  used,  since  this  is 

believed  to  best  represent  the  air  temperature.  In  Alberta, 

this  value  for  June  in  most  cases  is  about  20%  higher  than  a 

mean  of  all  maximum  and  minimum  temperatures  and  is  about 

30%  lower  than  the  extreme  maximum  temperature.  Another 

approximation  that  was  considered  but  not  used  was  the 

extreme  maximum  temperature.  The  typical  values  recorded  at 

many  stations  in  June  in  Alberta  (sometimes  during  the 

.  o  o 

period  of  record)  are  greater  than  32  C  (90  F).  The  use  of 
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these  temperatures  was  ruled  out  since  in  reality  such 
values  are  usually  only  attained  after  prolonged  warm  spells 
lasting  beyond  the  2  days  or  more  duration  used  here;  in 
such  cases  the  snow  cover  probably  disappeared  long  before 
this. 

3.2.2  Condensation  and  Convection  Melt 

Two  other  processes  that  must  be  considered  in 
determining  snowmelt  are  the  condensation  of  vapor  and 
convection  (U.S.  Army  Corps  of  Engineers,  1956).  The  term 
for  the  combined  convective-condensation  process,  Mce,  can 
be  divided  into  two  components;  a  condensation  component 
denoted  by  Me  and  a  convective  component  denoted  by  M c 
When  water  vapor  from  the  atmosphere  condenses  on  a 
snow  surface,  the  heat  of  condensation  of  water  is  absorbed 
by  the  snowpack.  The  theory  behind  the  condensation 
processes  producing  snowmelt  is  quite  complex;  generally, 
temperature  and  wind  velocity  are  considered  to  be  the 
principal  parameters  influencing  this  process.  The  snowmelt 
produced  by  condensation  may  be  estimated  from  the  derived 
relation 

Me  =  0.054  (za  zp'1/6  (ea  -  eo)  u,  (3.3) 

where  Me  is  the  snowmelt  due  to  condensation  in  inches  per 
day,  zs  and  z^are  the  heights  (ft.)  above  the  snow  surface 
of  the  air  vapor  pressure  and  wind  speed,  respectively,  ea 
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is  the  air  vapour  pressure  in  mb,  e0  is  the  snow  surface 
vapor  pressure  in  mb  (6.11  mb  at  a  melting  snow  surface), 
and  u^  is  the  wind  speed  in  miles  per  hour. 

Melt  produced  by  convection  results  mainly  from  heat 
transferred  from  warm  air  advected  over  the  snow  surface. 
The  theory  of  turbulent  transfer  in  the  atmosphere  is 
complex,  but  experiments  have  indicated  that  simple 
approximations  can  be  used  that  consider  chiefly  the 
temperature  gradient  of  the  air  above  the  snow  surface,  the 
wind  speed,  and  the  air  density  (taken  as  a  function  of  air 
pressure).  Convection  melt  may  be  estimated  from  the 
equation  (U.S.  Army  Corps  of  Engineers,  1956) 

Mc  =  0.00629  (A)  (z  zb)'l/6(T  -  T  )  ub  (3.4) 


where  Mc  is  the  snowmelt  due  to  convection  in  inches  per 
day,  p  and  pQ  are  the  air  pressure  at  the  site  and  at  sea 
level,  respect i vely ,  zaand  zb  are  the  heights  (ft.)  above 
the  snow  surface  of  air  temperature  and  wind  speed, 
respectively,  T&  and  TQ  are  the  air  and  snow  surface 
temperature  in°F,  with  the  snow  surface  temperature  being 
normally  taken  as  32  °F  ( 0  °C )  ,  and  is  the  wind  speed  in 
miles  per  hour.  For  practical  application,  in  hydrology 
Equations  3.3  and  3.4  can  be  further  simplified. 

The  ratio  (p/pQ)  varies  from  1.0  at  sea  level  to  0.7  at 
an  elevation  of  10  000  ft.  (3  km).  For  moderate  topographic 
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changes  a  value  of  0.8  is  usually  used  and  is  appropriate 

for  application  in  Alberta.  Condensat ion-convect ion  melt 

during  rain  is  assumed  to  occur  under  saturated  air,  hence 

vapor  pressures  can  be  assumed  to  be  adequately  represented 

by  dew  point  temperatures,  and  these  in  turn  by  air 

temperatures .  Fixing  heights  of  measurement  of  Ta  and  T0  in 
o 

F  at  10  ft.  (3  m)  and  the  height  of  wind  speed  at  50  ft.  (15 
m)  as  well  as  replacing  air  temperatures  by  the  mean  daily 
maximum  temperatures  and  wind  speed  by  the  maximum  of  the 
average  daily  wind  speed  results  in  the  reduction  of 
Equations  3.3  and  3.4  to 


Mce  =  0.0084  uM  (Tm  -  32) 


(3.5) 


where  Mce  is  the  snowmelt  due  to  convection-condensation  in 

inches  per  day;  ujy;  is  the  maximum  average  daily  wind  speed 

in  miles  per  hour;  T ^  is  the  mean  of  the  daily  maximum  air 

o 

temperatures  in  F  during  the  specific  month.  Thus 
condensation  of  vapor  and  convection  are  two  important 
sources  of  heat  energy  that  affect  snowmelt. 

3.2.3  Melt  from  Rainfall 

Rain  falling  on  the  snow  surface  at  temperatures  above 
o  o 

32  F  (0  C)  transfers  heat  to  the  snow  and  thus  produces 
melt.  The  amount  of  heat  given  up  to  the  snow  by  the  rain  is 
directly  proportional  to  the  quantity  of  rain  and.  to  its 
temperature  excess  above  that  of  the  snowpacK.  For  a  melting 
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snowpack,  each  centimetre  depth  of  rainfall  releases  one 

o 

langley  (cal  per  cm2)  for  each  C  above  freezing: 


-  T  )  P 
s  r 


(3.6) 


where  Hp  is  the  heat  released  by  the  rainfall  in  langleys, 

T  r  and  Ts  are  the  temperatures  of  the  rain  and  the  snowpack, 
respectively,  in°C,  and  Pr  is  the  depth  of  rainfall  in 
centimetres.  Substituting  English  units  (since  computations 
in  the  present  work  are  in  English  units)  in  the  above 
equation  and  setting  Ts=  32°F,  Equation  3.6  becomes 

H  g  (T  -  32)  (2.54)  P 

jyj  -  - E -  =  2. - 1 - - 

p  203.2  B  203.2  (0.9?)  ,3'71 

=  0.00716  (Tr  -  32)  Pr 


where  M  is  the  daily  snowmelt  due  to  rainfall  in  inches, 
Hp  is  the  heat  released  by  rainfall  in  langleys,  Tr  is  the 
temperature  of  the  rain  in°F,  and  Pr  is  the  rainfall  in 
inches.  Substituting  the  PMP  loading  for  Pr  ,  replacing  the 
wet  bulb  temperature  by  the  mean  daily  maximum  temperature, 
and  assuming  0.97  for  the  thermal  quality  of  the  snowpack 
results  in  the  equation 


Mr  =  0.007  P jyjp  (  Tjyj  -  32) 


( 3-.  8 ) 
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where  Mr  is  the  daily  snowmelt  from  rain  in  inches,  P  is 

MP 

the  daily  probable  maximum  precipitation  in  inches,  and 

o  ^ 

is  the  mean  of  the  daily  maximum  air  temperatures  in  F 
during  the  specific  month.  This  equation  provides  a  tool  for 
calculating  the  effect  of  rainfall  on  snowmelt. 

3.2.4  Ground  Conduction 

Ground  conduction  is  another  process  involved  in  the 
heat  transfer  associated  with  snowmelt.  The  melt  produced  by 
this  process  is  the  result  of  the  upward  transfer  of  heat 
from  ground  to  snowpack  due  to  thermal  energy  stored  in  the 
ground  during  the  preceding  summer  and  early  fall.  Heat 
transfer  by  ground  conduction  can  be  expressed  by  the 
following  relationship  (U.S.  Army  Corps  of  Engineers,  1956) 


dT 

dz 


(3.9) 


where  is  the  heat  transfer  by  ground  conduction,  K  is  the 

dT 

thermal  conductivity  of  soil,  and  —  is  the  temperature 

dz 

gradient  perpendicular  to  soil  surface.  The  thermal 
conductivity  of  a  soil  varies  directly  with  its  density  and 
its  moisture  content.  Generally,  the  snowmelt  from  ground 
conduction  is  exceedingly  small.  Wilson  (1941a)  found  that 
there  was  an  insignificant  heat  transfer  from  ground  to  snow 
after  about  30  days  of  continuous  snow  cover.  The-  U.S.  Army 
Corps  of  Engineers  (1960)  estimated  the  snowmelt  from  ground 
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conduction  to  be  approximately  0.5  mm  per  day  (0.02  in.  per 

day).  The  snowmelt  from  ground  heat,  M  ,  contributes  very 

-2g  -4 

little  to  the  melt  and  about  2.1  x  10  D  mm  (8.3  x  10  D 

in.)  to  the  total  melt,  where  D  is  the  duration  of  melt  in 
hours . 

Combining  the  five  components  of  snowmelt  (i.e., 
Equations  3.1,  3.2,  3.5,  3.8,  and  the  snowmelt  from  ground 
conduction),  the  total  melt  for  PMP  loading  at  a  point 
location  can  be  expressed  as 
(a)  daily  melt  in  inches 


KL™  =  M  +  M  ,  +  M  + 

PMP  rs  rl  ce 


M  + 

P 


M 


=  (Tm  -  32)  (0.029  +  0.0084  uM  +  0.00?  PMp)  +  3. 271 


(b)  hourly  melt  in  inches 


(3.10) 


MpMp  =  10-4  (Tm  -  32)  (12.0  +3.5  UM  +  2.9  PMp)  +  0.136 

(3.11) 


For  rain-on-snow  events  under  a  PMP  loading,  the  predominant 
heat  transfer  processes  are  due  to  convection  and 
condensation,  shortwave  (solar)  radiation,  rainfall,  and 
longwave  radiation;  heat  input  by  the  other  processes  is 
relatively  minor.  The  relative  importance  of  the  heat 


' 
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processes  can  be  examined  by  considering  the  following 

values  for  a  24-hour  PMP  loading  for  a  specific  location  in 

0  ° 

Alberta  (for  example):  =  65  F  (18  C),  P  =  12.0  in.  (305 

mm),  and  u  =  15  miles  per  hour.  Using  these  values, 
Equations  3.11  becomes 

MpMp  =  0.135  +  0 . 040  +  0.173  +  0.116  +  0.0008 
=  0.465  (in. /hr.) 


For  the  above  values,  the  total  snowmelt  equals  0.46  in.  per 

hour  (11.7  mm  per  hour).  Of  this  total,  convection  and 

condensation  account  for  37.2%  of  the  total  melt;  shortwave 

(solar)  radiation  for  29.0%;  rainfall  for  24.9%;  longwave 

radiation  for  8.6%;  and  ground  conduction  for  0.2%.  For 

rain-on-snow  events  not  subject  to  PMP  loading,  the  U.S. 

Army  Corps  of  Engineers  (1956)  and  Gray  (1970)  state  that 

only  convection  and  condensation  are  important  to  the  total 

melt,  but  from  the  values  obtained  by  using  Equation  3.11  it 

becomes  obvious  that  other  processes  are  also  important. 

From  Equation  3.11,  the  total  hourly  melt  can  be 

o 

estimated  by  Knowing  the  temperature  (T  )  in  F,  probable 
maximum  precipitation  (P.flJ  in  inches,  and  the  wind  speed 

Mr 

(u  )  in  miles  per  hour,  for  a  duration  (D)  in  hours.  This  is 

lYi 

expressed  by  variables  that  are  readily  available  and  are 
not  considered  to  be  subject  to  large  spatial  and  temporal 
variability.  The  parameters  uM  ,PMp  ,  and  T ^  can  -be  called 
meteorological  parameters,  since  they  are  functions  of  the 
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meteorological  situation.  As  such,  they  can  be  expressed  in 
terms  of  a  location  in  the  river  basin  or  watershed. 

The  total  snowmelt  computed  by  Equation  3.11  is  a  point 
estimate  representing  a  probable  maximum  rate  of  melt  for  a 
given  location.  This  point  estimate  is  probably 
representative  of  the  location,  since  the  meteor logi ca 1 
parameters  comprising  the  estimate  are  generally  considered 
to  be  representative  by  the  collecting  agency.  For 
hydrological  application,  point  estimates  need  to  be 
examined  for  the  basin  in  terms  of  (a)  the  spatial 
variability,  (b)  the  maximum  melt  rate,  and  (c)  obtaining 
estimates  for  watersheds  lacking  meteorological  data. 
Although  the  basic  meteorological  parameters  of  temperature, 
wind,-  and  precipitation  are  believed  to  vary  spatially 
slightly  throughout  the  plains  of  Alberta,  it  is  worthwhile 
to  verify  the  parameters  used  in  computation  of  the  snowmelt 
to  see  if  they  behave  similarly.  In  a  basin  a  number  of 
point  estimates  may  be  computed  using  Equation  3.11;  hence 
it  is  desirable  to  find  the  maximum  value  for  a  watershed 
with  and  without  available  meteorological  data.  One  approach 
is  to  correlate  the  point  values  for  available 
meteorological  data  in  order  to  obtain  a  relationship  and 
then  to  investigate  the  relationship  for  various  data 
values.  The  next  section  discusses  the  use  of  this  approach 
to  determine  relationships  of  the  meteorological  parameters. 

3.3  Relationships  of  Meteorological  Parameters 
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In  order  to  obtain  relationships  for  the  three 
meteorological  parameters,  the  daily  maximum  air 
temperature,  the  average  daily  wind  speed,  elevation,  and 
probable  maximum  precipitation  for  various  lengths  of  time 
were  examined,  analyzed,  and  correlated  for  all  first-order 
weather  stations  and  other  available  meteorological  data  for 
major  Alberta  river  basins.  These  data  were  obtained  from 
the  numerous  Atmospheric  Environment  Service  publications, 
as  discussed  in  Section  1.3.  Regarding  the  best  way  to 
stratify  and  represent  this  information,  since  one  of  the 
desired  results  is  a  streamflow  hydrograph  of  the  PMF ,  the 
most  logical  stratification  seemed  to  be  by  river  basin. 
Stratification  by  other  means,  for  example  by  PMP  or 
temperature  distribution,  would  imply  some  form  of 
interbasin  transfer,  which  was  not  one  of  the  objectives  in 
this  study. 

In  comparing  the  values  of  the  above  parameters,  it  was 
noted  that  they  varied  by  several  orders  of  magnitude, 
suggesting  the  need  to  use  a  logarithmic  form.  Analysis 
using  a  multilinear  regression  approach  was  carried  out 
using  both  the  original  and  those  values  converted  to 
logarithmic  form.  The  relationships  obtained  by  using  the 
logarithmic  values  had  much  higher  correlation  coefficients 
than  the  relationships  with  unconverted  data,  confirming  the 
original  observations  about  the  data.  A  large  number  of 
comparisons  were  made,  and  a  summary  of  the  relationships 
obtained  in  the  analysis  of  the  three  meteorological 
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parameters  for  the  various  basins  is  presented  in  Appendix 
II. 

3.3.1  Relationship  of  the  Air  Temperature 

The  mean  daily  maximum  air  temperature  ( Tjyr  )  was  found 
to  vary  with  location  (which  in  turn  can  be  expressed  in 
terms  of  latitude  and  longitude)  and  elevation,  and  can  be 
represented  in  the  functional  form 

TM  V  E>  (3.12) 

where  %  denotes  the  mean  of  the  daily  maximum  air 
temperature  in  °F,  L  is  the  latitude  in  degrees,  L  is  the 
longitude  in  degrees,  and  E  is  the  elevation  above  sea  level 
in  feet.  Because  very  little  temperature  data  are  available 
for  the  Red  Deer  River  Basin,  the  data  from  the  North 
Saskatchewan  River  and  Red  Deer  River  basins  were  combined, 
and  a  relationship  was  obtained.  This  expression  can  be 
written  mathematically  as 

1.3^76 

L0 

TM  =  12339*6  (  2.3768  0.27242^  (3.13) 

La  e 


where  T^  is  the  mean  of  the  daily  maximum  air  temperature  in 
F,  La  is  the  latitude  in  degrees,  L0  is  the  longitude  in 
degrees,  and  E  is  the  elevation  above  sea  level  in  feet. 
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This  relationship  has  a  high  multicorrelation  coefficient  of 
0.995,  which  seems  to  suggest  that  this  quantity  is  well 
defined  by  the  variables  considered.  Elevation  seems  to  be 
important  in  the  computation  of  the  air  temperature,  for 
without  this  parameter  the  relationship  had  a  correlation 
coefficient  about  30%  lower  (see  Appendix  II).  The 
computations  of  Equation  3.13  can  be  depicted  graphically; 
these  are  shown  in  Figure  3.4.  The  values  give  the  air 


Figure  3.4.  Temperature  Variation  (  °F )  for  the  Red  Deer 

River  Basin. 
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temperature  loading  for  the  Red  Deer  River  Basin  and  show 
that  the  greatest  decrease  in  air  temperature  occurs  in  the 
westward  direction  and  is  dependent  upon  elevation.  The  main 
advantage  of  such  a  relationship  is  that  estimates  of 
temperature  can  be  made  at  locations  where  no  meteorological 
data  is  available.  No  relationship  is  presented  for  the 
mountainous  areas,  because  very  little  meteorological  data 
are  available  for  these  areas. 

3.3.2  Relationship  of  the  Wind  Speed 

The  second  parameter  for  which  a  spatial  relationship 
was  developed  was  the  wind  speed.  The  magnitude  of  the 
largest  average  wind  speed  integrated  over  an  hour  was  used 
for  computing  snowmelt.  The  maximum  average  wind  speed, 
similar  to  the  maximum  air  temperature  was  found  to 
correlate  well  with  location  and  elevation  (Appendix  II), 
and  this  relationship  can  be  generally  expressed  as 

UM  ^  f(LA’  V  E)  (3.14) 

with  Ujyr  being  the  value  of  the  largest  average  wind  speed  in 
miles  per  hour  and  the  other  parameters  being  the  same  as  in 
the  previous  case.  The  functional  relationship  suggested  for 
the  Red  Deer  River  basin  has  a  correlation  coefficient  of 
0.800  and  can  be  expressed  as 
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36.01 

o  ho  /  L0  X  (3.15) 

UM  '  8,42  ^  T 01^4.461 

LA  E 

where  uM  the  average  wind  speed  in  miles  per  hour,  L  is  the 

A 

latitude  in  degrees,  L0  is  the  longitude  in  degrees,  and  E 
is  the  elevation  above  sea  level  in  feet.  Examination  of  the 
existing  wind  speed  data  showed  that  this  quantity  varied 
with  location  in  May  and  June.  The  approach  used  in  the 
generalized  snowmelt  equations  is  to  consider  wind  speed  as 
a  constant  (i.e.,  assign  a  value  of  5  miles  per  hour).  From 
the  examination  and  analysis  of  meteorological  data,  the 
values  that  emerge  for  average  wind  speeds  are  up  to  three 
times  larger  than  those  used  in  generalized  snowmelt 
equations.  For  the  Red  Deer  River  Basin,  the  spatial 
distribution  of  values  of  wind  speed  are  graphically  shown 
in  Figure  3.5;  these  appear  to  be  dependent  on  the 
elevation.  For  this  river  basin  the  average  wind  speed  seems 
to  vary  from  about  12  to  18  miles  per  hour,  with  most  of  the 
lower  values  computed  in  the  eastern  portions  of  the  river 
basin.  The  value  of  18  miles  per  hour  is  over  three  times 
the  5  miles  per  hour  value  commonly  used  in  the  generalized 
snowmelt  equations.  Maximum  observed  hourly  speeds  are  also 
available;  however,  these  were  not  used  because  the  values 
are  very  large  and  occur  in  a  much  shorter  time  period  than 
was  considered  in  this  work.  A  summary  of  wind  speed 
relationships  for  other  river  basins  is  provided  -in  Appendix 
II.  No  relationship  was  obtained  for  the  Peace  River  Basin 
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Figure  3.5.  Variation  of  the  Wind  Speed  (miles  per  hour)  for 

the  Red  Deer  River  Basin. 

because  of  the  lack  of  meteorological  data;  for  this  area, 
values  obtained  for  the  Athabasca  River  Basin  are 
recommended . 

3.3.3  Relationship  of  the  Precipitation 

The  third  and  last  parameter  examined  was  the  point 
values  of  PMP  as  obtained  by  the  statistical  approach.  The 
factors  that  can  be  identified  in  determining  the  PMP  are 
the  location  (which  can  be  expressed  in  terms  of  -latitude 
and  longitude),  duration,  storm  area,  and  the  elevation.  In 
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the  DAD  analyses  (Chapter  2  and  Verschuren  and  Wojtiw,  1980) 
the  PMP  was  observed  to  vary  with  the  duration  of  the 
rainfall  and  with  the  storm  area.  Although  the  decrease  in 
the  PMP  value  was  not  expected  to  be  large  for  the  storm 
area,  it  was  considered  in  the  analysis.  The  elevation  and 
location  were  considered  necessary  in  estimating  PMP.  A 
summary  of  a  number  of  combinations  of  these  parameters  is 
given  in  Appendix  II.  The  combination  of  all  these 
parameters  (a  sample  size  of  75  values  with  5  degrees  of 
freedom  was  used  for  the  regression)  gives  the  largest 
correlation  coefficient  in  a  multilinear  regression  analysis 
and  can  be  expressed  by 


^MP  ^  ^  *  Lo’ 


(3.16) 


where  P^p  is  the  PMP  in  inches,  l_A  is  the  latitude  in 
degrees,  LQ  is  the  longitude  in  degrees,  D  is  the  duration 
in  hours,  E  is  the  elevation  above  sea  level  in  feet,  and  A 
is  the  watershed  area  in  sq .  mi..  For  the  Red  Deer  River 
Basin,  the  mathematical  expression  has  the  following 
equation 

4.246  x  10-3  l33-455  d0.419 
PMP  =  l31.584  e3-257  a0.0?1 


(3.17) 
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and  produces  a  correlation  coefficient  of  0.940.  The 
watershed  area  does  not  contribute  significantly  to  the  PMP 
estimate  in  this  basin.  This  is  also  true  for  the  other 
basins  (Appendix  II).  The  other  three  variables  (location, 
elevation,  and  duration)  do  contribute  substantially,  and 
the  elimination  of  any  of  these  variables  would 
significantly  decrease  the  correlation  coefficient.  For  a 
number  of  locations  in  the  Red  Deer  River  Basin  this 
equation  was  computed  with  the  contoured  results  depicted  in 
Figure  3.6,  which  shows  a  gradual  decrease  in  the  24-hour 
PMP  from  over  305  mm  (12  in.)  in  the  headwater  of  the  basin 
to  about  230  mm  (9  in. )  in  the  eastern  portions  of  the 
basin.  This  relationship  provides  spatial  estimates  of  PMP 
for  a  watershed  within  a  basin  when  meteorological  data  are 
lacking,  since  the  four  remaining  unknown  parameters  (L.  , L 

r\  O 

,D  ,  and  E)  can  be  determined  from  topographical  maps.  This 
type  of  relationship  provides  a  graphical  or  numerical 
method  of  determining  the  maximum  PMP  in  the  basin.  Without 
such  a  relationship  it  is  not  certain  if  the  PMP  estimate 
obtained  by  either  the  statistical  or  the  physical  method 
has  been  computed  for  maximum  value. 

3.4  Maximum  Water  Loading 

The  maximum  water  loading  can  now  be  determined,  since 
this  is  a  combination  of  the  PMP  and  maximum  snowmelt  (that 
is,  the  summation  of  Equations  3.11  and  3.17).  These  two 
quantities  can  be  defined  in  terms  of  the  air  temperature, 
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Figure  3.6.  24-hour  Spatial  PMP  Variation  (inches)  in  the 

Red  Deer  River  Basin. 

wind  speed,  and  probable  maximum  precipitation  for  a  given 
duration,  which  in  turn  can  be  expressed  by  functional 
relationships  for  more  easily  obtainable  parameters,  i.e., 
location,  elevation,  and  drainage  area  of  the  watershed. 
Elevation  and  area  of  basin  can  be  predetermined  for  a  given 
location  from  topographic  and  drainage  maps.  Minimum 
elevation  is  used  since  it  gives  the  maximum  result.  Thus 
the  maximum  water  loading  lasting  for  D  hours  for  a  given 
watershed  area  reduces  to  a  simple  functional  relationship 
in  terms  of  the  location  (i.e.,  latitude  and  longitude), 
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elevation,  and  drainage  area: 

WM  <<  f(LA,  L0.  E,  A)  (3-  18) 

where  WM  is  the  maximum  water  loading  in  inches,  l_A  is  the 
latitude  in  degrees,  Lq  is  the  longitude  in  degrees,  E  is 
the  elevation  above  sea  level  in  feet,  and  A  the  storm  area 
in  sq .  mi..  The  numerical  equation  for  the  Red  Deer  River 
basin  can  be  obtained  easily  using  the  above  equations 
(Equations  3.11,  3.13,  3.15,  and  3.17)  but  is  not  presented 
here  because  of  complexity. 

The  use  of  a  functional  relationship  for  water  loading 
is  advantageous  in  that  spatial  estimates  can  be  obtained 
for  a  watershed  where  meteorological  data  are  lacking  and 
can  provide  an  estimate  of  the  maximum  value. 

The  concept  of  functional  relationship  discussed  in 
this  chapter  has  been  presented  in  a  general  form  to  suggest 
that  this  concept  can  be  applied  to  other  basins.  The 
relationships  depend  upon  the  quality  and  quantity  of  data 
available  and  on  the  homogeneity  of  meteorological 
conditions  and  topographical  features  within  the  basin. 
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CHAPTER  4 


THE  TIME-AREA  PROBABLE  MAXIMUM  FLOOD  MODEL 

4 . 1  Introduction 

As  was  noted  in  the  literature  review,  there  have  been 
few  estimates  of  PMP  and  PMF  for  Alberta  river  basins.  One 
of  the  reasons  for  this  is  the  lack  of  a  simplified  method 
or  model  that  can  be  modified  for  PMP  conditions  and  from 
which  a  PMF  hydrograph  can  be  estimated.  It  was  also 
believed  that  the  computations  of  the  slope  of  the  watershed 
would  significantly  influence  the  PMF  estimate. 

To  address  these  concerns  (and  as  an  objective  of  this 
thesis),  it  was  necessary  to  develop  a  numerical  model  for 
estimating  PMF.  The  resulting  model  had  to  incorporate  an 
estimate  of  the  shape  of  the  watershed.  A  model  with  this 
requirement  is  presented  in  this  chapter;  it  has  been  termed 
the  Time-Area  Probable  Maximum  Flood  Model  ( TAPMF  Model). 
Because  it  was  also  desirable  to  examine  two  commonly  used 
models  for  PMP  loading  and  to  compare  the  peak  discharges, 
the  SSARR  and  HYMO  models  are  presented  in  Chapter  5,  and 
the  results  from  all  three  models  are  compared  in  Chapter  6. 

4.2.0  The  TAPMF  Model 

The  TAPMF  model  is  graphically  depicted  in  Figure  4.1. 
Functional  relationships  are  used  to  obtain  estimates  of  PMP 
and  the  maximum  snowmelt  for  locations  within  the  watershed 
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Figure  4.1.  The  TAPMF  Model. 
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(as  presented  in  Chapter  3),  and  these  estimates  are 
combined  to  give  the  total  water  loading  for  the  watershed. 
For  simplicity,  the  maximum  value  in  the  watershed  is  used. 
Using  a  design  storm  hyetograph  and  a  time-area  histogram, 
an  inflow  hydrograph  is  obtained  from  this  loading,  which 
yields  the  total  water  available  for  runoff.  This  hydrograph 
is  routed  through  the  watershed  as  surface  flow  by  using  a 
modified  storage-outflow  relationship  and  the  continuity 
equation  to  produce  the  PMF  hydrograph.  The  runoff  is 
expressed  as  direct  surface  flow,  since  saturated  soil 
conditions  are  assumed  for  the  PMP  loading.  The  resulting 
flow  then  generates  an  outflow  hydrograph  termed  as  the  PMF. 
For  computations,  the  basin  is  divided  into  a  number  of 
watersheds,  so  that  the  drainage  is  toward  the  main  river 
system.  To  i 1 lustrate  the  development  of  the  model , 
watersheds  from  the  Red  Deer  River  Basin  are  used,  as 
defined  by  Water  Survey  of  Canada  (Figure  4.2). 

Although  the  TAPMF  Model  presented  in  this  thesis 
follows  the  overall  design  of  a  number  of  hydrological 
models,  not  all  the  components  can  be  found  in  a  single 
model,  nor  are  they  estimated  or  calculated  using  the  same 
procedures.  In  designing  the  TAPMF  Model,  the  author  has 
used  parts  of  models  that  seemed  appropriate  for  Alberta 
river  basins,  that  is,  components  that  seemed  physically 
suitable.  The  main  features  of  the  design  are  simplicity  of 
data  requi rements ,  usage,  and  potential  for  modification. 
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Figure  4.2.  Map  of  the  Watersheds  of  the  Red  Deer  River 

Basin  (After  Muzik,  1975). 


4.2.1  Determination  of  the  Inflow  Hydrograph 

Watershed  simulation  generally  utilizes  some  form  of 
synthesizing  technique  to  predict  the  temporal  and  spatial 
variations  of  a  flood  wave  as  it  traverses  the  watershed. 
Hydrologic  techniques  are  usually  applied  to  the  synthesis 
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of  runoff  hydrographs  from  gauged  and  ungauged  watersheds. 

In  these  techniques  an  equation  of  continuity  is  used  with 
either  an  analytic  or  an  assumed  relationship  between 
storage  and  discharge  within  the  system.  A  number  of 
techniques  exist  for  hydrologic  synthesizing  of  the 
watershed . 

The  TAPMF  Model  uses  the  time-area  histogram  technique, 
also  Known  as  Clark's  method  (Viessman  et  al.f  1977).  This 
method  was  selected  because  by  construction  of  isochrones  it 
directly  incorporates  the  basin  shape  in  its  application. 

The  shape  of  the  basin  is  a  parameter  that  is  believed  by 
the  author  to  contribute  significantly  to  the  peak  discharge 
and  the  PMF  hydrograph.  Although  basin  shape  is  considered 
indirectly  in  the  SSARR  and  HYMO  models,  in  both  of  these 
models  it  is  not  easy  to  isolate  the  contribution  of  this 
parameter  to  the  peak  discharge. 

In  the  time-area  histogram  technique,  an  inflow 
hydrograph  is  needed  (Figure  4.1)  The  construction  of  this 
hydrograph  requires  three  elements:  (a)  creation  of  a  design 
storm  hyetograph,  (b)  separation  of  the  watershed  into 
isochrones,  and  (c)  construction  of  a  time-area  histogram. 

4.2.2  Design  Storm  Hyetograph 

The  design  storm  hyetograph  (Figure  4.1)  can  be 
considered  the  water  loading  for  a  given  time  period.  The 
total  water  loading  for  a  watershed  is  a  combination  of  the 
PMP  loading  and  the  maximum  snowmelt.  (The  latter  was 
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discussed  in  Chapter  3  and  can  be  estimated  for  a  watershed 
by  functional  relationships.)  As  part  of  the  PMP  loading, 
the  precipitation  distribution  (the  variation  of  the  depth 
of  precipitation  with  time)  affects  the  shape,  peak,  and  lag 
time  of  the  resultant  flow.  The  time  distribution  of 
precipitation  with  large  recorded  depths  seems  to  be  ignored 
in  the  literature,  however,  this  distribution  plays  an 
important  part  in  PMF  estimates.  In  the  few  PMF  studies 
(e.g.,  Muzik,  1975),  the  usual  approach  was  to  distribute 
the  total  precipitation  equally  in  a  number  of  daily 
periods,  which  resulted  in  a  uniform  and  low- i ntensi ty 
loading  on  the  watershed. 

As  an  alternative,  the  author  grouped  Alberta 
rainstorms  that  produced  depths  of  100  mm  (4  in.)  or  more  in 
6-hour  intervals  and  stratified  these  for  different  storm 
durations.  The  maximum  mass  distribution  for  durations  of 
24-72  hours  at  the  maximum  recorded  depth  is  shown  in  Figure 
4.3.  This  can  provide  an  indication  of  when  the  intense 
rainfall  occurred.  Durations  greater  than  72  hours  were  also 
examined;  however,  these  are  not  included  here  because  the 
author  was  unable  to  define  a  representat i ve  curve  for 
rainstorms  with  these  durations,  because  of  the  limited 
available  data  in  Alberta. 

For  rainstorm  durations  of  48  hours  or  more,  small 
percentages  of  rain  occur  in  the  first  18  hours,  and  the 
precipitation  is  produced  gradually  with  time,  as-  is 
indicated  in  Figure  4.3.  These  storms  lack  the  high- 
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Figure  4.3.  Maximum  Mass  Distribution  at  the  Recorded 
Maximum  Depth  for  Various  Durations. 


intensity  characteristics  evident  in  some  thunderstorms  and 
rainstorms  with  durations  less  than  48  hours.  Rainstorms 
with  durations  of  48  hours  or  more  are  termed  the  cold  lows 
(Chapter  1)  and  are  capable  of  producing  precipitation  over 
150  mm  (6  in.)  at  the  maximum  depth.  For  shor t -durat ion 
rainstorms  (36  hours  and  less)  most  of  the  precipitation 
occurs  at  the  beginning  of  the  storm  and  is  produced  in  a 
very  short  period  of  time.  This  suggests  that  for-  rainstorms 
with  durations  of  36  hours  and  less,  shor t -durat ion 


intensity  rainfall  data  could  be  used  for  obtaining  the 
precipitation  distribution.  The  use  of  short-time  intensity 
data  is  further  advocated  because  the  maximum  discharge 
would  result  from  the  greatest  amount  of  precipitation 
occurring  in  the  shortest  time.  The  available  intensity  data 
were  examined  for  stations  in  the  vicinity  of  the  Red  Deer 
River  Basin.  The  largest  recorded  maximum  depths  were  29  mm 
(1.14  in.)  for  a  1-hour  duration,  49  mm  (1.94  in.)  for  a  6- 
hour  duration,  and  73  mm  (2.89  in.)  for  a  12-hour  duration. 
Although  intensities  for  durations  in  minutes  are  available, 
hourly  values  are  used  here  to  coincide  with  the  time  period 
used  in  the  models  and  for  practical  considerations 
(computing  time  and  the  results). 

Three  parameters  need  to  be  considered  in  the  design  of 
the  precipitation  distribution:  the  total  PMP  loading, 
duration,  and  intensity.  Expression  of  these  parameters 
needs  to  be  such  that  the  product  of  the  duration  and 
intensity  is  equal  to  the  total  PMP  loading.  Various 
combinations  of  these  parameters  are  possible,  the  simplest 
being  division  of  the  total  PMP  loading  by  the  duration, 
which  produces  a  low  average  intensity  compared  to  the  1- 
hour  i ntens i ty-dur at i on  data.  Such  an  average  intensity 
would  in  most  cases  coincide  with  a  frequent  event  (one  in  2 
years  or  one  in  5  years)  and  not  an  extremely  rare  event  (as 
the  PMP  is  considered  to  be).  A  more  logical  approach,  and 
one  that  is  advocated  here,  is  to  use  the  1 argest-  recorded 
1-hour  intensity  (since  this  is  the  shortest  time  used)  and 
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distribute  this  value  over  the  shortest  duration  period 
possible  for  total  PMP  loading.  For  example,  for  a  24-hour 
PMP  loading  of  343  mm  (13.5  in.)  with  an  intensity  of  38  mm 
(1.5  in.)  per  hour,  the  storm-design  distribution  would 
cons ist  of  rainfall  i ntens i ty  of  38  mm  (1.5  in.)  per  hour 
for  the  first  9  hours  and  a  zero  intensity  for  the  remaining 
15  hours.  For  validity,  this  approach  can  be  compared  to 
analyzed  short-duration  intensity  data,  depicted  in  Figure 
4.4  for  the  Edmonton  Municipal  Airport  (with  45  years  of 
record).  Unfortunately,  few  intensity-frequency-duration 
graphs  are  available  for  Alberta;  hence  it  was  necessary  to 
select  the  station  nearest  to  the  watershed  that  had  a  long 
record  of  data.  In  this  case  the  Edmonton  Municipal  Airport 
was  chosen  (Figure  4.4).  In  terms  of  the  shortest  time 
period  considered  (1  hour),  the  38  mm  (1.5  in.)  per  hour 
intensity  suggests  about  a  1  in  100  year  event,  while  a  6- 
or  12-hour  duration  comparison  suggests  that  the  38  mm  (1.5 
in.)  per  hour  intensity  may  occur  once  in  several  thousand 
or  even  a  million  years.  The  use  of  the  largest  recorded  1- 
hour  intensity  seems  to  provide  a  practical  aspect  to  the 
storm-design  hyetograph. 

4.2.3  Separation  of  Watershed  into  Isochrones 

The  second  requirement  for  devising  an  inflow 
hydrograph  for  the  TAPMF  Model  is  the  separation  of  the 
watershed  into  isochrones.  The  watershed  is  f i rst- divided 
into  a  number  of  time  zones  separated  by  isochrones,  or 
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Figure  4.4.  Intensi ty- Frequency-Durat ion  Curves  for  Edmonton 
Municipal  Airport  (Data  from  Atmospheric  Environment 

Servi ce ) . 


lines  of  equal  travel  time  from  the  watershed  outlet,  as 
shown  in  Figure  4.1.  The  maximum  isochrone  is  believed  to  be 
an  approximation  of  the  watershed  shape  that  affects  the 
peak  flow.  Numerous  authors  have  attempted  to  obtain  a  bulk 
parameter  for  the  watershed  shape,  and  although  some  have 
had  moderate  success  with  particular  watersheds,  no 
universal  parameter  has  been  found.  The  importance  of  the 
shape  of  the  watershed  on  the  outlet  hydrograph  can  be  seen 
in  Figure  4.5,  which  illustrates  two  different  watershed 
shapes.  This  figure  shows  that  the  larger  the  value  of  the 
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Figure  4.5.  Effect  of  Shape  on  the  Outflow  Hydrograph. 

maximum  isochrone,  the  greater  the  peak  discharge.  The  time- 
area  approach  is  believed  to  provide  a  reasonable 
representat ion  of  the  shape  of  the  watershed.  The  areas 
between  the  isochrones  can  be  determined  and  plotted  against 
the  travel  time  (i.e.,  time-area  histogram). 

To  divide  the  watershed  into  equal  time  zones,  an 
estimate  of  the  speed  of  water  flow  in  the  streams  of  the 
watershed  is  needed.  The  velocity  of  water  in  a  stream  or 
river  has  been  approximated  by  a  logarithmic  relationship 
(Water  Survey  of  Canada,  1977),  that  is,  a  straight  line  is 
obtained  when  logarithmic  values. of  velocity  are  plotted 
against  the  logarithmic  values  of  discharge.  Figure  4.6 
illustrates  this  dependency  (from  graphs  produced  by  Inland 
Water  Survey)  for  a  number  of  streams  in  the  Red  Deer  River 
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Basin  watersheds.  From  this  figure,  a  value  of  stream 


Figure  4.6.  D i scharge- Ve loci ty  Distributions  for  a  Number  of 
Streams  in  the  Red  Deer  River  System  (Data  from  Water  Survey 

of  Canada ,1977). 


velocity  can  be  estimated,  and  this  in  turn  can  be  used  to 
calculate  the  distance  water  travels  for  a  specified  time 
interval.  From  this,  isochrones  for  the  watershed  can  be 
drawn.  The  water  travel  time  on  the  river  can  be  similarly 
obtained  by  considering  di scharge- ve loci ty  distributions  for 
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the  river.  This  is  illustrated  for  the  Red  Deer  River  system 
in  Figure  4.7,  In  most  models,  the  velocity  is  often 


Figure  4.7.  Di scharge- Veloci ty  Distribution  at  a  Number  of 
Locations  on  the  Red  Deer  River  (Data  from  Water  Survey  of 

Canada ,  1977 ) . 


approximated  by  a  velocity-discharge  relationship  with  a 
fixed  upper  limit. 

4.2.4  Time-area  Histogram 

In  creating  an  inflow  hydrograph  for  the  TAPMF  Model, 
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the  third  requirement  is  to  obtain  a  time-area  histogram. 
From  the  watershed  isochrones,  the  area  for  each  time  zone 
can  be  planimetered ;  when  plotted  against  time,  a  time-area 
histogram  results  (Figure  4.1).  The  reservoir  inflow 
hydrograph  ordinates  (I)  for  any  selected  design  hyetograph 
can  now  be  determined.  Each  block  of  rain  in  Figure  4.1  is 
applied  to  the  entire  watershed;  the  runoff  from  each 
subarea  reaches  the  outflow  at  lagged  intervals  defined  by 
the  time-area  histogram.  The  simultaneous  arrival  of  the 
runoff  from  areas  A^,  k  ,  etc.,  for  various  design  storm 
intensities  ,etc. ,  can  be  determined  by  properly 

lagging  and  adding  contributions,  as  expressed  by  the 
equation 


+  + 


+  R1Ai 


.(4.1) 


where  I.  are  the  inflow  hydrograph  values  in  cfs,  r  are  the 
i  i 

design  storm  intensity  values  in  inches  per  hour,  and  a  are 

i 

the  time-area  histogram  values  in  acres.  By  applying 
Equation  4.1  it  is  possible  to  compute  the  inflow  hydrograph 
shown  in  Figure  4.1.  The  time-area  histogram  method 
(although  a  little  tedious  to  apply)  is  easier  to  interpret. 

4.2.5  Flow  Separation 

Anticipating  no  evapotranspi rat  ion  losses,  the  total 
generated  runoff  can  be  partitioned  into  the  components  of 
base,  surface,  and  subsurface  flow  runoff.  For  PMF  estimates 
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the  maximum  surface  runoff  would  occur  when  base  and 
subsurface  flows  are  assigned  zero  flows  (equivalent  to 
saturated  soil  conditions).  Hence,  the  total  generated 
runoff  is  expressed  as  surface  runoff. 

4.2.6  Routing  Method 

The  basic  routing  method  used  in  the  TAPMF  model  is 
similar  to  that  described  by  Wilson  (1941),  expanded  upon  by 
Rockwood  (1958)  and  used  in  the  SSARR  Model  (U.S.  Army  Corps 
of  Engineers,  1972).  In  this  method,  routing  through 
watershed,  river  system,  and  reservoir  components  of  the 
model,  relies  on  the  use  of  the  law  of  continuity  in  the 
storage  relationship 

(  n  2  ■■  "+1)  At  -  (-S-g . ati)  At  =  Sn+1  -  Sn  (4.2) 


where  I n  is  the  inflow  at  the  beginning  of  the  time  period, 

In+lis  the  inflow  at  the  end  of  the  period,  0  is  the 

n 

outflow  at  the  beginning  of  the  period,  0  is  the  outflow 

n+ 1 

at  the  end  of  the  period,  Sn  is  the  storage  at  the  beginning 
of  the  period,  S^is  the  storage  at  the  end  of  the  period, 
and  At  is  the  duration  in  hours. 

The  inflow  at  time  t  ( I  )  can  be  expressed  by 


I 


t 


dS 

dt 


(4.3) 
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In  natural  lakes,  where  storage  is  a  function  of  outflow  at 
any  given  elevation,  K  represents  the  propor t iona 1 i ty  factor 
between  storage  and  outflow: 


Differentiating  with  respect  to  time,  the  following 
relationship  is  obtained: 


dS  _  t/  t  do  N 

te  -  K 


(4.5) 


Substituting  this  expression  into  Equation  4.3  and  rewriting 
the  terms  results  in 


dO 

dt 


(4.6) 


which  represents  the  basic  form  of  the  storage  equation. 
Watershed  routing  can  be  evaluated  through  successive 
increments  of  lake-type  storage.  A  watershed  can  be 
considered  a  series  of  small  lakes  that  represent  the 
natural  delay  of  runoff  from  upstream  to  downstream  points. 
During  the  routing  computation,  the  computer  operates  on 
each  specified  increment  of  storage  consecutively,  one 
period  at  a  time.  The  equation  used  by  the  computer  (in  the 
TAPMF  Model  and  also  in  SSARR  Model)  can  be  derived  by 
considering  Equation  4.2,  the  storage  routing  equation. 
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Setting  Im  ^ In  +  In+1^  /  2 

and  substituting  this  in  Equation  4.2  gives 

(y  At  .  (  °n  *  °n+1)At  =  Sn+1  -  Sn 


Dividing  by  At  ,  and  subtracting  0n  from  both  sides,  the 
equation  becomes 


^  '  °n 


f-3n+l  ~  °nN  A  /  °n+l  °ri\ 
v  At  '  '  2 - ' 


Multiplying  both  sides  by  -7-q - — ~ — y  and  rearranging  the 

v  n+ 1 “  urr 

terms  results 


I  -  0 
m  n 


°n+l  -  °ru  ,  Sn+1  ~  Sn 

°n+ 1  -  °n 


At 


•)  ( 


+  At  /  2) 


Wi  th 


K  = 


(Sn+l  -  Sn>  /  <°„+l  -  °„> 


the  above  equation  can  be  re-written  as 


°ml  -  °n 

At 


-  °n 

K  +  At/2 


Solving  for  0^+^the  equation  in  the  TAPMF  model  becomes 


' 
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0 


m 


0 


n 


n+1 


K  +  At/2 


)  At 


0 


n 


(4.7) 


where  0  is  the  outflow  at  the  end  of  the  period,  I  is  the 
n+l  m 

mean  inflow,  is  the  outflow  at  is  the  beginning  of  the 
period,  At  is  the  time  duration  of  computational  period,  and 
K  is  the  storage  coefficient.  For  given  values  of  inflow  and 
outflow  at  the  beginning  of  the  period,  the  outflow  at  the 
end  of  the  period  (0^+^)  is  found  by  Equation  4.7.  In  the 
computer  program,  the  value  of  At  / ( K  +  At/2)  is  evaluated 
for  a  specifi eu  condition,  and  mu  1  tipi i ed  by  the  difference 
between  Im  and  0*  to  obtain  the  change  in  outflows.  The 
outflow  ( O2 )  computed  from  the  first  incremental  routing  is 
then  saved  and  converted  to  initial  outflow  (01)  for  the 
next  period  routing. 

The  coefficient  K  must  be  estimated  and  can  be  a  few 
hours  to  a  month  or  more.  For  watersheds  where  hydrograph 
data  are  not  avai Table,  K  can  be  estimated  from  an  empirical 
equation  by  equating  this  quantity  with  the  basin  lag  time. 

A  number  of  equations  are  available.  The  equation  that  is 
used  in  the  TAPMF  Model  is  the  one  developed  by  the  Soil 
Conservation  Service  (1969): 


1  1900  Y°  *  ^  (4.8) 

where  t1  is  the  lag  time  in  hours,  1  is  the  length  to  the 
divide  in  feet,  Y  is  the  average  watershed  slope  in  percent, 


- 
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and  S  is  the  potential  maximum  retention  in  inches.  The 
potential  maximum  retention  is  equivalent  to  ( 1 000/CN ) -  1 0 , 
where  CN  is  a  curve  number.  For  saturated  soil  conditions, 
CN  =  100,  and  hence  S  becomes  equal  to  zero.  Substituting  S 
=  0.0  reduces  Equation  4.8  to: 


10.8 

1900  Y° ' ^ 


(4.9) 


For  a  watershed  in  the  Red  Deer  River  Basin  with  an  area  of 
2470  Km2  (954  sq .  mi.)  an  average  slope  of  5.7  m  per  Km 
(30.1  ft.  per  mi)  or  0.5701  expressed  as  a  percentage,  and  1 
equal  to  85  Km  (53  mi . ) ,  a  lag  time  of  16  hours  is 
calculated  from  Equation  4.9.  This  lag  time  can  be  used  as 
an  estimate  of  the  storage  coefficient  in  Equation  4.7. 


4.2.7  Operations  and  Data  Requirements 

The  above  components  are  then  organized  and  coordinated 
to  produce  a  river  basin  model.  The  system  configuration 
(similar  to  Figure  5.7  for  SSARR  Model)  essentially 
describes  to  the  computer  the  physical  layout  and 
relationships  of  all  components  of  the  system.  The 
configuration  is  in  upstream  to  downstream  order  of  all 
watersheds,  laKes,  reservoirs,  channel  reaches,  and 
confluence  points  for  a  particular  basin.  The  hierarchical 
procedure  involves,  first,  watershed  routing,  then 
consecutive  channel  routing,  and  finally  combining  until  all 
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operations  are  complete.  The  computer  model  then  integrates 
the  two  hierarchies  and  performs  the  model  simulation  of  the 
river  basin  system. 

The  input  data  needed  for  operation  of  the  TAPMF  Model 
include:  (a)  nonvariable  character i s t i cs  data  which  describe 
physical  features  such  as  drainage  area,  (b)  time  variable 
data  which  include  physical  data  expressed  as  time  series; 
for  example,  probable  maximum  precipitation,  snowmelt 
loading,  storm  hyetograph,  and  time-area  histogram,  and  (c) 
miscellaneous  job  control  and  time  control  data  which 
specify  such  items  as  total  computation  period,  routing 
intervals,  and  special  computer  instructions  to  control 
input-output  alternatives. 

Results  obtained  using  the  TAPMF  Model  are  compared 
with  those  from  the  SSARR  and  HYMO  models  in  Chapter  6. 
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CHAPTER  5 


THE  SSARR  AND  HYMO  MODELS 

5.1.0  Introduction 

An  objective  in  this  study  was  to  compare  the  predicted 
PMF  peak  discharge  from  the  TAPMF  Model  (presented  in 
Chapter  4)  with  that  obtained  from  two  commonly  used  models 
in  Alberta  (SSARR  and  HYMO).  Before  the  comparison  is 
presented  in  Chapter  6,  the  PMP  loading  for  the  SSARR  Model 
is  discussed  in  Section  5.2  and  for  the  HYMO  Model  is 
discussed  in  Section  5.3. 

The  SSARR  (Streamflow  Synthesis  and  Reservoir 
Regulation)  Model  was  one  of  the  earliest  continuous 
streamflow  simulation  models  using  a  lumped  parameter 
representation;  its  primary  strength  lies  in  its  verified 
accuracy.  Successful  tests  have  been  conducted  in  several 
large  drainage  basins,  including  the  Columbia  River  Basin 
( NW  United  States)  and  the  Mekong  River  Basin  (SE  Asia) 

(U.S.  Department  of  Commerce,  1970).  This  model  has  also 
been  used  by  Alberta  Environment  in  hydrological  studies  on 
several  river  basins  in  Alberta  (Alberta  Environment,  1980). 
For  these  reasons,  this  model  was  chosen  for  comparison  in 
this  work.  The  success  of  any  model  is  usually  dependent 
upon  the  applicability  of  the  relationships  to  a  watershed 
or  basin;  these  relationships  are  examined  in  some  detail  in 
Section  5.2.  For  PMF  estimates,  the  sensitivity  of  a  model 


124 


i 


125 


response  to  parameter  changes  is  important,  since  some 
parameters  may  affect  the  shape  of  the  resulting  hydrograph 
and  the  peak  estimate  of  the  PMF .  The  parameterization  and 
sensitivity  of  the  SSARR  and  HYMO  models  for  PMP  loading  are 
emphasized  here,  since  this  is  important  to  an  understanding 
and  interpretation  of  the  estimates  of  PMF. 

5.2.0  SSARR  Watershed  Model 

The  SSARR  Model  is  a  widely  used,  continuous  streamflow 
simulation  model  designed  for  large  basins.  It  was  developed 
by  the  U.S.  Army  Corps  of  Engineers  (1972)  for  streamflow 
and  flood  forecasting  and  for  reservoir  design  and  operation 
studies.  It  is  a  mathematical  hydrologic  model  of  a  river 
basin  system  throughout  which  streamflow  can  be  synthesized 
by  evaluating  snowmelt  and  rainfall. 

Applications  of  the  model  begin  with  a  subdivision  of 
the  drainage  basin  into  homogeneous  hydrologic  units  of  a 
size  and  character  consistent  with  subdivides,  channel 
confluences,  reservoir  sites,  diversion  points,  soil  types, 
and  other  distinguishing  features.  The  streamflows  can  be 
computed  for  all  significant  points  throughout  the  river 
system.  An  SSARR  watershed  model  flow  chart  is  shown  in 
Figure  5.1. 

Rainfall  data  can  be  supplied  at  any  number  of  stations 
in  the  watershed.  The  part  that  will  run  off  is  divided  into 
the  base  flow,  subsurface  or  interflow,  and  surface  runoff. 
The  separation  is  based  on  indexes  and  on  the  intensity  of 
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Figure  5.1.  SSARR  Watershed  Model  Flow  Chart  (After  Alberta 

Environment,  1977b). 
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the  direct  runoff.  Each  component  is  simply  delayed 
according  to  different  processes,  and  all  are  then  combined 
to  produce  the  final  watershed  outflow  hydrograph.  This 
runoff  can  then  be  routed  through  stream  channels  and 
combined  with  other  watershed  hydrographs,  all  of  which 
become  part  of  the  resultant  hydrograph. 

Routing  through  channels  requires  an  assumption  of 
short  stream  reaches.  Streamflows  are  synthesized  on  the 
basis  of  rainfall  and  snowmelt  runoff.  Snowmelt  can  be 
determined  on  the  basis  of  the  precipitation  depth, 
elevation,  air  and  dew  point  temperatures,  albedo, 
radiation,  and  wind  speed.  Snowmelt  options  include  the 
temperature  index  method  or  the  energy  budget  method. 

Input  includes  the  precipitation  depths;  the  watershed- 
runoff  indexes  for  subdividing  flow  among  the  three 
processes;  initial  reservoir  elevations  and  outflows; 
drainage  areas;  bounds  on  usable  storage  and  allowable 
discharge  from  reservoirs;  total  computation  periods; 
routing  intervals;  and  other  special  instructions  to  control 
plots,  prints,  and  other  input-output  alternatives. 

5.2.1.  Rainfal 1 -Runoff  Relationships  for  PMF  Estimates 

Rain  falling  on  a  watershed  must  either  be  (a)  runoff, 
(b)  retained  in  the  soil  system,  (c)  intercepted  and 
evaporated  from  trees  and  other  plants,  (d)  evaporated  from 
pond,  lake,  and  stream  surfaces,  (e)  evaporated  from  the 
soil,  (f)  returned  to  the  atmosphere  by  transpiration  from 
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trees  and  other  plants,  or  (g)  percolated  into  the 
groundwater  system  and  thereby  lost  to  the  surface  water 
system.  If  percolation  to  the  groundwater  is  considered 
negligible;  this  rainfall  input  can  be  grouped  into  three 
classes:  (1)  runoff,  (2)  soil  moisture  increases,  and  (3) 
evapotranspi rat  ion  losses. 

The  percentage  of  water  available  for  runoff  is  found 
from  empirically  derived  relationships  of  the  Soil  Moisture 
Index  (SMI)  versus  Runoff  Percent  ( ROP ) ,  with  intensity  as  a 
third  variable.  In  the  SSARR  Model  this  is  entered  into  the 
program  in  the  form  of  a  table  for  each  watershed.  These 
values  are  usually  developed  by  trial  and  error.  Before  a 
design  curve  for  PMF  estimates  is  provided,  two  examples  of 
SMI-ROP  relationships  are  cited.  The  first  example  is  a 
three  variable  SMI  relationship  developed  for  the  Bird  Creek 
Basin  near  Sperry,  Oklahoma  (Figure  5.2)  (U.S.  Army  Corps  of 
Engineers,  1972).  This  basin  is  character ized  by  hot  summer 
temperatures,  long  dry  periods,  high  evapotranspi rat  ion 
rates,  and  high  intensity,  shor t -durat i on  rainfall.  The  SMI- 
ROP  relationship  was  developed  by  trial  and  error.  For 
saturated  soil  conditions  (equivalent  to  SMI  greater  than 
100  mm  (4  in.))  a  large  percentage  of  the  precipitation  is 
runoff;  for  very  low  precipitation  intensity  (0.25  mm  (0.01 
in.)  per  hour)  about  30%  is  runoff;  while  for  high- i ntensi ty 
values  (25.4  mm  (1.0  in.)  per  hour  and  greater)  the  runoff 
is  100%. 

The  second  example  is  the  relationship  calibrated  by 
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Figure  5.2.  Three  Variable  Soil  Moisture  Index  Relationships 
Developed  for  Bird  Creek  Basin  (After  U.S.  Army  Corps  of 

Engineers ,  1 972 ) . 

Alberta  Environment  (1980)  for  the  Little  Red  Deer  River 
watershed  in  the  Red  Deer  River  Basin,  as  illustrated  in 
Figure  5.3.  This  watershed,  because  it  is  located  much 
farther  north  than  the  basin  in  the  previous  example,  is 
characterized  by  a  moderate  climate  with  moderate  summer 
temperatures,  intermediate  evapotranspi rat  ion  rates,  and 
low- i ntens i ty  duration  rainfall.  The  main  difference  between 
these  two  examples  is  that  in  the  latter  case  the  lower 
intensity  of  the  rainfall  would  produce  a  greater  runoff, 
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Figure  5.3.  Three-variable  Soil  Moisture  Index  Relationship 
for  the  Little  Red  Deer  River  and  a  Design  Curve  for  PMP 
Loading  (After  Alberta  Environment,  1980). 

probably  because  of  the  lack  of  long  dry  periods  in  central 
Alberta . 

The  SMI  is  an  indicator  of  the  relative  soil  wetness 
and  is  used  to  determine  runoff.  When  the  soil  moisture  is 
depleted  to  approximately  the  permanent  wilting  point,  the 
SMI  is  a  small  number  associated  with  little  or  no  runoff. 
When  precipitation  recharges  soil  moisture,  the  value  of  the 
SMI  increases  until  it  reaches  a  maximum  value  considered  to 
represent  its  field  capacity  or  water  holding  capacity.  At 
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this  value  the  runoff  would  approach  100%.  The  SMI  is 
depleted  only  by  the  evapotranspi rat  ion  index  (ETI).  The  ETI 
can  be  specified  in  table  form  as  monthly  versus  average 
daily  potential  evapotranspi rat  ion  (in  inches  per  day). 
Maximum  runoff  is  characterized  by  a  zero  value  of 
evapotranspi rat  ion ,  which  means  there  would  be  no  losses  to 
evaporation  or  transpiration. 

For  PMF  estimates  the  maximum  runoff  is  desired,  and 
this  can  be  achieved  by  assuming  (a)  that  the  soil  is 
saturated  and  (b)  that  the  precipitation  occurs  at  a  high 
intensity.  To  represent  these  conditions,  all  SMI  values  are 
assigned  an  ROP  value  of  100%  and  ETI  is  given  a  value  of 
zero.  This  is  equivalent  to  the  relationship  shown  by  the 
solid  line  (labelled  PMP  loading)  in  Figure  5.3. 

5.2.2  Flow  Separation 

In  the  SSARR  Model,  the  rainfall  that  will  run  off  is 
divided  into  three  parts:  a)  base  flow,  b)  subsurface  or 
interflow,  and  c)  surface  runoff.  Because  of  antecedent 
conditions  (i.e.,  saturated  soil),  the  maximum  runoff  for 
PMF  estimates  would  occur  when  there  were  no  base  and 
subsurface  flows.  This  is  equivalent  to  the  1  -  to- 1  broken 
line  in  Figure  5.4  and  is  specified  for  the  SSARR  Model  in 
the  form  of  a  table.  The  solid  Tine  depicts  a  calibrated 
curve  that  has  been  used  in  the  Red  Deer  River  Basin  for 
natural  flows  during  the  summer  (Alberta  Environment,  1980). 
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Figure  5.4.  Surface-Subsurface  Separation  Curve  (After 

Alberta  Environment,  1980). 

5.2.3  Routing  of  Surface  Flow 

In  the  SSARR  Model,  each  component  of  surface  flow  is 
computed  as  input  rates,  expressed  in  inches  per  period. 

Each  period  value  is  converted  to  the  equivalent  inflow  rate 
in  cfs,  based  on  the  drainage  area  and  the  length  of  the 
period  in  hours.  Each  component  of  inflow  is  routed  through 
a  specified  number  of  increments  of  storage.  These 
increments  can  be  considered  a  series  of  small  lakes  that 
delay  runoff.  Routing  is  accomplished  by  solution  of  the 
basic  storage  equation  in  finite  time  periods  for  each 
increment  of  storage  routing  (Section  4.2.6).  Inflows  are 
computed  as  mean  values  for  the  period  for  watershed 
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routing.  The  time  of  storage  per  increment  can  be  a  few 
hours  to  a  month  or  more.  The  routing  specifications 
determine  the  time  delay  to  the  outflow  point  of  the 
watershed  and  the  shape  of  the  outflow  hydrograph.  In  this 
model,  virtually  unlimited  flexibility  in  time  distribution 
of  runoff  can  be  achieved  by  the  incremental  storage  routing 
technique.  The  lag  time  and  unit  peak  relationships 
available  in  the  SSARR  Model  for  various  combinations  of 
mu  1 1 i - i ncrement  routing  coefficients  are  illustrated  in 
Figure  5.5.  The  intersection  of  the  number  of  phases  and  the 
time  of  storage  per  phase  in  this  figure  give  the  lag  time 
and  the  peak  discharge.  In  this  model,  the  time  of  storage 
per  phase  and  the  number  of  phases  need  to  be  determined. 

The  first  of  these  two  parameters  is  usually  computed  by  an 
empirical  relationship  similar  to  Equation  4.9,  while  the 
second  parameter  can  be  obtained  through  calibration  of  the 
mode  1 . 

The  typical  routing  specifications  for  surface  flow 
that  have  been  calculated  for  a  539  km2  (208  sq .  mi.) 
watershed  in  Oregon  are  four  phases  and  2.5  hours  storage 
per  phase  (U.S.  Army  Corps  of  Engineers,  1972).  For  the  Red 
Deer  River  watershed  at  Sundre,  the  values  that  have  been 
calibrated  and  calculated  for  a  2471  km2  (954  sq .  mi.) 
watershed  are  two  phases  and  16  hours  storage  per  phase 
(Alberta  Environment,  1980). 


. 


134 


Figure  5.5.  Watershed  Multiple  Phase  Storage  Routing 
Available  in  SSARR  (After  U.S.  Army  Corps  of  Enqineers 

1972). 


5.2.4  Runoff  from  Snowmelt 

Runoff  from  snowmelt  may  contribute  significantly  to 
streamflow  in  Alberta,  and  hence  an  applicable  model  for 
this  province  must  consider  this  component.  In  the  SSARR 
Model,  extensive  routines  are  programmed  for  computation  of 
runoff  from  snowmelt.  Because  snowmelt  runoff  poses  complex 
problems  for  the  hydrologist,  this  aspect  has  been  ignored 
by  many  modellers  in  their  design. 

Unlike  rainfall,  snowmelt  is  not  generally  measured 
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quantitatively;  it  must  be  estimated  indirectly  from 
observations  of  meteorological  parameters  (e.g.,  air 
temperature,  depth  of  snow,  wind  speed)  to  determine  the 
snowmelt  rate  and  the  water  equivalent  of  the  snowpack. 

In  the  SSARR  Model  snowmelt  can  be  calculated  by  use  of 
the  generalized  snowmelt  equation  for  a  partly  forested  area 
(U.S.  Army  Corps  of  Engineers,  1956).  This  equation  requires 
measurements  of  air  temperature,  dewpoint  temperature,  wind 
speed,  insolation  (solar  radiation  on  horizontal  surface), 
snow  surface  albedo,  and  forest  canopy  cover,  for  estimating 
the  snowmelt  rate.  The  SSARR  Model  has  been  modified  by 
Alberta  Environment  (1977b)  to  perform  sp 1 i t -watershed 
computations,  that  is,  the  snow-covered  and  snow-free  areas 
are  treated  as  two  separate  watersheds,  each  with  its  own 
characteristics  and  parameters.  Since  the  number  of 
watersheds  is  doubled,  this  approach  is  more  cumbersome  and 
requires  more  data  for  simulation;  however,  it  does  allow 
for  separate  accounting  of  the  SMI  in  the  snow- free  and 
snow-covered  areas.  The  various  hydrologic  options  available 
in  the  split  watershed  approach  are  illustrated  in  Figure 
5.6.  These  options  are  specified  for  the  model  in  either 
tabular  form  or  as  part  of  the  data. 

5.2.5  Channel  Routing 

The  time  rate  of  change  of  streamflow  in  a  river  reach 
is  evaluated  by  first  dividing  the  reach  into  a  series  of 
increments.  The  increments  should  be  small  enough  so  that 
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Figure  5.6.  Split  Watershed  Option  of  the  Snow  Cover 
Depletion  Method  Available  in  the  SSARR  Model  (After  U.S. 

Army  Corps  of  Engineers,  1972). 

wedge  1  storage  is  negligible  compared  to  prismatic  storage. 
Inflows  are  computed  at  point  values  in  time  for  channel 
routing.  Outflow  from  each  increment  serves  as  inflow  to  the 
next  downstream  increment.  Time  of  storage  for  channel 
routing  may  vary  inversely  or  directly  as  a  power  function 
of  discharge  and  is  expressed  as  the  following  relationship: 


Linsley  et  al .  ( 1949)  showed  that  the  storage  in  a  channel 
reach  may  be  considered  to  be  the  sum  of  two  portions:  prism 
storage,  or  the  water  below  an  imaginary  line  drawn  parallel 
to  the  channel  bottom;  and  wedge  storage,  or  the  water 
between  that  line  and  the  actual  water  surface  profile. 
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rp  _  KTS 

s  (5.1) 

where  Tg  is  the  time  of  storage  per  increment  in  hours,  KTS 
is  a  constant  determined  by  trial  and  error  or  estimated 
from  physical  measurements  of  flow  and  cor respondi ng  routing 
times,  Q  is  the  discharge  in  cfs,  an  n  i s  a  coefficient 
usually  between  -1  and  +1.  A  value  of  n  =  0.20  was  found 
reasonable  for  most  streams  in  the  Columbia  Basin  (U.S.  Army 
Corps  of  Engineers,  1972);  n  =  0.35  was  found  for  the  Mekong 
River  (U.S.  Department  of  Commerce,  1970);  and  n  =  0.13  was 
determined  for  the  Red  Deer  River  Basin  (Alberta 
Environment,  1980).  Although  the  last  appears  to  be  low  with 
respect  to  the  other  two  values,  it  has  been  obtained  by 
Alberta  Environment  from  calibration  for  the  Red  Deer  River 
Bas i n . 

Routing  of  flow  is  accomplished  by  an  iterative 
solution  of  the  storage  equation.  A  backwater  mode  exists 
for  cases  in  which  elevation  and  discharge  are  affected  by 
backwater  from  a  downstream  time-variant  source.  Examples  of 
such  occurrences  are  river  estuaries  affected  by  tidal 
fluctuations,  river  reaches  upstream  from  a  junction  with  a 
major  tributary,  and  the  upstream  reaches  of  a  reservoir  or 
a  lake  whose  outflow  is  affected  by  the  elevation  of  another 
lake  just  downstream.  This  model  utilizes  a  three-var i able 
relationship  between  upstream  elevation,  downstream 
elevation  or  flow,  and  discharge  from  the  upstream  location. 

The  above  components  are  organized  and  coordinated  to 
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produce  a  river  basin  model  that  simulates  the  hydrologic 
response  of  the  physical  system.  The  system  configuration 
essentially  describes  to  the  computer  the  physical  layout 
and  relationships  of  all  components  of  the  system.  The 
configuration  is  in  upstream  to  downstream  order  of  all 
watersheds,  lakes,  reservoirs,  channel  reaches,  and 
confluence  points  for  a  particular  basin.  The  hierarchical 
procedure  involves,  first,  watershed  routing,  then 
consecutive  channel  routing,  and  finally  combining  until  all 
operations  are  complete.  The  computer  model  then  integrates 
the  two  hierarchies  and  performs  the  model  simulation  of  the 
river  basin  system. 

An  example  of  a  configuration  for  the  watersheds  of  the 
Reo  Deer  River  Basin  from  the  headwaters,  through  Sundre,  to 
Red  Deer  is  depicted  in  Figure  5.7.  In  this  figure,  four 
watersheds  are  simulated:  headwater  watershed  to  Sundre 
(drainage  area  of  2471  km2  or  954  sq .  mi.),  Little  Red  River 
watershed  (drainage  area  of  2392  km2  or  924  sq .  mi.),  James 
River  watershed  (drainage  area  of  821  km2  or  317  sq.  mi.), 
and  Medicine  River  watershed  (drainage  area  of  655  km2  or 
253  mi2).  Because  of  the  snowmelt  computations,  each  of  the 
above  watersheds  is  treated  as  two  separate  watersheds  (one 
for  the  snow-free  area  and  the  other  for  the  snow-covered 
area).  For  example,  in  Figure  5.7  the  number  301  is  assigned 
to  the  snow-covered  area  for  the  Little  Red  Deer  River 
watershed,  while  the  number  302  refers  to  the  snow-free  area 
for  the  same  watershed. 
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Headwater  Basin 
Sundre/Red  Deer 
(95^  sq.mi.) 


Red  Deer  Basin 
Outflow 


Figure  5.7.  A  Watershed  Configuration  Used  in  the  SSARR 
Model  for  the  Red  Deer  River  Basin. 


5.2.6  PMF  Estimates  from  the  SSARR  Model 


The  various  parameters  and  relationships  discussed 
previously  are  first  maximized  and  then  used  in  the  SSARR 
Model  to  produce  a  PMF  hydrograph  and  a  peak  discharge.  In 
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the  above  computations  the  soil  is  considered  saturated, 
with  no  evapotranspiration  losses,  and  all  available  water 
is  translated  as  surface  runoff.  An  underlying  assumption  is 
that  the  surface  area  of  precipitation  and  the  area  of  the 
watershed  are  interchangeable.  Usually  these  two  quantities 
are  different,  however,  under  maximum  conditions  they  are 
commonly  assumed  to  be  similar.  Hence  these  terms  can  be  and 
are  used  interchangeably  in  PMP  and  PMF  studies. 

The  computer  program  of  the  SSARR  Model  allows  the 
results  of  the  simulation  to  be  produced  either  in  tabular 
or  graphical  form.  Because  the  tabular  results  are 
extensive,  graphical  depiction  is  usually  desirable.  A 
graphical  example  of  a  PMF  hydrograph  for  the  Red  Deer  River 
Basin  at  Sundre  is  depicted  in  Figure  5.8.  The  PMF 
hydrograph  is  the  result  of  a  24-hour  precipitation  loading 
of  340  mm  (13.4  in.)  on  a  247  1  Km 2  (954  sq.  mi.)  drai nage 
area  at  Sundre.  The  precipitation  loading  is  distributed  at 
a  maximum  intensity  of  30.5  mm  (1.2  in.)  per  hour  for  11 
hours.  Snowmelt  is  calculated  using  the  generalized  snowmelt 
equation  for  a  partly  forested  area.  In  future  discussions 
of  the  watershed  at  Sundre  it  can  be  assumed  that  the  above 
values  of  precipitation  loading  apply.  Any  deviation  from 
these  values  will  be  mentioned  in  the  discussion.  The 
largest  recorded  hydrograph  at  Red  Deer  is  shown  on  the  same 
diagram  (Figure  5.8)  for  comparison.  The  Red  Deer  station 
was  used  because  it  is  the  only  station  with  long-term 
discharge  data  in  the  Red  Deer  River  Basin.  The  hydrograph 
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Figure  5.8.  PMF  Hydrograph  at  Sundre  and  the  Largest 
Recorded  Hydrograph  at  Red  Deer. 
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recorded  at  Red  Deer  (drainage  area  of  11  600  Km2  or  4480 
sq.  mi.)  occurred  from  26  June  to  2  July  1915  and  resulted 
in  a  maximum  instantaneous  discharge  of  1933  m3  per  second 
or  68  250  cfs  (a  ratio  of  0.17  m3  per  second  per  Km2  or  15.2 
cfs  per  sq.  mi.),  a  maximum  daily  discharge  of  1586  m3  per 
second  or  56  000  cfs  (a  ratio  of  0.14  m3  per  second  per  Km2 
or  12.5  cfs  per  sq .  mi.),  and  a  record  gauge  height  of  5.8  m 
(19.05  ft.)  at  Red  Deer.  Historical  reconstruct  ion  of  this 
record  rainstorm  suggests  that  probably  the  rainfall  (which 
is  unknown)  associated  with  this  event  also  produced  record 
hydrograph  flows  in  the  entire  Red  Deer  River  Basin,  but 
because  of  a  lack  of  recording  stations  in  the  basin  at  that 
time,  this  is  uncertain. 

The  PMF  results  can  be  summarized  by  comparing  the 

ratio  of  the  discharge  divided  by  the  drainage  area  (cfs  per 

sq .  mi.)  to  the  drainage  area  ( sq .  mi.),  as  is  shown  in 

Figure  5.9.  This  ratio  (discharge  divided  by  the  drainage 

area)  can  be  interpreted  as  normalized  discharge.  In  Figure 

5.9,  the  recorded  and  computed  results  are  compared  for  the 

Red  Deer  River  Basin.  The  largest  recorded  instantaneous 

flows  for  the  various  recording  stations  in  the  basin  are 

depicted  by  triangles.  These  flows  can  be  enveloped  by  a 

0.5 

relationship  of  the  form  Q  =  1000  A  ,  where  Q  is  the 
instantaneous  maximum  discharge  (cfs)  and  A  is  the  drainage 
area  ( sq .  mi.).  Also  shown  in  the  figure  are  the  PMF 
estimates  obtained  from  the  SSARR  Model  for  some  of  the 
watersheds  (displayed  by  stars). 
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5.2.7  Sensitivity  of  SSARR  Model. 

The  sensitivity  of  the  SSARR  Model  to  parameter  changes 
for  PMP  loading  are  discussed  next  in  terms  of  a  number  of 
hypothetical  simulations.  Because  the  model  allows  the  user 
considerable  flexibility  to  simulate  hydrologi ca 1 ly 
different  watersheds,  good  judgment  and  logic  are  essential 
to  properly  develop  parameters  and  relationships.  Various 
parameters  and  relationships  were  examined  for  the  watershed 
at  Sundre  on  the  Red  Deer  River.  The  sensitivity  of  the 
model  was  tested  by  holding  all  parameters  constant  at 
maximum  conditions  and  then  changing  only  the  desired 
parameter  or  relationship.  The  maximum  conditions  used  in 
the  simulations  are  a  24-hour  precipitation  loading  of  340 
mm  (13.4  in.)  (with  a  maximum  intensity  of  30.5  mm  (1.2  in.) 
per  hour)  on  a  2471  Km2  (954  sq .  mi.)  drainage  area.  Other 
watersheds  in  the  Red  Deer  River  Basin  were  also  tested,  but 
those  results  are  not  included  here  since  they  showed  the 
same  shape  and  peak  discharge  in  the  resultant  hydrograph  as 
the  watershed  at  Sundre. 

One  of  the  first  relationships  examined  was  the  three- 
parameter  SMI-ROP  relationship.  This  relationship  for 
natural  flows  determines  the  volume  of  runoff  and  affects 
the  shape  of  the  hydrograph.  Two  quantities  seem  to  be 
important:  the  SMI  value  and  the  average  intensity  of  the 
loading.  The  SMI  can  be  any  value  from  0  to  127  mm  (5  in.) 
depending  upon  the  antecedent  conditions.  The  upper  value 
was  chosen  as  an  arbitrary  upper  limit  for  the  tests 
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conducted  in  this  work  because  others  (U.  S.  Army  Corps  of 
Engineers,  1972;  Alberta  Environment,  1980)  have  used  this 
value.  For  PMP  loading,  because  saturated  soil  conditions 
are  assumed,  the  SMI  is  assigned  the  upper  limit  value  of 
127  mm  (5  in.).  Using  this  upper  value  a  number  of  different 
average  intensities  were  tested.  The  results  were  identical 
to  those  shown  in  the  PMF  hydrograph  depicted  in  Figure  5.8. 
This  was  because  of  the  allocation  of  the  total  flow  to 
surface  flow  in  the  surface-separation  relationship.  Hence, 
for  upper  limit  values  of  SMI,  the  SMI-ROP  relationship  has 
no  effect  on  PMF  estimates  in  the  SSARR  Model. 

A  second  relationship  that  was  tested  w'as  the  daily 
evapotranspi rat  ion .  Values  of  evapotranspi rat  ion  from  0  to 
10  mm  (0.4  in.)  per  day  were  used,  and  these  had  no  effect 
on  the  PMF  hydrograph.  Furthermore,  the  resultant 
hydrographs  were  identical  to  that  shown  in  Figure  5.8. 
Although  this  index  may  be  important  for  long-term 
simulations,  it  is  not  believed  to  be  critical  for  PMP 
loading . 

A  third  relationship  investigated  was  the  surface- 
subsurface  separation  curve.  For  PMF  estimates  this 
relationship  was  represented  by  the  1  -  to- 1  broken  line  in 
Figure  5.4  and  corresponded  to  saturated  soil  conditions  and 
maximum  runoff.  For  natural  flow  conditions,  such  a  curve 
would  rarely  be  attained.  Instead,  for  those  conditions  the 
expected  curve  would  be  similar  to  the  solid  line  curve 
shown  in  Figure  5.4  and  is  usually  determined  by 
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calibration,  through  a  tr i a  1 -and-error  procedure.  This 
relationship  can  be  examined  hypothetically  for  a  water 
loading  with  different  surface-subsurface  curves,  which  will 
provide  insight  into  the  effects  of  these  curves  on  the 
shape  and  peak  discharge  of  the  resultant  hydrograph.  The 
relationship  was  examined  by  decreasing  the  surface  curve 
with  each  simulation.  As  the  surface  component  of  the 
relationship  is  decreased,  the  subsurface  component  is 
increased,  since  the  total  sum  of  surface  and  subsurface  is 
constant.  Three  different  sur face-subsur f ace  curves  were 
tested,  and  the  resulting  hydrographs  are  illustrated  in 
Figure  5.10.  The  effects  of  varying  the  sur f ace-subsur f ace 
separation  relationship  on  streamflow  depends  on  the 
magnitude  of  the  surface  to  subsurface  components.  When  most 
of  the  flow  is  allocated  to  surface  and  little  to  the 
subsurface  flows,  a  higher  and  earlier  peak  results  in  the 
hydrograph.  The  largest  peak  in  magnitude  occurs  when  100% 
is  allocated  to  surface  flow  and  0%  to  subsurface  flow,  as 
is  the  case  for  PMP  loading  (Curve  A).  The  other  extreme  is 
allocating  100%  of  the  flow  to  the  subsurface  and  0%  to 
surface,  and  a  straight-line  hydrograph  results.  In  between 
these  two  extremes  are  the  curves  depicted  by  B  and  C, 
equivalent  to  what  may  be  called  unsaturated  soil 
conditions.  The  variation  of  the  sur face-subsurf ace 
separation  relationship  would  only  be  important  to  PMF 
estimates  if  unsaturated  soil  conditions  exist.  - 

In  addition  to  the  three  relationships  above,  a  number 
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Figure  5.10.  Effect  on  Streamflow  of  Varying  the  Surface- 
.  Subsurface  Separation  Relationships. 
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of  pouting  parameters  were  examined.  The  effects  on  the 
outflow  hydrograph  shape  and  timing  response  can  be  observed 
by  changing  the  number  of  routing  phases  for  each  component 
of  flow  and  the  time  of  storage  per  phase  (T  ).  Because  two 
parameters  are  involved  ( i .e. ,  the  number  of  routing  phases 
and  the  time  of  storage  per  phase) ,  there  are  a  number  of 
possibilities  for  examining  these  parameters.  One  way  to 
display  any  possible  effect  is  to  Keep  one  parameter 
constant  while  varying  the  second,  and  vise  versa.  The 
simulations  discussed  here  are  hypothetical,  since  the  two 
parameters  are  predetermi ned  by  computation  and  calibration, 
which  results  in  a  fixed  value  for  each  parameter.  The 
hypothetical  simulation  is  useful  in  that  it  provides 
insight  into  the  sensitivity  of  these  parameters  for  PMF 
est imates . 

In  the  first  of  these  simulations,  all  coefficients  and 
relationships  were  Kept  constant  while  the  number  of  routing 
phases  for  surface  flow  was  allowed  to  vary  from  the 
calibrated  value  of  2.  The  results,  displayed  in  Figure 
5.11,  indicate  two  effects  arising  from  an  increase  in  the 
number  of  routing  phases:  (a)  a  delay  in  the  time  of  peaKing 
and  (b)  a  decrease  in  the  peaK  discharge.  Thus  the  largest 
peaK  discharge  can  be  obtained  using  the  smallest  possible 
number  of  phases.  In  all  of  these  examples  the  total  volume 
of  runoff  was  the  same. 

In  the  second  simulation,  the  time  of  storage  per  phase 
(Tg)  was  varied  from  the  16-hour  calibrated  value.  This 
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Figure  5.11.  Effect  on  Streamflow  of  Varying  the  Number  of 

Surface  Routing  Phases. 
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variation  can  also  result  in  significant  changes  in  timing 
and  in  the  hydrograph  character i st i cs .  To  examine  this 
parameter,  the  surface  phase  time  was  varied  while  all  other 
parameters  were  Kept  constant.  The  variation  of  the  surface 
phase  time  is  illustrated  in  Figure  5.12,  where  two 
different  times  were  selected.  The  values  selected  varied  by 
50%  from  the  16-hour  surface  phase  time.  Sixteen  hours  was 
the  surface  phase  time  calibrated  for  a  2471  Km2  (954  sq . 
mi.)  watershed  area  in  the  Red  Deer  River  Basin  and  is  the 
value  suggested  for  PMF  estimates.  As  expected,  an  increase 
in  surface  phase  time  results  in  lower  peak  discharges  and 
greater  lag  times. 

The  last  parameter  examined  for  the  SSARR  Model  was  the 
coefficient  n  in  Equation  5.1.  The  value  of  n  (0.13)  was 
obtained  from  calibration  for  the  Red  Deer  River  Basin  and 
can  be  considered  a  fixed  parameter.  This  coefficient  was 
examined  for  n  =  0.26  to  determine  the  effect  on  the  shape 
and  peak  discharge  of  the  resultant  hydrograph  for  PMF 
estimates.  The  resultant  hydrographs  for  both  values  of  the 
coefficient  were  nearly  the  same  in  shape  and  in  peak 
discharge,  suggesting  that  changes  only  in  this  parameter 
are  insignificant  for  PMF  estimates. 

In  summary,  for  a  given  watershed  area  and  mass 
precipitation  distribution  with  the  same  intensity  loading 
and  for  saturated  soil  conditions,  two  parameters  affect  the 
PMF  peak  discharge  for  PMP  loading:  the  number  of  surface 
routing  phases  and  the  time  of  storage  per  phase.  These  two 
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Figure  5.12.  Effects  on  Streamflow  of  Surface  Time  of 

Storage . 
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parameters  are  interrelated  (Figure  5.5),  and  an  effect  on 
one  produces  an  effect  on  the  other.  The  other  relationships 
examined  have  little  effect  on  PMF  estimate  in  the  SSARR 
Model.  The  values  used  and  suggested  for  a  PMF  estimate  for 
a  24-hour  precipitation  loading  of  340  mm  (13.4  in.)  on  a 
2471  Km2  (954  sq .  mi.)  drainage  watershed  at  Sundre  in  the 
Red  Deer  River  Basin  are:  (1)  SMI  equal  to  127  mm  (5  in.) 
for  all  rainfall  intensities,  (2)  evapotranspi rat  ion  value 
of  0,  (3)  100%  of  the  total  flow  to  be  allocated  to  surface 
flow,  (4)  number  of  routing  phases  is  2,  (5)  the  time  of 
storage  per  phase  is  equal  to  16  hours,  and  (6)  the 
coefficient  n  is  equal  to  0.13. 

5.3.0  HYMO  Model 

HYMO  (Williams  and  Hann,  1973)  is  a  prob 1 em-or i ented 
computer  language  consisting  of  a  main  program  and  16 
subroutines  that  when  combined  for  watersheds  result  in  a 
hydrologic  model,  hence  the  name  HYMO.  This  model  can  be 
used  for  planning  flood  prevention  projects,  forecasting 
floods,  and  research  studies.  The  HYMO  language  was  designed 
to  transform  rainfall  data  into  runoff  hydrographs  and  to 
route  these  hydrographs  through  streams  and  valleys  or 
reservoirs.  The  procedures  used  in  this  model  were  selected 
by  Williams  and  Hann  (1973)  because  of  their  accuracy, 
general  applicability,  practicality  of  inputs,  and 
computational  efficiency. 

The  HYMO  language  is  based  on  hydrologic  techniques 
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suggested  by  the  Soil  Conservation  Service  (SCS,  1969)  for 

the  prediction  of  flood  hydrographs  on  ungauged  basins.  The 

main  features  in  this  model  include  (a)  ability  to  subdivide 

the  drainage  area  into  various  subareas,  (b)  ability  to 

account  for  the  effect  of  various  land  uses,  (c)  ability  to 

estimate  the  time  of  concentration  for  each  subwatershed, 

which  in  turn  is  a  function  of  the  flood  travel  time,  (d) 

ability  to  account  for  the  effects  of  channel  storage  in 

each  river  reach,  and  (e)  ability  to  produce  flood 

hydrographs  at  any  point  in  the  watershed.  The  main 

conceptual  and  computational  features  include  the  peak  flow 

rates  (q  )  and  the  unit  hydrograph.  The  peak  flow  rates  from 
1? 

individual  watershed  subareas  can  be  computed  by  the 
equation 


where  qp  is  the  peak  flow  rate,  B  is  a  statistically  based 
watershed  parameter  computed  within  the  model  as  a  function 
of  watershed  characteristics,  A  is  the  watershed  area  ( sq . 
mi.);  Q  is  the  volume  of  runoff  (in.),  and  tp  is  the  time  to 
peak  (hours).  The  unit  hydrographs  are  divided  into  three 
parts  for  computation  (Figure  5.13). 

In  the  first  part  (from  the  beginning  of  the  rise  to 
the  inflection  point,  tQ),  the  hydrograph  is  computed  by  the 
two-parameter  gamma  distribution  equation 
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Figure  5.13.  Dimensionless  Unit  Hydrograph  Used  in  the  HYMO 
Model  (After  Williams  and  Hann,  1973). 


.  n-1  (l-n)(t  /  t  -  1) 

q  =  q_  (-f)  e  P  (5-3l 

P  XP 

where  q  is  the  flow  rate  at  time  t,  q  is  the  peak  flow 

P 

rate,  t  is  the  time  to  peak  (hours),  and  n  is  a 
dimensionless  parameter  determined  from  watershed 
character i s  t ics . 

For  the  second  part,  from  t0  to  t  ^ ,  the  recession 
depletion  equation  is  used  to  compute  the  hydrograph: 


q 


e(to  -  /  K 


(5.4) 
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where  q  is  the  flow  rate  at  the  inflection  point,  t0  is  the 
time  at  the  inflection  point  (hours),  and  K  is  the  recession 
constant  (hours)  such  that  t^  =  tQ  +  2K. 

In  the  third  part,  from  ti  tooo,  the  recession 
depletion  equation  becomes 

q  =  q1  e(tl  “  /  K1  (5.5) 

where  q  is  the  flow  rate  at  t  ,  and  K  =  3K  is  a  second 
1  11 

recession  constant  (hours).  The  dimensionless  shape 

parameter,  n,  is  a  function  of  K/t  ,  as  shown  in  Figure 

P 

5.14;  therefore,  the  entire  unit  hydrograph  can  be  computed 
if  K  and  tp  are  Known.  To  compute  K  and  tp  for  ungauged 
watersheds,  HYMO  uses  the  equations 

K  =  16.1  A0,24  s-0.84  (5.6) 


and 


tp  =  6.54  A0*39  S'0’-50  (5.7) 

where  A  is  the  watershed  area  ( sq .  mi.),  and  S  is  the 
difference  in  elevation  in  feet,  divided  by  flood-plain 
distance  in  miles,  between  the  watershed  outlet  and  the  most 
distant  point  on  the  watershed.  Calculating  the  peak  flow 
using  equation  5.2,  the  volume  of  runoff  (Q)  in  inches  is 
computed  by  utilizing  the  SCS  rai nf a  1 1 -runoff  relationship 
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Figure  5.14.  Relationship  Between  Dimensionless  Shape 
Parameter  n  and  Recession  Constant/Time  to  Peak  (After 

Williams  and  Hann,  1973). 


(Soil  Conservation  Service,  1969).  This  can  be  obtained  by 
the  following  equations 


(  for  P >  0 . 25)  (5.8) 


and 


(5.9) 


where  P  is  the  rainfall  (inches),  S  is  the  rainfall 
retention  parameter,  and  CN  is  the  soil -cover  complex  number 
determined  from  actual  soil  types  and  land  cover  for  each 
subwatershed.  The  SCS  rai nf a  1 1 -runoff  relationship  is 
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expressed  by  a  set  of  numbered  curves  (Figure  5.15).  In  the 


Figure  5.15.  Graphical  Solution  of  Rainfall -runoff  Equation 
(After  Soil  Conservation  Service,  1969). 


SCS  method  the  change  in  S  (actually  in  CN)  is  based  on  the 
soil  antecedent  moisture  condition,  determined  by  the  total 
rainfall  preceding  the  storm.  For  best  estimates  of  peak 
flood  flows,  the  watershed  should  be  subdivided  into 
subcatchment  areas,  each  of  which  is  homogeneous  with 
respect  to  soil  type,  topography,  vegetation  cover,  land 
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use,  and  drainage  pattern. 

The  variable  storage  coefficient  (VSC)  flood-routing 
method  (Williams,  1969)  is  used  in  the  HYMO  method.  This 
method  uses  an  iterative  solution  and  is  free  of  convergence 
problems.  The  VSC  routing  equations  are 


(5.10) 


2  At 


C 


1 


2  T±  +  At 


(5.11) 


C 


2  At 


2 


2  T2  +  At 


(5.12) 


T1  =  (  1800  (V,  +  Vn  ) 

1 1  U1 


L 


L  x  SLPq  1/2 


(5.13) 
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1800  (V  +  V  ) 

12  U2 
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(5.14) 


159 


where  I  is  the  inflow  rate  (cfs),  0  is  the  outflow  rate 
(cfs),  I  is  the  average  inflow  rate  (cfs),  C  is  the  storage 
coefficient,  T  is  the  travel  time  through  the  reach  (hours), 
L  is  the  reach  length  (ft.),  V  is  the  velocity  (ft.  per 
second),  SLPo  is  the  normal  slope,  and  D  is  the  depth  (ft.). 
In  the  above  equations  the  subscripts  1  and  2  refer  to  the 
beginning  and  end  of  the  time  interval,  At.  Subscripts 
II  , 1 2  iOi  »  and  0  2  are  also  used  in  conjunction  with  the 
normal  depth,  D,  and  velocity,  V,  to  define  the  spatial  and 
temporal  values  of  these  quantities.  Since  T2  and  C2  are 
dependent  upon  0^  ,  an  iterative  technique  is  required  to 
solve  the  routing  equations.  In  Equation  5.10,  I  and  0^ 
are  Known,  and  C  can  be  computed  from  Equation  5.12.  Hence 
only  O2  and  C2  are  unknown.  0^  can  be  used  as  a  first 
approximation  of  0^  .  The  normal  depth  and  velocity  for  the 
approximate  value  of  0^  are  entered  into  Equation  5.14  for 
computing  .  Then  Equation  5.11  is  used  to  compute  C2  • 

The  second  approximation  of  O2  is  obtained  from  Equation 
5.10.  This  iterative  process  continues  until  there  is  a 
difference  of  about  0.1%  or  less  between  successive  O2 
values . 

HYM0  uses  the  storage- indication  method  to  route  floods 
through  reservoirs.  This  method  has  been  widely  used  and 
accepted  because  it  is  practical  and  accurate.  The  model 
requires  rating  curves  along  a  valley  to  adequately  describe 
the  hydraulics  of  the  stream  and  valley.  The  model  also  uses 
Manning's  equation  (Viessman  et  at.,  1977)  to  compute  the 
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normal  flow-rating  curves  that  are  used  in  the  VSC  flood¬ 
routing  method.  A  normal  flood  plain  slope  is  determined  for 
each  valley  section  by  plotting  a  profile  of  the  flood 
plain.  The  normal  channel  slope  is  determined  by  plotting  a 
profile  of  the  flood  plain  with  channel  distances.  Since 
flow-rating  curves  are  needed  in  this  model  but  are  not 
always  available,  this  tends  to  make  the  HYMO  Model  less 
effective.  Also,  the  user  needs  to  describe  the  valley 
section  to  the  computer,  which  may  not  always  be  possible. 
This  model  has  an  advantage  in  that  it  is  a  user -or i ented 
model,  with  a  smaller  data  requirement  than  the  SSARR  Model. 
Because  of  its  simplicity,  however,  it  lacks  powerful 
methods  for  sophisticated  hydrological  simulation,  including 
snowmelt  computations. 

5.3.1  PMF  Estimates  from  the  HYMO  Model 

A  PMF  estimate  can  be  computed  with  the  HYMO  Model  by 
(1)  assigning  the  curve  number  (CN)  a  value  of  100,  which 
represents  saturated  soil  conditions  and  (2)  distributing 
the  PMP  loading  to  produce  the  maximum  discharge.  Snowmelt 
can  be  simulated  by  adding  the  melt  rate  to  the 
precipitation  mass  distribution.  The  precipitation 
distribution  can  be  specified  in  the  model  by  the  time 
period  (DT)  and  the  mass  distribution.  To  investigate  the 
choice  of  the  time  period,  three  different  time  periods  were 
tested  for  a  24-hour  water  loading  of  645  mm  (25.4  in.)  over 
a  2471  km2  (954  sq.  mi.)  watershed  area  with  a  slope  of  5.69 
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m  per  Km  (30.1  ft.  per  mi.)  in  the  Red  Deer  River  Basin.  The 
water  loading  consisted  of  a  PMP  loading  of  340  mm  (13.4 
in.)  and  a  snowmelt  of  305  mm  (12  in.).  The  precipitation 
loading  of  340  mm  (13.4  in.)  is  distributed  at  a  maximum 
intensity  of  30.5  mm  ( 1 . 2  in.)  per  hour  for  11  hours.  The 
resultant  PMF  hydrographs  are  depicted  in  Figure  5.16  for  DT 
=1,3,  and  6  hours.  For  the  different  DT  values  the  results 
show  only  a  slight  difference  in  the  shape  of  the  hydrograph 
and  in  the  peak  discharge.  Hence,  using  a  time  period  of  1, 

3 ,  or  6  hours  will  result  in  comparable  hydrographs  and  peak 
discharges.  Because  the  largest  peak  discharges  were 
obtained  with  the  1-  and  3-hour  time  periods,  these  periods 
were  used  in  analysis  in  this  work. 

5.3.2  Sensitivity  of  the  HYMO  Model 

In  the  HYMO  Model  a  number  of  parameters  are  needed  for 
simulation  of  the  PMF  hydrograph,  and  these  quantities  are 
either  recorded,  measured,  or  calculated.  For  a  particular 
watershed,  specific  values  of  these  parameters  are  obtained 
that  are  called  fixed  quantities,  that  is,  they  do  not  vary 
within  measuring  tolerance.  Even  though  these  parameters  are 
fixed,  it  is  still  possible  to  examine  the  model's 
sensitivity  to  parameter  variation  by  considering 
hypothetical  simulations.  These  simulations  are  advantageous 
in  that  they  can  provide  insight  into  variations  of  the 
magnitude  of  the  peak  discharge.  The  same  watershed  as  was 
described  in  sections  5.2.6  and  5.3.1  is  used  here  for  the 
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Figure  5.16.  Variation 
Water  Loading  of  645  mm 
mi  .  )  Watershed  Area 


of  the  Time  Period  for  a  24-hour 
25.4  in.)  over  a  2471  Km2  (954  sq . 
in  the  Red  Deer  River  Basin. 
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hypothetical  simulations. 

The  first  parameter  is  the  curve  number.  The  HYMO  Model 
uses  the  SCS  rai nfal 1 -runoff  relationship  (Soil  Conservation 
Service,  1969)  to  obtain  an  estimate  of  the  direct  runoff, 
and  hence  an  approximation  of  CN  is  required  in  the  model.  A 
value  of  CN  =  100  gives  the  maximum  direct  runoff  and 
corresponds  to  the  soil  being  saturated.  For  natural  flows, 
CN  is  usually  lower  than  100;  therefore,  it  is  important  to 
investigate  the  effect  of  using  lower  values  than  CN  =  100. 
The  effect  is  examined  by  using  three  different  CN  numbers 
(100,  80,  and  65,  corresponding  to  curves  A,  B,  and  C  in 
Figure  5.17).  Curves  A,  B,  and  C  show  that  the  peak 
discharge  decreases  with  a  decrease  in  the  CN  value 
(equivalent  to  progressing  from  saturated  to  unsaturated 
conditions).  The  maximum  peak  discharge  occurs  with  CN  = 

100;  hence,  this  value  is  used  for  the  other  simulations  of 
the  HYMO  Model . 

The  second  parameter  investigated  was  the  average  slope 
of  the  watershed,  represented  by  K  and  t  .  These  two 
quantities  can  be  computed  by  equations  5.6  and  5.7  in  the 
HYMO  Model.  For  any  watershed,  the  average  slope  is  a 
predetermi ned  fixed  value.  It  is  possible  to  investigate  the 
effect  of  this  fixed  value  on  the  PMF  estimate  by 
considering  hypothetical  simulations  in  which  all  parameters 
are  kept  constant  while  the  average  slope  is  varied  with 
each  simulation.  To  test  any  possible  effect,  two  average 
slopes  (ranging  from  0.95  to  5.69  m  per  km)  were  simulated 
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Figure  5.17.  Effect  on  the  Resultant  Hydrograph  by 
Variation  of  the  Curve  Number  in  HYMO  Model. 
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for  the  2471  Km2  (954  sq .  mi.)  watershed  with  average  slope 
of  5.69  m  per  Km  (30.1  ft.  per  mi.)  in  the  Red  Deer  River 
Basin.  The  2471  Km2  watershed  has  one  of  the  steepest 
average  slopes  of  the  watersheds  in  the  Red  Deer  River 
Basin.  Thus  the  average  slopes  selected  in  the  simulation 
are  arbitrary  lower  values  than  the  estimated  value  of  5.69 
m  per  Km  for  this  watershed.  The  effect  of  the  hypothetical 
simulations  is  depicted  in  Figure  5.18.  The  resultant 
hydrograph  shows  that  the  largest  peaK  discharge  coincides 
with  the  largest  average  slope,  and  similarly  the  smallest 
peaK  discharge  coincides  with  the  smallest  average  slope. 
Increases  in  the  average  slope  produce  increases  in  the  peaK 
discharge.  This  figure  shows  that  in  the  HYMO  Model,  average 
slope  not  only  affects  the  shape  but  also  the  peaK  of  the 
hydrograph.  Two  watersheds  with  similar  areal 
characteristics  but  different  slopes  therefore  would  not 
have  the  same  peaK  discharge. 

The  third  parameter  was  the  area  of  the  watershed.  The 
watershed  area,  similar  to  the  average  slope,  is  a 
predetermined  fixed  value  in  the  HYMO  Model.  For  this 
parameter  it  was  desirable  to  investigate  the  effect  the 
watershed  area  might  have  on  the  PMF  estimate.  In  this 
investigation  the  watershed  at  Sundre  was  again  used,  and 
two  different  watershed  areas  were  hypothetically  simulated, 
with  all  other  conditions  being  the  same.  The  results,  as 
shown  in  Figure  5.19,  i nd i ca te  a  1 ogar i t hmi c  i ncrease  i n 
area  with  peaK  discharge.  Two  watersheds  with  similar 
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Figure  5.18.  Effect  of  Average  Slope  of  the  Watershed  on  the 

Resultant  Hydrograph. 
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watershed  parameters  but  of  different  areas  would  also  have 
different  peak  discharge. 

In  summary,  for  a  watershed  with  the  same  intensity 
loading,  three  parameters  can  significantly  affect  the  PMF 
peak  discharge  in  the  HYMO  Model:  the  curve  number,  average 
slope,  and  the  watershed  area.  For  PMF  estimates,  all  three 
of  these  parameters  are  fixed  quantities,  and  once 
determined,  they  do  not  change.  Because  these  fixed 
quantities  affect  the  resultant  PMF,  it  is  important  that 
they  are  determined  as  accurately  as  possible. 


F i gure 


id 


co 


o 


Average  Slopet  30*1  feet  per  mile 
Curve  Number  i  100 


AREA  (SO.  Ml.) 


5.19.  Area  Versus  Peak  Discharge  Relationship  for  PMF 
Estimates  in  the  HYMO  Model. 
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CHAPTER  6 

COMPARISON  OF  MODELS 

In  the  last  two  chapers,  three  models  (TAPMF,  SSARR, 
and  HYMO)  were  presented  for  estimating  the  PMF  for  gauged 
watersheds.  An  objective  in  this  study  was  to  compare  the 
results  from  these  models  in  order  to  determine  which  model 
would  be  the  most  appropriate  for  estimating  the  peak  PMF 
discharge  in  Alberta.  The  PMF  can  be  regarded  as  an  extreme 
event  and  probably  has  not  yet  occurred,  hence  only 
comparisons  between  model  results  can  be  made.  A  graphical 
approach  was  selected  for  comparing  the  results.  In  such 
comparisons,  the  normalized  discharge  (expressed  as  m3  per 
second  per  km2  or  cfs  per  sq .  mi.)  is  usually  plotted 
against  the  drainage  area  of  the  watershed  (km2  or  sq.  mi.). 
The  normalized  discharge  can  be  regarded  as  the  discharge 
per  unit  area.  Figure  6.1  presents  recorded  maximum 
discharges  and  a  number  of  PMF  estimates  for  the  Red  Deer 
River  Basin. 

6.1  Recorded  Maximum  Discharges 

The  available  maximum  recorded  discharges  are  indicated 
by  triangles  in  Figure  6.1.  The  enveloping  curve,  a  curve  in 
which  all  discharge  values  are  equalled  or  exceeded,  is 
commonly  used  in  hydrological  studies  to  identify  a 
relationship  between  discharge  and  area.  This  curve  is 
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Figure  6.1.  Comparison  of  Recorded  Maximum  Discharge  and  PMF 
Estimates  for  the  Red  Deer  River  Basin. 
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obtained  graphically  and  may  be  either  a  curved  or  a 

straight  line.  For  simplicity,  the  latter  is  most  frequently 

used.  The  enveloping  curve  for  the  maximum  recorded 

discharges  in  Figure  6.1  can  be  approximated  by  a 

0.5  . 

relationship  of  the  form  Q  =  1000A  ;  this  straight  line 

coincides  with  the  100-year  flood  estimates  given  by  Neill 
(  1965)  . 

6.2  Muz i K' s  PMF  Estimate 

The  PMF  estimate  obtained  by  Muzik  (1975)  (from  Figure 
1.15)  at  Raven  in  the  Red  Deer  River  Basin,  is  designated  by 
the  letter  "X"  in  Figure  6.1.  This  estimate  was  obtained  by 
imposing  a  292  mm  (11.5  in.)  PMP  loading  distributed  over  a 
48-hour  period,  so  that  228  mm  (9  in. )  were  distributed  the 
first  day  and  the  balance  of  64  mm  (2.5  in.)  was  distributed 
the  second  day.  Since  no  mass  distribution  was  given  with 
this  estimate,  the  228  mm  precipitation  loading  was  probably 
uniformly  distributed  over  a  24-hour  period  (giving  an 
average  intensity  equivalent  to  9.65  mm  (0.38  in.)  per 
hour).  The  resulting  instantaneous  discharge  of  4  814  m3  per 
sec  (170  000  cfs)  obtained  by  Muzik  is  consistent  with 
results  obtained  by  the  author  using  the  SSARR  Model  (not 
presented  here)  for  a  similar  average- i ntensi ty  PMP  loading 
at  Raven  in  the  Red  Deer  River  Basin.  It  would  seem  from  the 
literature  (Muzik,  1975)  that  this  estimate  is  strictly  due 
to  a  PMP  loading  without  the  snowmelt  loading. 
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6.3  PMF  Estimates  from  TAPMF,  SSARR,  and  HYMO  Models 

PMF  estimates  for  a  number  of  watersheds  in  the  Red 
Deer  River  Basin  as  obtained  from  the  three  models  (TAPMF, 
SSARR,  and  HYMO)  are  also  presented  in  Figure  6.1.  For  a 
given  drainage  area,  all  three  models  produced  greater 
normalized  discharge  than  either  the  maximum  recorded 
discharge  or  the  PMF  discharge  estimate  by  Muzik  (1975).  The 
estimates  from  the  three  models  seem  to  be  comparable  in 
value.  Any  differences  can  be  attributed  to  the  data 
required  and  the  mathematical  approximations  used  in  the 
models.  For  example,  in  the  HYMO  Model  a  value  of  the 
average  slope  is  used  for  simulation,  but  this  parameter  is 
not  needed  in  the  other  two  models.  In  the  SSARR  Model, 
snowmelt  is  estimated  using  the  generalized  snowmelt 
equation  for  a  partly  forested  area,  while  in  the  HYMO  and 
TAPMF  models  the  snowmelt  is  simulated  by  adding  the  melt 
rate  (Chapter  3)  to  the  precipitation  mass  distribution. 
Another  example  of  a  difference  between  the  models  is  in  the 
technique  used  for  routing  of  floods.  In  the  HYMO  Model,  the 
variable  storage  coefficient  flood-routing  method  (Williams, 
1969)  is  used  for  routing  of  floods  through  streams  and 
valleys  (Section  5.3.0).  This  method  accounts  for  the 
variation  in  water  surface  slope  during  a  flood.  In  the 
SSARR  and  TAPMF  models,  routing  is  accomplished  by  solution 
of  the  basic  storage  equation  in  finite  time  periods  for 
each  increment  of  storage  routing  (Sections  4.2.6  and 
5.2.3) . 
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It  was  originally  expected  by  the  author  that  the 
results  from  the  SSARR  Model  would  be  vastly  different  from 
those  of  the  HYMO  Model,  since  the  former  model  is 
considered  by  many  users  to  be  a  more  superior  hydrological 
model.  The  comparison  in  Figure  6.1  seems  to  suggest  that, 
to  the  contrary,  either  of  the  three  models  can  be  used  for 
PMF  estimates.  It  was  also  believed  that  the  computations  of 
the  shape  of  the  watershed  would  significantly  influence  the 
PMF  estimate.  To  test  this  concept  the  TAPMF  Model  was 
developed.  Comparison  of  the  results  from  the  TAPMF  Model 
with  the  results  obtained  from  the  other  two  models  seems  to 
dismiss  this  belief.  Because  the  results  were  comparable,  it 
is  not  practical  (due  to  time  and  resource  constraints)  to 
estimate  the  PMF  for  the  other  watersheds  in  Alberta  by 
using  all  three  models.  Instead,  it  was  decided  to  use  the 
model  requiring  the  least  amount  of  data  and  with  minimal 
computational  costs  for  estimating  the  PMF  for  all  other 
watersheds  in  Alberta.  The  model  that  best  met  these 
requirements  was  the  HYMO  Model ;  hence  this  model  was  used 
for  other  simulations  in  this  thesis.  The  HYMO  Model 
requires  less  data,  or  at  least  data  that  are  more  easily 
accessible,  for  simulation  than  either  the  SSARR  and  the 
TAPMF  models.  The  HYMO  Model  also  costs  approximately  the 
same  as  the  TAPMF  Model  and  less  than  the  SSARR  Model  for 


computer  simulation. 
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6.4  Ratio  of  PMF  to  Maximum  Recorded  Discharge 

A  second  comparison  of  the  three  models'  results 
involved  examining  the  ratio  of  the  PMF  estimate  (Or)  to  the 

maximum  recorded  discharge  (Q  )  for  each  watershed.  This 

MR 

ratio  gives  the  number  of  times  the  PMF  estimate  is  greater 

than  the  recorded  discharge.  Figure  6.2  shows  that  the  PMF 

* 

estimates  from  the  HYMO  Model  are  from  over  10  to  several 
100  times  the  recorded  values.  The  recorded  discharges  are 
the  largest  available  values,  and  for  many  stations  these 
values  are  small  because  large  discharges  have  not  yet  been 
recorded  due  to  the  short  record  period. 

For  completeness,  the  PMF  estimate  from  the  HYMO  Model 
was  compared  to  recorded  flows  using  a  frequency  analysis 
for  the  Red  Deer  River  at  the  Red  Deer  station  (Figure  6.3). 
This  station  was  chosen  because  of  its  long  record  length 
compared  to  the  other  stations  in  the  basin  and  also  because 
this  station  had  the  largest  recorded  discharge  in  the  Red 
Deer  River  Basin.  A  frequency  analysis  is  usually  displayed 
on  Gumbel  extreme-va 1 ue  logarithmic  probability  paper,  where 
discharge  data  from  a  small  period  (say,  50  years  or  so)  is 
extrapolated  to  return  periods  of  thousands  or  millions  of 
years.  The  plotting  positions  on  Figure  6.3  were  determined 
by  ranking  the  extreme  yearly  values  of  the  discharge,  from 
the  highest  to  the  lowest  (m  =  1,2,...,n)  and  calculating 
the  probability  (e.g.,  P  =  m  /  (n  +  D).  In  Figure  6.3, 
extrapolating  a  straight  line  through  the  points  for  the  Red 
Deer  station  suggests  that  the  PMF  estimate  of  21  041  m3  per 
second  (743  482  cfs)  would  have  a  return  period  of  9  x  106 
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Figure  6.2.  Ratio  of  the  PMF  to  Recorded  Discharge  for 
Watersheds  in  the  Red  Deer  River  Basin. 
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years.  As  a  general  rule,  however,  a  frequency  analysis  is 
not  recommended  by  statisticians  when  the  return  periods  to 
be  estimated  are  greater  than  twice  the  record  length 
( Pugsley ,  1981 ) . 
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Figure  6.3.  Frequency  Analysis  of  the  Annual  Maxima 
Discharges  of  the  Red  Deer  River  at  Red  Deer. 
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CHAPTER  7 


PMF  EMPIRICAL  RELATIONSHIPS  FOR  ALBERTA  RIVER  BASINS 

7.1.0  Introduction 

For  hydrological  applications  a  relationship  is  desired 
to  relate  peak  PMF  discharge  to  other  parameters  that  can  be 
easily  measured  or  computed.  A  parameter  that  has  been 
considered  by  many  authors  (numerous  examples  are  given  in 
Gray,  1970)  is  the  drainage  area.  An  equation  using  this 
parameter  is  called  a  Q-A  relationship.  The  general  form  of 
a  Q-A  relationship  is 

Am 

Q  =  cx  A  (7.1) 

where  Q  is  the  peak  discharge,  A  is  the  drainage  area  of  the 
watershed,  cl  is  a  coefficient  believed  to  be  dependent  on 
various  meteorological,  hydrological,  geological,  and 
watershed  parameters,  and  m  is  a  constant  exponent  that 
governs  the  slope  of  the  relationship.  The  exponent  can  be 
either  a  positive  or  a  negative  value,  thus  allowing  an 
infinite  number  of  combinations  relating  the  discharge  to 
area.  Many  different  combinations  have  been  suggested  by 
various  authors,  and  the  two  that  are  frequently  advocated 
are  the  relationships  with  m  =  0.5  and  m  =  1.0.  These  two 
forms  of  Equation  7.1  are  examined  using  the  HYMO  Model  for 
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gauged  watersheds  in  the  five  major  river  basins  in  Alberta 
to  define  the  coefficient  c^  for  estimating  the  PMF.  Once  a 
coefficient  is  defined  for  gauged  watersheds  in  a  basin,  it 
may  be  possible  to  use  this  value  to  obtain  PMF  estimates 
for  ungauged  watersheds  in  the  same  basin. 

The  work  discussed  in  this  chapter  is  divided  into 
three  sections.  In  the  first  section  (7.2.0)  the  square-root 
approximat ion  (or  m  =  0.5  in  Equation  7.1)  is  examined.  Some 
authors  (e.g.,  McKay  and  Stichling,  1961)  have  suggested 
that  this  approximation  defines  the  enveloping  curve  of  the 
Q- A  relationship.  A  multilinear  regession  analysis  of  peak 
discharge  and  area  gave  a  value  of  m  =  0.51  for  Alberta. 

This  value  suggests  that  a  square-root  approximation  may  be 
a  valid  consideration  to  define  the  enveloping  curve  in  the 
Q-A  relationship. 

In  the  second  section  (7.3.0)  the  second  form  of 
Equation  (7.1)  (called  the  Rational  Formula)  is 
investigated.  The  Rational  Formula  has  been  used  for 
estimates  of  peak  discharge.  The  formula  is  obtained  by 
setting  c^  =  KCI  and  m  =  1.0  in  Equation  7.1  so  that 

Q  =  K  C  I  A  (72) 

where  Q  is  the  peak  discharge,  K  is  the  conversion  constant 
(equivalent  to  0.278  if  Q  is  in  m3  per  second  and  645.3  if  0 
is  in  cfs),  C  is  the  runoff  coefficient,  and  I  is  the 
average  rainfall  intensity  lasting  for  period  of  time  t 
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(time  of  concentration).  Equation  7.2  is  investigated 
(Section  7.3.0)  to  define  the  spatial  distribution  of  the 
runoff  coefficient  for  PMF  estimates  in  Alberta  river 
bas i ns . 

In  the  third  section  the  normalized  discharge  and  the 
runoff  coefficient  are  examined.  Because  these  vary  within  a 
basin  and  from  basin  to  basin  a  further  objective  of  this 
study  was  to  investigate  the  spatial  distribution  of  the 
normalized  discharge  and  the  runoff  coefficient  for  the 
river  basins  in  Alberta.  This  work  is  presented  in  Section 
7.4.0. 

7.2.0  PMF  Relationship  (m  =  0.5) 

The  coefficent  c^  is  believed  by  many  authors  (e.g., 

McKay  and  Stichling,  1961)  to  be  a  function  of  many 

variables,  including  the  climate,  geology,  and  hydrology  of 

a  given  watershed.  For  PMF  estimates,  because  of  the 

requirement  for  saturated  soil  conditions,  a  number  of  these 

variables  become  less  important  and  contribute  little  to  the 

estimation  of  the  coefficient.  Slope  and  the  intensity  of 

the  precipitation  loading  were  identified  as  important 

parameters  affecting  the  PMF  estimates.  Parameters  such  as 

maximum  watershed  width  and  length  are  to  some  degree 

0.5 

covered  by  the  A  approximation  and  hence  are  not  expected 
to  significantly  influence  the  coefficient. 

Three  types  of  intensities  can  be  identified:  (1)  the 
intensity  for  precipitation  loading,  (2)  the  intensity  for 
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snowmelt  loading,  and  (3)  the  sum  of  the  precipitation  and 
snowmelt  loading,  i.e.,  the  intensity  for  water  loading.  Of 
these  three,  the  intensity  of  the  precipitation  loading  was 
selected  for  analysis  in  this  work  for  simplicity.  Either  of 
the  other  two  intensities  could  also  have  been  used. 

7.2.1  Relationship  of  Coefficient  c^  with  Slope 

A  watershed  characteristic  that  can  be  easily  estimated 
is  the  average  slope.  The  average  slope  of  watersheds  in  the 
Red  Deer  River  Basin  was  plotted  against  the  coefficient  c 
(Figure  7.1)  to  investigate  any  possible  relationship 
between  the  two  parameters.  The  coefficient  was  calculated 
from  Equation  7.1,  using  the  PMF  discharge  (data  given  in 
Appendix  II)  from  the  HYMO  Model  and  the  drainage  area  for 
each  watershed.  The  results  of  these  computations,  as 
displayed  in  Figure  7.1,  indicate  a  logarithmic  increase  in 
the  coefficient  with  an  increase  in  slope.  These  results 
further  suggest  that  the  slope  may  contribute  to  the 

estimation  of  coefficient  c  . 

1 

7.2.2  Relationship  of  Coefficient  c_^  with  Intensity 

The  second  parameter  believed  to  affect  the  coefficient 
c^  is  the  intensity  of  the  precipitation  loading.  To 
investigate  this  aspect,  three  different  simulations  using 
the  HYMO  Model  were  carried  out  for  the  Sundre  watershed  of 
the  Red  Deer  River  Basin.  In  these  simulations  all 
parameters  were  kept  constant,  and  only  the  intensity  of  the 
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Figure  7.1.  Relationship  between  Slope  and  Coefficient  c1 . 
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precipitation  loading  was  varied  from  simulation  to 
simulation.  The  constant  parameters  were  (a)  the  drainage 
area  (2471  Km2  or  954  sq .  mi.),  (b)  the  average  slope  (5.7  m 
per  Km  or  30.1  ft.  per  mi.),  and  (c)  24-hour  total  water 
loading  (645  mm  or  25.4  in.).  The  results  of  these  three 
simulations  are  depicted  in  Figure  7.2. 

The  first  simulation  is  depicted  by  Curve  A,  the 
resultant  hydrograph  of  a  water  loading  intensity  of  43.2  mm 
(1.70  in.)  per  hour.  This  intensity  is  composed  of  a  12.7  mm 
(0.5  in.)  per  hour  (for  24  hours)  snowmelt  loading  and  a 
30.5  mm  (1.20  in.)  per  hour  (for  9  hours)  precipitation 
loading.  The  resulting  peak  discharge  of  21  055  m3  per 
second  (743  482  cfs)  giving  c1  equal  to  24  071. 

For  the  second  simulation,  Curve  B  depicts  the 
hydrograph  for  the  same  watershed  but  with  a  water  loading 
intensity  of  35.6  mm  (1.40  in.)  per  hour,  composed  of  a  12.7 
mm  (0.5  in.)  per  hour  (for  24  hours)  snowmelt  loading  and  a 
22.9  mm  (0.90  in.)  per  hour  (for  9  hours)  precipitation 
loading.  A  peak  discharge  of  19  829  m3  per  second  (700  183 
cfs)  and  the  coefficient  c1  equal  to  22  669  were  obtained. 

The  third  simulation  is  shown  by  Curve  C  for  a  water 
loading  intensity  of  26.9  mm  (1.06  in.)  per  hour,  composed 
of  a  12.7  mm  (0.5  in.)  per  hour  (for  24  hours)  snowmelt 
loading  and  a  14.2  mm  (0.56  in.)  per  hour  (for  24  hours) 
precipitation  loading.  A  peak  discharge  of  16  735  m3  per 
second  (590  910  cfs)  and  a  coefficient  value  of  19  131  were 
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Figure  7.2. 
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riation  of  the  PMF  Discharge  with  Changes 
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The  three  curves  in  Figure  7.2  show  that  the  PMF  peak 
discharge  increases  with  an  increase  in  the  intensity  of  the 
precipitation  loading  (even  though  the  total  water  loading 
remained  constant  for  the  24-hour  period).  The  increase  is 
suspected  to  be  logarithmic.  To  verify  this,  the 
relationship  between  the  coefficient  c^  and  the 
precipitation  loading  intensity  was  plotted  (Figure  7.3). 
Three  different  average  slopes  were  simulated,  since  it  was 
suspected  that  slope  might  also  affect  the  results.  Figure 
7.3  confirms  that  the  coefficient  c  is  logarithmically 
related  to  (a)  the  intensity  of  the  precipitation  loading 
and  (b)  the  average  slope.  Of  these  two  quantities,  the 
precipitation  loading  intensity  is  the  dominant  contributor 
to  the  coefficient  c^,  since  increases  in  intensity  produce 
larger  increases  in  the  coefficient  than  similar  increases 
for  the  slope.  Because  the  precipitation  loading  intensity 
makes  up  the  water  loading  intensity,  the  above  statments 
can  also  be  applied  to  water  loading.  Varying  the  snowmelt 
loading  intensity  will  also  have  an  effect  (probably  similar 
to  that  observed  for  the  precipitation  loading  intensity), 
but  this  was  not  examined  in  this  work. 

7.3.0  PMF  Relationship  (m  =  1.0) 

The  second  Q-A  relationship  investigated  is  Equation 
7.2,  or  Equation  7.1  with  m  =  1.0  and  c^  =  KCI.  The 
coefficient  C  (called  the  runoff  coefficient)  is  divided 
into  three  parts:  (1)  K,  the  conversion  constant,  (2)  C,  the 
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7.3.  Relationship  of  the  Water  Loading  Intensity  with 
the  Coefficient  c1  . 
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runoff  coefficient,  and  (3)  I,  the  water  loading  intensity. 
In  Equation  7.2  intensity  is  considered  a  variable  parameter 
and  not  part  of  the  coefficient  c^»  as  was  the  case  in 
Section  7.2.0.  The  water  loading  intensity  is  considered 
here  since  the  runoff  coefficient  is  defined  as  the  ratio 
between  surface  runoff  following  a  rainfall  on  a  snowpacK 
and  the  total  volume.  The  coefficient  C  was  examined  in 
relation  to  the  average  slope  using  the  HYMO  Model,  and  the 
resulting  scatter  for  the  Red  Deer  River  Basin  is  depicted 
in  Figure  7.4.  The  results  in  the  figure  indicate  a 
logarithmic  relationship  between  the  two  parameters.  The 
graphs  for  the  other  Alberta  river  basins  are  given  in 
Appendix  IV.  These  graphs  also  indicate  a  logarithmic 
re  1  at i onshi p . 

7.4.0  PMF  Estimates  for  Alberta  River  Basins 

The  normalized  discharge  was  thought  to  vary  within  a 

basin;  therefore,  this  aspect  was  investigated  by  using  the 

HYMO  Model  with  the  PMF  estimates  obtained  for  the  various 

Alberta  river  basins  depicted  in  Figure  7.5.  Also  given  in 

this  figure  are  the  maximum  recorded  discharges  for 

available  watersheds.  The  enveloping  curve  for  the  maximum 

recorded  discharges  seems  to  be  well  approximated  by  the 

0  5 

equation  Q  =  2000A  .  The  enveloping  curve  for  the  PMF 

estimates  can  be  approximated  by  an  equation  Q  =  1.35  x  105 
0  29 

A  f  suggesting  a  relationship  with  a  smaller  exponent 
value  (m  =  0.29)  than  the  square-root  value  obtained  for  the 
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Figure  7.4.  Variation  of  the  Runoff  Coefficient  C  with 
Average  Slope  for  the  Red  Deer  River  Basin. 
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Figure  7.5.  Relationship  Between  Normalized  Discharge  and 
Drainage  Area  for  Alberta  River  Basins. 
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recorded  discharges.  An  equation  with  a  square-root 
approximation  (i.e.  Q  =  3.16  x  104  A  is  also  drawn  on 
Figure  7.5  to  show  the  deviation  of  this  equation  from  the 
one  with  an  exponent  of  m  =  0.29.  The  separation  of  the 
normalized  discharge  by  basin  is  presented  in  APPENDIX  IV, 
together  with  the  maximum  recorded  discharges  for  the 
correspondi ng  drainage  areas.  For  individual  river  basins 
(APPENDIX  IV),  the  results  seem  to  be  represented  by  an 
enveloping  curve  of  the  form  Q  c<  f(A  ).  The  large  spread 
in  the  values  of  the  normalized  discharge  is  believed  to  be 
due  to  the  spatial  variation  of  the  average  slope.  Also 
presented  in  APPENDIX  IV  is  the  graphical  representat i on  of 
the  relationships  between  slope  and  PMF  discharge  and  slope 
and  runoff  coefficient  for  the  different  basins. 

7.4.1  Spatial  Distribution  of  the  Normalized  Discharge 
The  normalized  discharges  were  examined  in  terms  of 
their  spatial  distribution  in  Alberta.  The  normalized 
discharge  can  be  computed  by  dividing  the  peak  PMF  discharge 
by  the  drainage  area  from  the  data  provided  in  APPENDIX  III. 
The  computer  contoured  results  of  these  computations  are 
depicted  in  Figure  7.6.  The  computer  contouring  system  used 
in  this  analysis  was  developed  for  the  Kansas  Geological 
Survey  and  is  presented  by  Sampson  (1975)  in  the  user's 
manual,  "SURFACE  II  Graphics  System".  SURFACE  II  is  a 
computer  software  system  for  creation  of  displays  of 
spatially  distributed  data.  The  form  of  graphic  display 
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Figure  7.6.  Spatial  Distribution  of  Normalized  Discharge  for 

Alberta  River  Basins. 
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shown  in  this  work  (Figure  7.6)  is  a  contour  map,  a  plot  of 
two  coordinates  (values  of  spatial  coordinates  x  and  y)  on 
which  values  of  the  third  variable  (variable  to  be 
contoured)  are  defined  by  lines  of  equal  values.  In  Figure 
7.6,  the  contoured  lines  depict  the  normalized  discharge.  A 
large  number  of  computer  commands  are  available  in  the 
SURFACE  II  graphic  system.  The  user  first  generates  a  grid 
matrix  of  normalized  discharge  values  from  irregularly 
spaced  data  points  by  using  a  distance  weighting  function  in 
the  averaging  process.  Sample  data  points  used  in  the 
estimation  procedure  are  weighted  so  their  influence 
declines  with  distance  from  the  point  being  estimated. 
Following  this,  a  contour  map  from  the  grid  matrix  is 
generated  by  the  graphic  system.  Smoothing  is  done  by 
piecewise  Bessel  intepolation  (a  spline  fitting  method) 
(Sampson,  1975)  within  the  grid  cell,  which  results  in  a 
smooth  path  of  contour  line  through  the  grid. 

The  results  in  Figure  7.6  show  that  watersheds  in  the 
foothills  and  along  the  continental  divide  generally  have 
high  normalized  discharge  values  ranging  from  several 
hundred  to  slightly  over  1000  cfs  per  sq .  mi..  The  high 
values  coincide  to  watershed  with  steep  average  slopes. 
Smaller  (values  less  than  200  cfs  per  sq .  mi.)  of  the 
normalized  discharge  were  obtained  for  the  northern  and 
eastern  sections  of  the  province. 
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7.4.2  Spatial  Distribution  of  the  Coefficients 

Another  aspect  examined  was  the  spatial  distribution  of 
the  coefficients  c ^  and  C  using  the  SURFACE  II  Graphic 
System.  These  results  are  depicted  in  Figures  7.7  and  7.8., 
respectively.  Generally,  high  coefficient  values  are  found 
for  the  foothills  and  along  the  continental  divide 
coinciding  with  watersheds  of  high  average  slopes.  Low 
coefficient  values  were  obtained  for  watersheds  in  the 
eastern  and  northern  sections  of  the  province,  coinciding 
with  watersheds  with  low  average  slopes. 

7.4.3  Topic  for  Future  Research 

The  results  obtained  in  Section  7.2.1  suggest  that  the 
average  slope  may  contribute  to  the  estimation  of  the  peak 
discharge  and  needs  to  be  considered  in  a  Q-A  relationship. 
To  examine  this  aspect  a  preliminary  investigation  was 
conducted  using  a  multilinear  regression  analysis  with  the 
three  parameters  (peak  discharge,  average  slope,  and 
drainage  area  as  given  in  APPENDIX  A. 3)  expressed  in 
logarithmic  form.  Using  this  analysis  the  following  equation 
was  obtained  for  Alberta  with  a  mu  1 1 i cor re  1  at i on  coefficient 
of  0.97: 

Q-  505  A0’78  S°-5?  <7'3) 

where  Q  is  the  peak  discharge  (cfs),  A  is  the  drainage  area 
(sq.  mi.),  and  S  is  the  average  slope  (feet  per  mile). 
Equation  7.3  was  computed  for  average  slopes  of  1.0,  10.0, 
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Figure  7.7.  Spatial  Distribution  of  the  Coefficient  for 

Alberta  River  Basins. 
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Figure  7.8.  Spatial  Distribution  of  the  Runoff  Coefficient  C 

for  Alberta  River  Basins. 
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and  100.0  feet  per  mile,  and  these  computations  are  depicted 
in  Figure  7.9  by  the  solid  lines.  The  data  used  in  the 
analysis  is  also  shown  in  Figure  7.9,  with  the  position  of 
the  various  values  represented  by  the  dots.  The  numerical 
values  next  to  the  dots  give  the  average  slope  for  a  given 
discharge  and  area.  Equation  7.3  can  be  considered  as  a 
preliminary  analysis  on  this  topic  and  further  work  is 
suggested . 


Discharge  (cfs) 


Area  (sq.  mi.) 


Figure  7.9.  Relationship  between  Peak  Discharge,  Average 
Slope,  and  Drainage  Area  for  Alberta  River  Basins. 
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CHAPTER  8 

DISCUSSION  AND  CONCLUSIONS 


8.1.0  Di scussion 

A  number  of  topics  dealing  with  estimation  of  PMP  and 
PMF  were  examined  by  the  author.  A  comprehensive 
climatological  analysis  of  rainstorms  was  conducted  using  a 
novel  approach  in  which  the  isohyetal  patterns  and  rainfall 
intensities  from  rainstorms  were  digitized  in  spatial  blocks 
and  analyzed.  Such  an  analysis  is  the  first  of  its  Kind  for 
Alberta,  and  is  lacking  for  many  provinces  in  Canada. 

PMP  was  estimated  using  two  common  approaches,  the 
meteorological  and  statistical  techniques,  and  included  an 
extensive  analysis  and  graphical  representat i on  of  the  (a) 
dewpoint  temperature,  and  (b)  enveloping  PMP  curves  for 
Alberta  river  basins.  Both  techniques  gave  similar  estimates 
suggesting  that  either  approach  can  be  used  for  Alberta 
river  basins. 

PMP  estimates  were  also  used  to  obtain  PMF  estimates.  A 
unique  concept  employed  the  spatial  relationship  of  PMP 
estimates  to  determine  the  location  of  the  largest  PMP 
loading  in  a  basin.  In  the  study,  the  author  found  that  such 
a  relationship  was  well  correlated  within  a  river  basin  and 
mainly  dependent  on  position  and  elevation  of  the  location 
for  which  the  estimate  is  required. 
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The  rain-on-snow  event  is  believed  to  produce  the 
largest  water  loading  on  watersheds  in  Alberta.  The 
important  processes  contributing  to  snowmelt  were  defined, 
and  an  equation  was  developed  to  estimate  maximum  snowmelt. 
Previously,  convection  and  condensation  were  identified  as 
necessary  to  the  total  melt  in  the  rain-on-snow  event; 
however,  analysis  showed  that  shortwave  radiation  and 
rainfall  are  also  important. 

To  investigate  the  effect  of  watershed  shape  on  PMF , 
the  TAPMF  Model  was  developed.  PMF  estimates  obtained  from 
this  model  were  comparable  to  those  obtained  by  the  SSARR 
and  HYMO  models.  The  lack  of  a  significant  difference  in  the 
estimate,  attibutable  to  the  watershed  shape,  is  probably 
due  to  the  dominance  of  area  on  the  shape  parameters.  The 
sensitivity  of  the  SSARR  and  HYMO  models  was  also 
investigated  hypothetically  and  the  results  show  the 
importance  of  accurately  determining  the  required  model 
parameters.  All  three  models  gave  about  three  times  larger 
than  a  previous  PMF  estimate  in  the  Red  Deer  River  Basin; 
this  is  mainly  due  to  the  difference  in  the  precipitation 
intensity  loading. 

Two  empirical  Q-A  relationships  were  investigated  with 
exponents  equal  to  0.5  (a  square-root  approximation)  and  1.0 
(Rational  Forumla).  The  Rational  Formula  was  selected 
because  the  coefficient  is  dimensionless  and  hence  easily 
i nterpretable .  Also,  the  Q-A  relationship  obtained  by  other 
investigators  have  exponent  values  less  than  1.0  and  hence 
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the  exponent  equal  to  1.0  can  be  considered  an  upper  limit 
in  the  relationship. 

The  spatial  distributions  of  the  coefficients  in  the 
two  Q-A  relationships  and  the  normalized  discharge  are  also 
presented.  The  Q-A  relationship  with  an  exponent  equal  to 
0.29  best  represents  the  data  for  all  the  river  basins  in 
Alberta,  while  for  many  individual  river  basins  a  square- 
root  relationship  best  describes  the  enveloping  curve  of 
PMF.  The  square-root  relationship  seems  to  vary 
logar thmi ca 1 1 y  with  slope  and  intensity. 

This  study  presents  PMP  and  PMF  estimates,  previously 
lacking,  for  Alberta  river  basins. 

A  number  of  conclusions  were  reached  during  the  course 
of  this  study  pertaining  to  the  climatological, 
meteorological,  PMP,  and  PMF  estimates;  these  are  summarized 
below. 

8.2.0  Conclusions 

8.2.1  Climatology  and  Meteorology 

1.  Rainstorms  with  depths  of  50  mm  (2  in.)  and  more  are 
frequent  in  Alberta  (average  of  11  occurrences  per  year), 
while  rainstorms  with  depths  of  150  mm  (6  in.)  and  more  are 
uncommon  (average  of  an  occurrence  per  two  year). 

2.  The  majority  of  rainstorms  occur  in  June  and  July,  with 
only  a  small  percentage  occurring  in  April  and  September. 
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3.  The  greatest  frequency  of  rainstorms  is  in  the 
mountainous  area  in  the  southern  part  of  the  province  (i.e., 
Waterton  Lakes  National  Park)  with  an  average  of  two 
rainstorms  a  year  with  depths  50  mm  (2  in.)  and  more.  In 
this  area  a  rainstorm  with  a  depth  of  100  mm  (4  in.)  or  more 
occurs  on  the  average  once  every  3  years,  while  one  with  a 
depth  of  150  mm  (6  in.)  or  more,  occurs  once  every  10  years. 

4.  The  frequency  of  rainstorms  decreases  eastward  and 
northward  from  the  continental  divide  in  Alberta. 

5.  Severe  rainstorms  (rainstorms  with  depths  of  150  mm  (6 
in.)  and  more)  occur  in  four  main  regions  (or  belts)  of  the 
province:  the  first  extends  through  southern  Alberta,  just 
south  of  Calgary;  the  second  is  in  central  Alberta  from 
south  of  Edson  eastward  to  the  Edmonton  region;  the  third  is 
from  Lesser  Slave  Lake  nor theastward  to  Fort  McMurray;  and 
the  fourth  is  around  the  Fort  Vermilion  area. 

6.  In  Alberta,  heavy  precipitation  events  (depths  of  150  mm 
(6  in.)  or  more)  are  produced  by  a  cold  low  system. 

8.2.2  PMP  and  PMF  Estimates 

Several  conclusions  can  be  drawn  from  the  analyses  of 
the  depth- area-dur at  ion  curves  and  the  estimates  of  PMP  and 
PMF. 

1.  DAD  analyses  using  the  physical  approach  show  that  only  a 
few  rainstorms  contribute  to  the  estimate  of  the  PMP.  The 
storms  that  contribute  to  the  PMP  are  not  obvious,  hence  the 
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procedure  of  plotting  DAD  curves  is  an  important  and 
necessary  part . 

2.  The  maximum  estimates  of  the  PMP  occur  in  June  for  the 
basins  in  the  southern  portion  of  the  province,  while  in 
central  and  northern  Alberta  these  occur  in  July. 

3.  Estimates  of  PMP  vary  spatially  within  each  river  basin 
in  Alberta.  Generally,  two  characteristic  variations  can  be 
defined:  the  first  and  the  most  predominant  is  a  decrease 
west  to  east  in  the  basin,  while  the  second  is  a  decrease 
from  south  to  north. 

4.  For  many  of  the  river  basins  in  Alberta,  estimates  from 
the  statistical  technique  produced  results  similar  to  those 
from  the  physical  approach. 

5.  Maximum  water  loading  on  a  watershed  can  be  expected  in 
late  May  or  early  June  in  Alberta,  since  the  coincidence  of 
severe  rainstorms,  snow  cover,  and  air  temperatures  above 
freezing  seems  to  occur  during  this  period. 

6.  For  rain-on-snow  events  under  a  PMP  loading,  the 
predominant  heat  transfer  processes  are  due  to  convection 
and  condensation,  shortwave  and  longwave  radiation,  and 
rainfall.  Maximum  snowmelt  can  be  expressed  as  a  function  of 
meteorological  parameters  (temperature,  PMP,  and  wind 
speed),  which  were  observed  to  be  well  correlated  with  the 
location,  elevation,  duration,  and  rainstorm  area. 

7.  Sensitivity  of  model  response  to  parameter  changes  is  an 
important  feature  of  the  SSARR  and  HYMO  models,  since  some 
parameters  affect  the  shape  of  the  resulting  hydrograph  and 
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the  peak  PMF  discharge. 

8.  The  three  models  for  PMF  produced  comparable  results, 
indicating  that  any  one  can  be  used  effectively  for 
watersheds  in  Alberta. 

9.  The  Q-A  relationship  expressed  as  Q  =  c  A  defines  the 
enveloping  curves  for  the  maximum  recorded  flows  and  PMF 
estimates  for  watersheds  in  Alberta. 

10.  The  coefficients  c1  in  Q  =  c-^A  and  C  (the  runoff 
coefficient)  in  Q  =  KCIA  in  the  Q-A  relationship  are 
logarithmically  related  to  the  average  slope  of  the 
watersheds . 

11.  High  normalized  discharge  and  runoff -coefficient  values 
are  mainly  observed  in  the  foothills  and  along  the 
continental  divide  where  watersheds  exhibit  high  average 
slopes.  Lower  values  of  the  these  parameters  were  obtained 
for  watershed  in  the  eastern  and  northern  portions  of  the 
province  (coinciding  with  watersheds  of  lower  average 
slopes ) . 
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APPENDIX  A. 1 


This  appendix  contains  tables  used  in  estimating  the 
PMP  by  the  physical  approach. 

Table  A. 1.1  presents  values  of  precipitable  water  (mm) 
between  the  1000  mb  surface  and  various  pressure  levels  up 
to  300  mb  in  a  saturated  pseudoadi abat i c  atmosphere  as  a 
function  of  the  1000  mb  dew  point.  Table  A. 1.2  lists 
similiar  values  for  layers  between  the  1000  mb  surface, 
assumed  to  be  at  zero  elevation,  and  various  heights  up  to  8 
Km. 
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Table  A  .  1 . 1 .  Preci pi  table  water  (mm)  between  1  OOP  mb  surface 
and  i ndi cated  pressure  ( mb )  i n  a  saturated  pseudoadi abat ic 
a tmosphere  as  a  f unct i on  of  the  1 000  mb  dew  point  (  °C) 


mb 
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R 

10 

1? 

13 

14 

15 

16 

17 

1£L 

L9 

-2D 

21 

2? 

23 

24 

990 

0 
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1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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7 

7 

940 
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3 
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4 

5 

5 

6 

6 

7 

7 
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9 

10 

10 

1  1 

920 
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9 

10 

10 

1  1 
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13 

14 
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65 
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15 
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21 
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40 
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48 

52 

57 

62 

67 

71 

T ABLE  A . 1 . 2  Precipi table  water  (mm)  between  1 000  mb  surface 
and  indicated  height  (m)  above  that  surface  i n  a  saturated 
pseudoadi abat i c  atmosphere  as  a  funct ion  of  the  1 000  mb  dew 


point  X  Cj. 
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33 
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42 
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APPENDIX  A. 2 


The  following  pages  describe  a  number  of  relationships 
of  air  temperature,  wind  speed,  and  precipitation  and  the 
mu  1 1 i cor re  1  at i on  coefficients  associated  with  each 
relationship  for  river  basins  in  Alberta.  The  notations  used 
in  these  relationships  are  defined  as  follows; 


Tm 

i  s 

mean 

la 

i  s 

the 

Lo 

i  s 

the 

E 

i  s 

the 

p 

r  mp 

i  s 

the 

A 

i  s 

the 

D 

i  s 

the 

Um 

i  s 

the 
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I .  North  Saskatchewan  River  Basin 


1.  Relationships  of  Air  Temperatures 


1.1 

T 

m 

-  130. 3A  +  0.2503  La  -  0.6627  L0 

1.2 

T 

m 

-  1AA.87  -  1.26A7  La  -  0.003169  Lq  -  0.003A65  E 

1.3 

log  T 

3  m 

•  3.2707  4  O.A3616  log  La  -  1.0663  log  Lo 

1.A 

log  T 

3  m 

-  A. 0913  -  2.3768  log  La  4  1 . 3A76  log  Lo 

-  0.272A2  log  E 

2.  Relationships  of  Wind  speed 


2.1 

U 

m 

-  17.365  + 

0.8865  La  - 

0.A679 

Lo 

2.2 

U 

m 

-  72.A73  - 

5.7A1A  La  4 

2.A667 

L  -  0.01A0  E 

0 

2.3 

log 

U 

m 

-  -12.513 

4  12.06A  log 

LA  ‘  3' 

.5232  log  L0 

2. A 

log 

U 

m 

-  0.925A  * 

3A.005  log 

La  ♦  36, 

.010  log  Lo 

-  A.A61A  log 

3.  Relationships  of  2A-hr  PMP 

3.1  P  -  -51.265  -  3-7902  L.  +  2.5287  -  0.009  E 

mp  «  ° 

3.2  log  Pmp  «  -20.21  +  3-6685  LA  +  7-2609  L0 

3.3  log  P  -  -11.502  -  26.182  log  LA  +  32.876  log  Lo 

-  2.8908  log  E 

3.1j  P  -  7O.A96  -  0.96661  T  +  0.51  A^42  U 

mp  mm 

3.5  log  Pmp  -  1A.591  -  7.7558  log  Tm  +  0.61750  log  Um 

A.  Relationships  of  D-hr  PMP 

A.1  P  -  5.2039  -  7- 3106  L.  +  3.8323  L  -  0.01588  E 
mp  «  0 

4  0.09558  D  -  0.0001AA  A 

A. 2  log  P  -  -2.372  -  31-58A  log  L.  4  33-A55  log  L 

mp  m  w 

-  3.2866  log  E  4  O.A1885  log  D  -  0.071215  log  A 


(r-0.591) 

(r=0. 762) 
(r-0.936) 

(r-0.999) 

(r=0. 375) 
(r-0.703) 
(r-0.691) 

(r*0.800) 

(r«0.9AA) 

(r-0.793) 

(r-0.906) 
(r-0.920) 
(r*0. 920) 


(r-0.909) 

(r-0.9A5) 


•  ^ 
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| | .  South  Saskatchewan  River  Basin 


1.  Relationships  of  Air  Temperatures 

1.1  T  -  313-61  +  0.581*5  L.  -  2.1*25  L 

m  n 

12  T  -  169.8  -  0.8261  L.  -  0.1*081  L  -  0.001*069  E 
m  M 

1.3  log  Tm  «  9-1855  +  0.1*31*8  log  -  3-9388  log  Lo 
1.1*  log  T  -  3-921  0.6511  log  L.  -  0.1003  log  L 

-  0.2191*  lo 

2.  Relationships  of  Wind  Speed 

2.1  U  -  -221.72  -  1.897  L.  +  3-0739  L  "  0.001*62  E  (r-0.9^9) 

m  H 

2.2  log  U  -  -1*1.372  -  6.3183  log  L  +  28.088  log  Lq 

m 

-  1.21*56  log  E  (r*0.959) 


(r-0.902) 

(r-0.967) 

(r-0.894) 

E  (r-0.978) 


3.  Relationships  of  D-hr  PMP 

3.1  P  -  -77.295  -  0.7653  L.  +  1.115  L0  +  0.1373  0 

J  mp  M 

-  0.0001*1  A  (r-o.91  A) 

3  2  P  -  -1*9.068  -  0.711*1*  L.  +  0.8216  L  +  0.00076  E 

3  mp  M  u 

+  0.1373  0  -  0.0001*1  A  (r*0.915) 

3.3  log  P  -  -13-169  -  2.2101*  log  L.  +  8.6522  log  LA 

mp 

-  0.08197  log  E  +  0.5029  log  D  -  O.O676  log  A  (r-0.936) 

3.1*  log  Pmp  -  -11.161  -  2.1137  log  l-A  +  7-1*51*7  log  L0 

■*■  0.5029  log  D  -  0.0676  log  A  (r-0.935) 


III.  Athabasca  River  Basin 


1.  Relationships  of  Air  Temperatures 


1.1 

T 

-  81 

.112 

+  0.1099 

L. 

-  0.1635  L 

m 

A 

0 

1.2 

T 

m 

-  99 

GO 

a\ 

-  0.9638 

la 

+  0.2459  Lq 

-  0.0033 

1.3 

log 

T 

m 

-  2. 

2064 

+  0.0966 

log 

La  -  0.262 

log  Lq 

1.4 

log 

T 

m 

-  4. 

GO 

r— 

-  1.751 

log 

La  +  0.6004 
A 

log  Lo 

-  0.1 A98  log 

2.  Relationships  of  Wind  Speeds 

2.1  U  -  87.153  -  0.1793  L.  -  0.6213  L  +  0.00148  E 

m  mo 

2.2  log  Um  -  49.124  -  26.131  log  1_A  +  1.277  log  Lo 

-  1.616  log 

3.  Relationships  of  D-hr  PMP 


3-1 

P  “ 

mp 

117.57  -  1.0002  La  -  0.4849  lQ  +  0.146  0 

-  0.00037 

3-2 

P  - 

mp 

104.69  -  0.5232  La  -  0.6245  Lq  +  0.0013  E 

+  0.146  D  -  0.00037 

3.3 

log  P  ■ 
mp 

19.2  -  4.841  log  La  -  5.08  log  Lq 

+  0.597  log  D  -  0.0706  log 

3-4 

log  P  ■ 
3  mp 

36.72  -  22.365  log  La  +  3-687  log  Lq 

-  1.5362  log  E  +  0.597  log  0  -  0.0706  log 

(r«0.4l9) 

(r-0.733) 

(r*0.4l9) 

(r*0.784) 

(r-0.570) 

(r-0.702) 

(r-0.960) 

(r-0.961) 

(r-0.962) 

(r«0.968) 


. 
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Peace 

River 

Basins 

Relat 

ionships  of  Air  Temperatures 

1.1 

T 

m 

-  62.546  +  0.5792 

la 

-  0*222  L 

0 

1.2 

T 

m 

-  63.091  -  0.5347 

la 

+  0.3654  L  -  0.0038  E 

O 

1.3 

log  T 

3  m 

-  1.7907  +  0.478 

log 

La  -  0.381  log  lo 

1.4 

log  T 

3  m 

-  1.778  +  0.0558 

log 

Lft  -  0.0436  log  Lq 

-  0.039  log  E 

2.  Relationships  of  Wind  Speeds 

(none  available  because  of  insufficient  wind  data) 

3.  Relationships  of  D-hr  PMP 

3.1  p  -  -251.19  -  0.0134  L,  +  2.19  L  -  0.094  D 

mp  H  w 

-  0.00028  A 

32  P  -  -27*.29  +  0.643  L.  +  2.037  L  +  0.0022  E 

mp  mu 

+  0.094  D  -  0.00028  A 

3  3  log  P  -  -59.416  +  1.0882  log  L.  +  27-93  log  L 

mp  M 

♦  0.508  log  D  -  0.0687  log  A 

3.4  log  P  -  -63.818  +  3.8176  log  L  +  27-310  log  L 
mp  " 

+  0.2773  log  E  +  0.508  log  D  -  O.O687  log  A 

V.  Combined  Basins  in  Alberta 

(North  Saskatchewan.  South  Saskatchewan.  Red  Deer.  Athabasca, 

1.  Relationships  of  Air  Temperatures 

1.1  Tm  -  112.31  *  1-2464  LA  +  0.3902  Lo  -  0.0089  E 

1  2  log  T  -  6.818  +  6.192  log  LA  -  7-631  log  LQ 

m 

1.3  log  Tm  -  6.468  +  4.233  log  LA  *  5-391  log  Lq 

-  0.259  log  E 

Relationships  of  Wind  Speeds 

2.1  Um  -  38.802  -  1.3324  lA  +  0.433  Lo  -  0.0018  E 

2.2  log  Um  -  3-4247  -  6.255  log  LA  ♦  4.765  log  Lq 

-  0.3965  log  E 


(r*0.869) 

(r-0.950) 

(r-0.869) 

(r-0.889) 

(r-0.955) 

(r-0.958) 

(r-0.965) 

(r»0.966) 

and  Peace) 

(r-0.469) 

(r-0.325) 

(r-0.328) 

(r-0.664) 

(r«0.65l) 


2. 


APPENDIX  A. 3 


The  following  pages  contain  information  calculated  and 
used  in  the  HYMO  Model  for  estimation  of  PMF .  The  notations 
used  are  defined  as  follows: 

(a)  Kellerhals  reach  number  (Kellerhals  et  al.,  1972) 

(b)  Water  loading  (inches) 

(c)  Total  intensity  (inches  per  hour) 

(d)  Drainage  area  ( sq .  mi.)  (Water  Survey  of  Canada) 

(e)  Average  slope  (feet  per  mile)  (Kellerhals  et  al.,  1972) 

(f)  PMF  (cfs) 

(g)  Instantaneous  daily  recorded  discharge  (cfs)  (Water 
Survey  of  Canada ) 


230 


( a ) 

(b)  ( 

c ) 

(d) 

(e) 

(f ) 

(g) 
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i  R i ver 

Basin 

9 

22.00 

2.85 
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2.69 
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10 

22.00 
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766062. 

37500. 

1  1 

22.00 
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8 
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4310. 

27 
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2340. 
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8500. 

29 

24.76 
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3.64 

124976. 

6270. 

30 

24.76 
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3.06 

124803. 

8670. 
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24.76 
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742. 
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108740. 

17360. 
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15 
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35 
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2.48 
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13.20 
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145000. 
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2.48 
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1.85 
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204500. 

39 

27.72 

2.48 
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1  .21 

20468. 

2310. 

40 
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32.21 

501838. 

18000. 

41 

27.72 
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6.34 

447946. 

39100. 

42 
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269846. 

3860. 

44 
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24000. 
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27.72 
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38 

27.72 

2.48 

21300. 

1.00 

963260. 

120000. 
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River 
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(a ) 

(b)  (c) 

(d) 

(e) 

(f ) 

(g) 

60 

24.96  2.04 
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0.69 

61819. 

14100. 

61 

24.96  2.04 
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13.20 
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74 
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40.66 
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15200. 

77 
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31.15 

278040. 
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78 
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79 
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(a) 

(b) 

(c) 

(d) 

(e ) 

(f ) 

(g) 

Red  Deer  River  Basin 

48 

25.40 

1.70 

954. 

30.10 

743482. 

23100. 

49 

25.40 

1  .70 

4480. 

6.90 

1203303. 

68250. 

50 

25.40 

1  .70 

9580. 

2.00 

838583. 

43000. 

55 
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1  .70 
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4.60 

216936. 

7400. 

52 
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1.70 
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143673. 
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119 
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1.70 

291  . 

4.90 

114335. 

316. 

57 

25.40 
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1400. 

6.34 

465673. 

2130. 

54 
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1.70 

754. 

2.69 

160094. 

4810. 

53 
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924. 

11.09 

463708. 

5310. 

51 

25.40 

1  .70 

16800. 

1  .58 

1057918. 

46200. 

APPENDIX  A. 4 


Graphs  of  HYMO  Model  results  obtained  for  Alberta  river 
basins.  These  results  include:  (1)  comparison  of  PMF 
normalized  discharge  with  drainage  area,  (2)  comparison  of 
the  runoff  coefficient  (C)  with  average  slope,  (3) 
comparison  of  the  ratio  of  PMF  to  maximum  recorded  discharge 
with  drainage  area,  and  (4)  comparison  of  PMF  normalized 
discharge  with  average  slope. 
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Q/fl  CFS/SQ.  MI. 
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PMF  Estimate  (HYMO  Model) 


FIGURE  A. 4.1  Comparison  of  the  PMF  Normalized  Discharge  with 
Drainage  Area  for  Alberta  River  basins. 


Q/fl  CFS/SQ.  MI. 
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FIGURE  A. 4. 2  Comparison  of  the  PMF  Normalized  Discharge  with 
Drainage  Area  for  the  Peace  River  Basin. 


Q/fl  CFS/SQ.  MI. 
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DRRINR&E  RRER  (  SQ.  MI.) 


FIGURE  A. 4. 3  Comparison  of  the  PMF  Normalized  Discharge  with 
Drainage  Area  for  the  Athabasca  River  Basin. 


Q/R  CFS/SQ.  MI. 
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FIGURE  A. 4.4  Comparison  of  the  PMF  Normalized  Discharge  with 
Drainage  Area  for  the  North  SasK.  River  Basin. 


Q/n  CFS/SQ.  MI. 
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FIGURE  A. 4. 5  Comparison  of  the  PMF  Normalized  Discharge  with 
Drainage  Area  for  the  Red  Deer  River  Basin. 
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DRAINAGE  ARE  A  (  SQ.  MI.) 


FIGURE  A. 4. 6  Comparison  of  the  PMF  Normalized 
Drainage  Area  for  the  South  SasK.  River 


Di scharge  wi th 
Bas i n . 


Q/fl  CFS/SQ.  MI. 
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DRAINAGE  AREA  [  SQ.  MI.) 


FIGURE  A. 4. 7  Comparison  of  the  PMF  Normalized  Discharge  with 
Drainage  Area  for  the  Bow  River  Basin. 


Q/fl  CFS/SQ.  MI. 
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-Jfc  PKF  Estimate  (HYMO  Model) 
&  Recorded  Maximum  Discharge 


DRRINR&E  RRER  (  SQ.  MI.) 


5 


105 


FIGURE  A. 4. 8  Comparison  of  the  PMF  Normalized  Discharge  with 
Drainage  Area  for  the  Oldman  River  Basin. 
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FIGURE  A. 4. 9  Comparison  of  the  Runoff  Coefficient  with  Slope 

for  Alberta  River  basins. 
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FIGURE 


o 


A. 4. 10  Comparison  of  the  Runoff  Coefficient  with 
Average  Slope  for  the  Peace  River  Basin. 
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FIGURE  A. 4. 11  Comparison  of  the  Runoff  Coefficient  with 
Average  Slope  for  the  Athabasca  River  Basin. 
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FIGURE  A. 4. 12  Comparison  of  the  Runoff  Coefficient  with 
Average  Slope  for  the  North  SasK.  River  Basin. 
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FIGURE  A. 4. 13  Comparison  of  the  Runoff  Coefficient  with 
Average  Slope  for  the  Red  Deer  River  Basin. 
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FIGURE  A. 4. 14  Comparison  of  the  Runoff  Coefficient  with 
Average  Slope  for  the  South  SasK.  River  Basin. 
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FIGURE  A. 4. 15  Comparison  of  the  Runoff  Coefficient  with 
Average  Slope  for  the  Bow  River  Basin. 
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FIGURE  A. 4. 16  Comparison  of  the  Runoff  Coefficient  with 
Average  Slope  for  the  Oldman  River  Basin. 
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FIGURE  A. 4.  1 
wi  th  D 


Q  (PMF)  *  PMF  Discharge 
Q  (MR)  *  Maximum  Recorded  Discharge 
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DRRINRGE  RRER  (  SQ.  MI.) 


Ratio  of  PMF  to  Recorded  Maximum  Discharge 
ainage  Area  for  Alberta  River  basins. 
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q  (PMF)  *  PMF  Discharge 
n  (  tvtr  ^  i  Maximum  Recorded  Discharge 


DRR I NRGE  RRER  (  SQ.  MI.) 


FIGURE  A. 4. 18  Ratio  of  PMF  to  Recorded  Maximum  Discharge 
with  Drainage  Area  for  the  Peace  River  Basin. 
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Q  ( PMF)  t  PMF  Discharge 
Q  (MR)  *  Maximum  Recorded  Discharge 
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DRAINAGE  AREA  (  SQ.  MI.) 


FIGURE  A. 4. 19  Ratio  of  PMF  to  Recorded  Maximum  Discharge 
with  Drainage  Area  for  the  Athabasca  River  Basin. 
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FIGURE  A. 4. 20  Ratio  of  PMF  to  Recorded  Maximum  Discharge 
with  Drainage  Area  for  the  North  SasK.  River  Basin. 
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Q  (MR)  i  Maximum  Recorded  Discharge 
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FIGURE  A. 4. 21  Ratio  of  PMF  to  Recorded  Maximum  Discharge 
with  Drainage  Area  for  the  Red  Deer  River  Basin. 
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FIGURE  A. 4. 22  Ratio  of  PMF  to  Recorded  Maximum  Discharge 
with  Drainage  Area  for  the  South  SasK.  River  Basin. 
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FIGURE 


Q  (PMF)  :  PMF  Discharge 
Q  (MR)  *  Maximum  Recorded  Discharge 
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A. 4. 23  Ratio  of  PMF  to  Recorded  Maximum  Discharge 
with  Drainage  Area  for  the  Bow  River  Basin. 
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Q  (PMF)  :  PMF  Discharge 
Q  (MR)  t  Maximum  Recorded  Discharge 
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FIGURE  A. 4. 24  Ratio  of  PMF  to  Recorded  Maximum  Discharge 
with  Drainage  Area  for  the  Oldman  River  Basin. 
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FIGURE  A. 4. 25  Comparison  of  the  PMF  Normalized  Discharge 
with  Average  Slope  for  Alberta  River  basins. 
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FIGURE  A. 4. 26  Comparison  of  the  PMF  Normalized  Discharge 
with  Average  Slope  for  the  Peace  River  Basin. 
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FIGURE  A. 4. 27  Comparison  of  the  PMF  Normalized  Discharge 
with  Average  Slope  for  the  Athabasca  River  Basin. 
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FIGURE  A. 4. 28  Comparison  of  the  PMF  Normalized  Discharge 
with  Average  Slope  for  the  North  SasK.  River  Basin. 
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FIGURE  A. 4. 29  Comparison  of  the  PMF  Normalized  Discharge 
with  Average  Slope  for  the  Red  Deer  River  Basin. 
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FIGURE  A. 4. 30  Comparison  of  the  PMF  Normalized  Discharge 
with  Average  Slope  for  the  South  SasK.  River  Basin. 
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FIGURE  A. 4. 31  Comparison  of  the  PMF  Normalized  Discharge 
with  Average  Slope  for  the  Bow  River  Basin. 
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FIGURE  A. 4. 32  Comparison  of  the  PMF  Normalized  Discharge 
with  Average  Slope  for  the  Oldman  River  Basin. 
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